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I.  EXECUTIVE  SUMMAKY 


pr^graa  la  nongerned  _w1th  ■  *•«*  major  problems  In  hydroacoustics; 


Thia. 

nameljC^  the  nature  of  the  low  wavenumber  characteristics  of  turbulent 
boundary  layer  pressure  fluctuations  including  their  relationship  with 
structural  vibration,  and  cavitation  noise  from  tip  vortices.  There  are 
four  major  tasks: 


a)  Development  of  signal  analysis  techniques  that  permit 
accurate  estimates  of  the  low  wavenumber  portion  of  the 
wavenumber^-frequency  pressure  spectrum. 

r 

b)  Measurements  of  the  pressure  spectrum  and  comparison  of 
experiment  with  a  theoretical  wave  guide  model. 

c)  Determination  of  the  structural  response  to  wall  layer 
pressure  fluctuations  and  the  influence  of  plate  vibration 
on  the  structure  of  turbulence  within  the  boundary  layer. 

d)  Measurement  of  the  spectral  characteristics  and  source 
distribution  of  cavitation  noise  in  tip  vortices  in  steady 
and  non-'steady  flow.  4r 


Both  computer  simulation  and  physical  measurements  are  involved.  All 
experimental  work  is  in  water,  utilizing  hydrodynamic  research  facilities 
at  the  St.  Anthony  Falls  Hydraulic  Laboratory. 

The  theoretical  and  algorithmic  aspects  of  our  research  on  the 
estimation  of  the  frequency-wavenumber  spectra  of  turbulent  boundary  layers 
(TBL)  proceeded  as  follows.  First,  a  statistical  model  of  the 
cross-spectral  density  matrices  (CSDM)  of  TBL  random  fields,  as  proposed  by 
Chase  and  Corcos  were  developed.  These  models  explicitly  included 
dependence  on  the  array  geometry,  measurement  of  time  and  signal- to-noise 
ratios  at  each  sensor  element,  and  formed  the  backbone  of  the  simulations 
that  were  carried  out  to  assess  the  performance  of  various  estimators. 
Next,  the  most  promising  "modern"  spectral  estimation  techniques  were 
evaluated  for  suitability  of  representing  the  TBL  spectral  models. 
Finally,  new  spectral  estimation  techniques  were  developed  that  were  aimed 
at  high-fidelity,  rather  than  high-resolution  spectral  estimation  of  the 
type  of  spectra  expected  from  the  TBL  randoa  fields.  Where  possible, 
theoretical  analyses  of  the  performance  of  the  estimators  were  made.  In 
all  cases,  extensive  simulations  were  carried  out  to  evaluate  both  the 
spectral  fidelity  and  variability  of  the  estimators. 

Capon's  Maximum  Likelihood  (ML)  prcved  to  be  the  most  stable  known 
modern  spectral  estimation  technique  'or  the  estimation  of  the  TBL 
wavenumber  spectra.  As  a  means  of  increasing  :he  effective  aperture  of  an 
array,  with  fixed  number  of  sensors,  ncnunlftrm  array  alternatives  were 
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Investigated  in  conjunction  with  the  MLM  estimator.  It  was  found  that 
array  geometries,  such  as  Minimum-Redundancy,  increase  resolution  in  the 
case  of  multiple  plane  waves,  but  do  not  offer  higher  fidelity  in  esti¬ 
mating  smooth  spectra.  This  was  attributed  to  the  increase  in  the  average 
level  of  the  array  sidelobes  when  a  narrower  mainlobe  is  generated  by  the 
extended  array  aperture.  Certain  limitations  of  this  method  were  iden¬ 
tified.  To  counter  these  deficiencies,  a  new  method,  denoted  the 
Covariance  Constrained  MLM  (CCMLM)  was  developed  to  alleviate  this  short¬ 
coming  and  in  the  process  improve  the  fidelity  of  the  MLM.  The  CCMLM  uses 
the  MLM  in  an  Iterative  fashion  to  generate  an  estimate  of  the  spectrum 
that  is  asymptotically  (large  number  of  iterations)  matched  to  the  measured 
CSDM.  A  practical,  stable  iterative  algorithm  was  designed  and  the  perfor¬ 
mance  of  the  method  was  evaluated  both  theoretically  and  through  extensive 
numerical  simulations.  The  CCMLM,  which  Includes  the  MLM  as  a  special 
case,  proved  to  be  the  best  technique  that  we  evaluated  or  developed  for 
the  estimation  of  spectra  such  as  those  suggested  by  Corcos  and  Chase. 

Research  on  boundary  layer  pressure  fluctuations  proceeded  on  two 
fronts:  theoretical  development  of  a  wave  guide  model  and  experimental 
verification  of  the  model  using  a  buoyant  body  facility.  The  experimental 
work  utilized  the  results  of  our  research  on  spectral  estimation.  The 
experimental  work  verified  the  wave  guide  model,  and  in  particular 
predicted  dispersion  characteristics  of  the  pressure  signal. 

The  major  purpose  of  the  experimentation  with  a  flexible  plate  in  a 
water  tunnel  was  to  Investigate  the  validity  of  the  assumptions  which 
underlie  the  weighted  Integral  approach  now  used  to  calculate  the  vibration 
response  of  water  loaded  plates  subjected  to  turbulent  boundary  layer  exci¬ 
tation.  The  most  Important  assumption  is  that  the  plate  vibration  has  so 
little  effect  on  the  turbulent  boundary  layer  that  the  power  spectral 
density  (spectrum)  of  the  pressure  field  on  the  wall  is  affected  only  to  a 
negligible  extent.  The  second  major  assumption  is  that  the  plate  response 
is  linear  so  that  the  superposition  can  be  used.  Measurement  of  both  the 
temporal  and  spatial  characteristics  of  the  pressure  field  adjacent  to  a 
rigid  wall,  as  well  as  the  vibratory  response  of  a  flexible  plate  to  the 
same  boundary  layer  excitation,  has  been  completed.  Work  is  continuing  on 
calculating  the  vibration  cross  spectra  using  dimensionless  pressure 
spectra  and  the  measured  plate  response  functions.  When  completed,  the 
results  will  be  compared  with  the  measured  vibration  cross  spectra. 

The  tip  vortex  research  consisted  of  three  parts:  a)  measurement  of 
cavitation  inception  and  fully  developed  cavitation  characteristics, 
b)  measurement  of  the  resulting  noise  intensity  spectral  characteristics, 
and  noise  source  distribution  within  the  cavitating  vortex  and  c)  LDV 
measurement  of  the  velocity  field,  construction  of  a  semi-empirical  vortex 
model  and  estimation  of  the  pressure  within  the  vortex  core. 

The  research  program  Included  the  most  exhaustive  examination  of 
cavitation  inception  in  tip  vortex  available  in  the  literature.  A  very 
subtle  interrelationship  between  viscous  effects  and  gas  content  effects 
has  been  identified  as  well  as  the  surprising  observation  that  cavitation 
inception  often  occurs  only  after  significant  levels  of  tension  have  been 
sustained.  This  infers  that  cavitation  inception  testing  may  be  facility 
dependent.  Unfortunately,  there  Is  not  enough  data  available  from 


different  facilities  to  verify  or  deny  this  possibility.  The  Important 
practical  finding  is  that  the  inception  process  at  prototype  Reynolds 
numbers  can  be  quite  different  from  the  inception  process  at  Reynolds 
numbers  typically  encountered  in  water  tunnel  testing  of  model  propellers. 

Several  different  types  of  noise  signatures  were  noted,  which  were 
related  to  different  types  of  cavitation  inception.  As  expected,  the 
levels  of  tip  vortex  cavitation  noise  are  substantially  lower  than  surface 
cavitation  when  it  occurs.  It  is  important  to  note  that  the  cavitation 
indices  for  vortex  cavitation  are  substantially  higher.  In  many  cases 
expansion  of  a  cylindrical  cavity  was  a  significant  noise  mechanism. 

A  hydrophone  array  was  utilized  to  identify  the  noise  source  location. 
In  spite  of  limitations  on  spatial  resolution  and  repeatability,  the 
results  of  this  study  clearly  indicate  that  cavitation  noise  sources  are 
concentrated  in  the  vortex  rollup  region  less  than  one  chord  length  from 
the  tip  of  the  hydrofoil. 

The  noise  intensity  associated  with  fully  developed  tip  vortex 
cavitation  is  substantially  lower  than  inception  noise,  or  surface 
cavitation  noise,  when  it  occurs.  At  high  angle  of  attack  and  under 
special  circumstances,  a  very  strong,  low  pitch  tone  occurs.  The  frequency 
of  the  tone  varied  with  cavity  core  radius  and  free  stream  velocity, 
discounting  the  thought  that  this  might  be  a  flow  Induced  vibration  of  the 
hydrofoil  itself.  Although  no  conclusive  statement  could  be  made  from  the 
present  set  of  data,  it  was  speculated  that  the  oscillation  of  the  twisted 
ribbon  type  cavity  is  responsible  for  this  tone. 

Although  these  tests  were  carried  out  in  a  highly  reverberant  water 
tunnel,  experiments  and  computer  simulations  demonstrated  that  the  shape 
of  the  measured  spectra  was  unaffected  by  reverberation. 

This  project  was  a  joint  effort  of  five  principal  Investigators  from 
three  different  departments.  All  the  experimental  work  was  carried  out  at 
the  St.  Anthony  Falls  Hydraulic  Laboratory.  In  addition  to  the  five 
faculty,  two  visiting  professors,  and  several  visiting  scholars  were 
Involved.  The  project  supported  partially,  or  totally,  four  Ph.D. 
dissertations  and  four  M.S.  theses.  In  addition,  several  undergraduates 
carried  out  special  research  projects. 


II.  INTRODUCTION 


This  program  was  a  multi-disciplinary  effort  in  hydroacoustics.  New 
digital  analysis  techniques  have  been  developed  to  measure  the  low 
wave-number  spectrum  of  turbulent  boundary  layer  pressure  in  an  unbiased 
manner  using  an  acoustic  array.  Experimental  research  on  turbulent 
boundary  layer  pressure  fluctuations  was  carried  out  in  two  facilities: 
a)  a  low  noise  recirculating  water  tunnel  and  b)  a  buoyant  rising  body 
facility.  Data  from  the  water  tunnel  research  was  utilized  as  input  to  an 
experimental  program  on  fluid-structure  interaction  in  a  turbulent  flow, 
including  the  effects  of  damping.  In  addition  to  the  research  on 
turbulent  boundary  layer  pressure  fluctuations,  cavitation  in  tip  vortices 
and  the  resulting  noise,  was  studied  in  a  low  noise  recirculating  water 
tunnel.  Additional  Laser  Doppler  Veloclmetry  measurement  of  the  structure 
of  trailing  vortices  was  made  in  a  second  water  tunnel  specially  modified 
for  the  purpose.  This  work  was  complemented  by  flow  visualization  studied 
in  a  flume  and  in  a  water  tunnel.  Details  of  each  part  of  the  program  is 
documented  in  subsequent  sections  of  this  report. 
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Higuchi,  and  R.  Taghavi,  Cavitation  and  Multiphase  Flow  Forum, 
American  Society  of  Mechanical  Engineers,  June  1985. 

4.  "Reverberation  Effects  on  Cavitation  Noise,"  R.E.A.  Arndt  and  H. 
Higuchi,  Cavitation  and  Multiphase  Flow  Forum,  American  Society  of 
Mechanical  Engineers,  June  1985. 

5.  "Cavitation  Inception  in  Trailing  Vortex  from  a  Hydrofoil,"  R.E.A. 
Arndt,  H.  Higuchi,  and  S.  Jiang,  13th  IAHR  Symposium,  Section  on 
Hydraulic  Machinery,  Equipment,  and  Cavitation,  1986. 

6.  "Assessment  and  Mitigation  of  Cavitation  Pitting  in  Hydraulic 
Turbines,"  R.E.A.  Arndt,  A.  Ferreira,  P.  R.  Rodrigue,  J.  P.  Sinclair, 
and  R.  P.  Voigt,  13th  IAHR  Symposium,  Section  on  Hydraulic  Machinery, 
Equipment  and  Cavitation,  1986. 

7.  "Tip  Vortex  Cavitation,"  R.E.A.  Arndt,  H.  Higuchi,  and  T.  Ikohagi, 
International  Symposium  on  Propeller  and  Cavitation,  Wuxi,  China, 
April  1986. 

8.  "Further  Studies  of  Tip  Vortex  Cavitation,"  R.E.A.  Arndt  and  H. 
Higuchi,  1986  Cavitation  and  Multiphase  Flow  Forum,  Fluids  Engineering 
Division,  ASME,  Atlanta,  Georgia,  May  1986. 
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"Characteristics  of  Tip  Vortex  Cavitation  Noise,"  R.E.A.  Arndt,  H. 
Higuchi,  and  M.  Rogers,  ASME  International  Symposium  on  Cavitation 
Noise,  Anaheim,  California,  Dec.  1986. 


Progress  Reports 


"Research  Program  in  Hydroacoustics  -  Progress  Reports"  submitted  to 
Office  of  Naval  Research 

-  March  1  through  October  15,  1983 

-  October  15,  1985  through  January  15,  1984 

-  January  15  through  April  15,  1984 

-  April  15  through  October  15,  1984 

-  October  16,  1984  through  April  15,  1915 

Presentations 

•  Acoustical  Society  of  America,  Minneapolis,  Minnesota,  April,  1984. 

"Hydroacoustic  Research  at  the  University  of  Minnesota," 
by  R.  Arndt,  invited  lecture. 

•  Hydraulic  Engineering  Research  Institute  (ICH),  Bucharest,  Romania, 
July,  1984. 

"Cavitation  Research  at  the  St.  Anthony  Falls  Hydraulic 
Laboratory,'*,  by  R.  Arndt. 

•  5th  ASCE-EMD  Specialty  Conference,  Laramie,  Wyoming,  August  1984. 

"Cavitation  Research  and  its  Relevance  in  the  Practice  of 
Hydraulic  Engineering,"  invited  lecture  by  R.  Arndt. 

•  American  Society  of  Civil  Engineers  Hydraulics  Division  Specialty 
Conference,  Coeur  d'Alene,  Idaho,  August,  1984 

"Cavitation  in  Various  Types  of  Shear  Flow,"  by  R.  Arndt. 

•  University  of  Illinois,  Hydrosystems  Graduate  Seminar  Series, 
Urbana,  Illinois,  Oct.  1984. 

"Cavitating  Jets  and  Related  Topics,”  by  R.E.A.  Arndt, 

Invited  Speaker. 


•  Central  Water  and  Power  Research  Station,  Puna,  India,  Jan.  1985 


Series  of  lectures  on  cavitation,  by  R.  Arndt. 


•  Acoustical  Society  of  American,  Minneapolis,  Minnesota,  Oct.  1984 


"Measured  Wall  Pressure  Cross  Spectra  Due  to  a  Turbulent 
Boundary  Layer,"  by  N.  Narayan  and  R.  Plunkett. 

"Measurements  of  Low-Wavenumber  Spectra  of  Wall  Turbulence  in 
a  Buoyant  Body  Facility,"  by  D.  Long,  M.  Wurtz,  and  R. 
Lambert, 

"Spectral  Estimation  of  the  Low-Wavenumber  Region  of  a 
Turbulent  Boundary  Layer:  Maximum  Likelihood  Method,"  by 
G.  Wakefield  and  M.  Kaveh. 

"Comparative  Examination  of  Tip  Vortex  and  Surface 
Cavitation  Noise,"  by  H.  Higuchl,  S.  Forbes,  and  R.E.A. 
Arndt. 


•  37th  Meeting  of  Division  of  Fluid  Dynamics,  American  Physical  Society, 
Providence,  Rhode  Island,  November  1984. 

"Estimation  of  the  Low-Wavenumber  Region  of  a  Turbulent 
Boundary  Layer  Using  the  Maximum  Likelihood  Method,"  by 
G.  Wakefield  and  M.  Kaveh. 

"Comparative  Examination  of  Tip  Vortex  and  Surface 
Cavitation  Noise,"  by  H.  Higuchl  and  R.E.A.  Arndt. 

"Measured  Wall  Pressure  Cross  Spectra  Due  to  a  Turbulent 
Boundary  Layer,"  by  N.  Narayan  and  R.  Plunkett. 

"Low  Wave  Number  Spectra  of  Wall  Pressure  on  a  Buoyant 
Underwater  Body,"  by  R.  Lambert,  D.  Long,  and  M.  Wurtz. 


•  IEEE  International  Conference  on  Acoustics,  Speech,  and  Signal 
Processing,  ICASSP  85. 

"Frequency-Wavenumber  Spectral  Estimation  of  Non-Planar 
Random  Fields,"  by  G.  Wakefield  and  M.  Kaveh. 


•  St.  Anthony  Falls  Hydraulic  Laboratory  Seminar  Series. 

"Cavitation  in  Vortical  Motion,"  by  R.  Arndt,  Nov.  1984. 

"The  Hydroacoustic  Program  at  St.  Anthony  Falls  Hydraulic 
Laboratory,  by  R.  Arndt,  March  1984. 


•  24th  AIAA  Aerospace  Sciences  Meeting,  Reno,  Nevada,  Jan.  1985. 

"An  Experimental  Investigation  of  Vortex  Roll-Up  for  an 
Elllptically-Loaded  Wing,"  by  H.  Higuchl,  J.  C.  Quadrelli, 
and  C.  Farell,  AIAA  Paper  86-0562. 


•  Henkel  Corporation  Seminar,  Minneapolis,  Minnesota,  May,  1985. 

"Cavitation  Research  at  the  St.  Anthony  Falls  Hydraulic 
Laboratory,  by  Roger  E.A.  Arndt. 


•  Joint  ASCE/ASME  Mechanics  Conference,  Albuquerque,  New  Mexico,  June 
1985. 

"Reverberation  Effects  on  Cavitation  Noise,"  by  H. 

Higuchi . 

"Sound  Source  Location  in  Cavitating  Tip  Vortices,"  by  H. 
Higuchi. 

"Tip  Vortex  Cavitation,"  by  R.E.A.  Arndt. 

"Cavitation  in  Various  Types  of  Shear  Flow,"  by  R.E.A. 

Arndt. 

"Utility  Survey  of  Cavitation  Erosion  in  Hydraulic 
Turbines,"  by  Roger  E.A.  Arndt. 


•  Megatrends  in  Hydraulic  Research  for  the  Future,  Fort  Collins, 
Colorado,  July  1985. 

"Cavitation  Research:  Synopsis  and  Perspective,"  invited 
lecture  by  R.  Arndt. 


•  21st  1AHR  Congress,  Hydraulic  Structures  Seminar,  Melbourne, 
Australia,  Aug.  1985. 

"Cavitation  in  Hydraulic  Structures,"  by  R.  Arndt. 


•  International  Symposium  on  Propeller  and  Cavitation,  Wuxi,  China, 
April  1986. 

"Tip  Vortex  Cavitation,"  by  R.  Arndt. 


•  Institute  of  High  Speed  Mechanics,  Sendai,  Japan,  April,  1986. 

"Cavitation  Research  at  St.  Anthony  Falls  Hydraulic 
Laboratory,"  invited  lecture  to  facalty  and  students,  by 
R.  Arndt. 


•  American  Society  of  Mechanical  Engineers  Cavitation  and  Multiphase 
Flows  in  Fluid  Machinery,  Atlanta,  Georgia,  May  1986. 

"Special  Comments  on  a  Recent  Symposium  on  Cavitation  held 
in  China,"  by  R.  Arndt. 


•  NSF  Romanian-U.S.  Workshop  in  Hydraulic  Engineering  Research, 
Bucharest,  Romania,  July  1986. 

"Cavitation  Research  in  Hydraulic  Engineering,"  by  R. 

Arndt. 


12th  American  Hydraulics  Congress  of  the  IAHR  Short  Course,  Sao  Paulo, 
Brazil,  Aug.  1986. 

"Cavitation,"  by  R.  Arndt 


American  Physical  Society,  Columbus,  Ohio,  Nov.  1986. 


"Viscous  Effects  on  Tip  Vortex  Roll-Up  and  Cavitation,"  by 
R.  Arndt. 


List  of  Honors,  Awards 

Robert  Lambert  received  the  first  John  C.  Johnson  Memorial  Education 
Award  from  the  Institute  of  Noise  Control  Engineering  (INCE)  "for 
excellence  in  teaching  of  acoustics  and  noise  control  engineering," 
Dec.  1984. 
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1.  Theoretical  end  Algorithmic  Aspects  of  Frequency-Wavenumber 
Spectral  Estimation  of  Turbulent  Boundary  Layers 

M.  Kaveh,  Principal  Investigator 


The  theoretical  and  algorithmic  aspects  of  our  research  on  the  estima¬ 
tion  of  frequency-wavenumber  spectra  of  turbulent  boundary  layers  (TBL) 
proceeded  as  follows.  First  a  statistical  model  of  the  cross-spectral 
density  matrices  (CSDM)  of  TBL  random  fields,  as  proposed  by  Chase  and 
Corcos  were  developed.  These  models  explicitly  included  dependence  on  the 
array  geometry,  measurement  time  and  signal-to-noise  ratios  at  each  sensor 
element  and  formed  the  backbone  of  the  simulations  that  were  carried  out  to 
assess  the  performances  of  various  estimators.  Next,  the  most  promising 
"modern"  spectral  estimation  techniques  were  evaluated  for  suitability  of 
representing  the  TBL  spectral  models.  Finally,  new  spectral  estimation 
techniques  were  developed  that  were  aimed  at  high-fidelity,  rather  than 
high-resolution  spectral  estimation  of  the  type  of  spectra  expected  from 
the  TBL  random  fields.  Where  possible,  theoretical  analyses  of  the  perfor¬ 
mance  of  the  estimators  were  made.  In  all  cases,  extensive  simulations 
were  carried  out  to  evaluate  both  the  spectral  fidelity  and  variability  of 
the  estimators. 

The  results  are  summarized  below.  Details  of  the  methods  and  their 
performance  can  be  found  in  the  attached  published  papers  and  the  Ph.D. 
dissertation  by  G.  Wakefield. 

1.  Traditional  approaches  for  the  estimation  of  the  TBL  wavenumber 
spectra  are  based  on  the  principle  of  beamforming  by  an  array  of  sensors. 
Unfortunately,  the  relatively  small  array  aperture  of  an  experimental  array 
does  not  allow  for  sufficient  selectivity  and  resolution  to  eliminate 
leakage  of  power  from  the  convective  and  sonic  regions  of  spectra  into  the 
low-wavenumber  region,  whose  characteristics  are  of  primary  interest.  Our 
estimation  results  based  on  the  beamformer  were  used  as  the  baseline  for 
evaluating  the  performance  of  other  methods. 

2.  The  maximum  Entropy  Method  is  perhaps  the  most  widely  used  tech¬ 
nique  for  high  resolution  spectral  estimation.  The  performance  of  this 
method  has  usually  been  evaluated  in  the  peak  regions  of  peaky  spectra. 
MEM  was  applied  to  the  TBL  models  with  mixed  results.  For  sufficiently 
high  signal-to-noise  ratios,  the  fidelity  of  the  estimates  were  found  to  be 
quite  good.  Statistical  performance  analyses,  however,  showed  the  MEM 
spectral  variance  in  the  low-wavenumber  regions  to  be  too  large,  rendering 
this  method  unreliable  for  this  application. 

3.  Capon's  Maximum  Likelihood  (ML)  proved  to  be  the  most  stable  known 
modern  spectral  estimation  technique  for  the  estimation  of  the  TBL  wave- 
number  spectra.  This  technique  produced  spectral  estimates  that  were  more 


14 


biased,  In  Che  low  wavenumber  regions,  Chan  esClmaCes  generated  by  Che  MEM. 
The  MLM  estimator  variance,  however,  was  substantially  lower  than  that  of 
the  MEM.  Furthermore  the  MLM  possesses  a  constrained  spatial  filtering 
Interpretation  that  Is  Intuitively  pleasing  and  devoid  of  the  Imposition  of 
a  parametric  model  on  the  spectra. 

4.  As  a  means  of  Increasing  the  effective  aperture  of  an  array,  with 
fixed  number  of  sensors,  nonuniform  array  alternatives  were  Investigated  In 
conjunction  with  the  MLM  estimator.  It  was  found  that,  whereas  array 
geometries  such  as  Minimum-Redundancy  Increase  resolution  In  the  case  of 
multiple  plane  waves,  they  do  not  offer  higher  fidelity  In  estimating  smooth 
spectra  such  as  those  of  the  TBL's.  This  was  attributed  to  Che  Increase  In 
the  average  level  of  the  array  sidelobes  when  a  narrower  mainlobe  Is 
generated  by  the  extended  array  aperture. 

5.  The  MLM  was  extended  to  operate  on  a  subspace  of  the  elgenstruc- 
ture  of  the  CSDM.  The  subspace  was  chosen  to  include  the  CSDM  eigenvectors 
whose  Fourier  transforms  primarily  contribute  to  the  low  wavenumber  region. 
This  Implied  the  elimination  of  the  eigenvectors  corresponding  to  several 
of  the  largest  eigenvalues  that  mainly  form  the  convective  specral  peak. 
This  technique  Improved  the  fidelity  of  the  MLM  In  the  low  wavenumber 
regions.  Determination  of  the  significant  eigenvalues  for  an  unknown 
spectral  model,  however,  proved  to  be  a  definite  shortcoming  of  this  tech¬ 
nique. 

6.  One  feature  of  the  MLM  is  that,  unlike  the  MEM  and  the  beamformer, 
the  inverse  Fourier  transform  of  the  estimated  spectrum  does  not  match  the 
CSDM  from  which  the  MLM  spectrum  was  computed.  A  new  method,  denoting  the 
Covariance  Constrained  MLM  (CCMLM),  was  developed  to  alleviate  this  short¬ 
coming  and  in  the  process  to  Improve  the  fidelity  of  the  MLM.  The  CCMLM 
uses  the  MLM  in  an  Iterative  fashion  to  generate  an  estimate  of  the 
spectrum  that  Is  asymptotically  (large  number  of  Iterations)  matched  to  the 
measured  CSDM.  A  practical,  stable  iterative  algorithm  was  designed  and 
the  performance  of  the  method  was  evaluated  both  theoretically  and  through 
extensive  numerical  simulations.  The  CCMLM,  which  includes  the  MLM  as  a 
special  case,  proved  to  be  the  best  technique  that  we  evaluated  or 
developed  for  the  estimation  of  spectra  such  as  those  suggested  by  Corcos 
and  Chase. 


2.  Wa 


ter/rrtfWB cy  Spectral  Estimation 


of  Vail  Pressure  Fluctuations 


R.  Lambert,  Principal  Investigator 


Measurements 


The  wall  pressure  fluctuations  due  to  turbulence  have  been  measured 
over  the  useful  range  of  frequencies  from  approximately  45  Hz  to  400  Hz  and 
wavenumbers  ranging  from  0  to  400  m“l.  The  measurements  were  obtained 
using  a  10  sensor  uniform  spaced  linear  array  with  a  sensor  spacing 
D  “  0.0158  M  center  to  center  and  aperture  of  0.14  meters.  The  data 
obtained  has  been  normalized  to  dimensionless  quantities  using  the 
following  flow  profile  parameters:  V  -  8.55  M/s;  6  ■  0.025  M;  6*  - 
0.125  «  -  0.00313  M.  " 

The  data  In  Figure  V2.1  through  V2.4  are  wavenumber  estimates  at 
various  selected  frequencies.  Figure  V2.1  shows  the  relative  performance 
of  the  three  spectral  estimation  techniques.  The  traditional  Blackman- 
Tukey  method  (BT);  the  Maximum  Likelihood  Method  (MLM);  and  the  Covariance 
Constrained  MLM  (CCMLM) .  The  BT  estimate  Is  dominated  by  a  large  main 
lobe  and  side  lobe  leakage  from  the  acoustic  peak  located  at  k  ■  t  w/c 
where  c  is  the  speed  of  sound  In  water.  In  Figure  V2.1  the  BT  estimate 
cannot  resolve  the  convective  peak  because  of  side  lobe  leakage.  (The  MLM 
estimate  on  the  other  hand  shows  a  significant  reduction  In  leakage  from 
the  acoustic  peak  clearly  resolves  the  convective  peak.)  The  CCMLM  esti¬ 
mate  further  Improves  the  fidelity  In  the  low  wavenumber  regions,  by 
approximately  12dB  at  this  frequency.  At  higher  frequencies,  however,  the 
CCMLM  Improvements  over  the  MLM  estimates  become  less  pronounced.  This 
results  because  the  MLM  Is  "doing  a  good  Job"  with  regards  to  matching  the 
Inverse  Fourier  transform  of  the  estimate  to  the  original  Cross  Spectral 
Density  Matrix  (CSDM). 

A  caution  In  using  the  CCMLM  is  that  as  the  estimate  approaches  Its 
error  reduction  limit,  anomalous  nulls  (or  very  small  spectral  levels) 
start  to  show  up  in  the  measured  spectrum  (see  Figure  V2.3).  It  is  also 
worth  mentioning  that  although  the  MLM  estimator  is  an  upper  bound  of  the 
true  spectrum  by  design,  the  CCMLM  Is  not,  and  It  Is  possible  for  the  esti¬ 
mate  to  be  significantly  below  the  true  spectrum  In  a  given  region  when 
nulling  occurs. 

The  wavenumber  spectra  have  been  assumed  to  be  predominantly  one  sided 
and  the  estimate  is  periodic  about  k  ■  tr/D  ■  0.0158  m-*.  This  assumption 
gives  approximately  a  doubling  In  thl  range  of  wavenumber  that  can  be  uti¬ 
lized  for  a  given  sensor  spacing.  Figure  V2.4  shows  the  CCMLM  estimate 
(which  In  this  case  is  the  MLM  estimate)  Just  prior  to  the  convective  peak 
aliasing  into  the  acoustic  peak. 
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Figure  V2.I.  Comparison  of  the  wavenumber-frequency  spectra  from  the  BT 
and  CCMLM  estimators  at  w*  =  0.19,  f  =  84.5  hz. 
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Figure 


Figure  V2.3.  MLM  and  a  nulling  CCMLM  estimate  at  «  =  0.58,  f  =  253  hz 
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The  buoyant  body  facility,  although  lacking  pump  or  turbine  noise,  was 
plagued  with  acoustic  noise  probably  generated  from  the  wake  and  amplified 
by  the  highly  reverberant  cylindrical  structure  of  the  buoyant  body  faci¬ 
lity.  Figure  V2.5  shows  the  normalized  auto  power  spectrum  averaged  over 
all  sensors.  The  spectrum  is  dominated  by  acoustic  energy  below  approxi¬ 
mately  a)*  *  0.3.  This  can  be  verified  by  the  phase  spectrum  of  the  first 
lag  shown  in  Figure  V2.6.  Note  that  the  highly  spatially  correlated 
acoustic  energy  causes  the  phase  to  be  nearly  equal  to  zero  below  u*  *  0.3. 
The  averaging  effect  of  finite  transducer  size  (d*  ■  3.0)  start  coming  into 
effect  at  frequencies  greater  than  approximately  w*  =  1.0,  see  Fig.  V2.5. 
Although  the  data  in  V2.5  and  V2.6  might  be  disappointing,  there  Is  no 
doubt  that  acoustic  noise  dominated  the  experimental  results.  However,  it 
should  be  noted  that  the  wavenumber  estimates  employing  CCMLM  distinguishes 
acoustic  energy  and  turbulent  energy  to  a  certain  extent.  The  limitation 
appears  to  be  side  lobe  leakage.  This  is  mentioned  so  that  the  data  on  the 
location  of  the  convective  peaks  employed  in  Part  B  should  not  be 
discounted  because  of  acoustic  contamination,  although  the  absolute  levels 
of  the  peaks  may  have  been  affected  somewhat  because  of  the  normalization 
and  calibration  techniques  used. 


Convective  Waveguide  Model-Wavenumber  and  Velocity  Predictions 

An  effort  was  made  early  on  in  this  research  program  to  study  and  be 
cognizant  of  any  possible  wavenumber /frequency  dispersion  effects  that 
could  possibly  be  observed  from  spectral  experiments  in  the  buoyant  body 
facility.  It  had  been  observed  earlier  that  phase  convective  velocity 
measurements  in  wind  tunnels  did  follow  a  characteristic  in  frequency  that 
was  dispersive.  In  fact,  in  both  wind  and  water  tunnels  the  phase  velocity 
appeared  to  vanish  at  zero  frequency,  peak  at  some  low  frequency,  and  then 
level  off  at  high  frequencies  at  some  fixed  percentage  of  the  mid-stream 
fluid  velocity,  usually  between  0.6  and  0.8.  If  the  peaks  in  the 
wavenumber/frequency  spectrum,  which  lie  on  the  convective  ridge,  could  be 
observed  with  some  fidelity,  then  any  dispersive  effects  in  the  buoyant 
body  facility  also  could  be  validated. 

Narayan  and  Plunkett  observed  from  phase  measurements  in  a  water  tun- 
jj  nel  that  the  convective  phase  velocity  did  appear  to  peak  at  low  frequen¬ 

cies  and  then  level  off  at  higher  frequencies  at  about  0.6  of  the  free 
stream  velocity  (see  Fig.  21  in  TR-AEM-85-1-RP  dated  14  February  1985).  A 
similar  observation  in  a  wind  tunnel  had  been  made  by  Blake,  and  his  data 
taken  at  values  of  ffi  ■  ox5*/V  <  0.50  exhibited  a  peak  and  a  definite 
1  trend  towards  zero  (see  Journal  of  Fluid  Mechanics,  44,  1970,  pp.  637-660). 

i  Blake' 8  data  are  reproduced  in  Fig.  V2.7,  attached. 

In  order  to  help  understand  these  and  related  phenomena,  a  set  of  nine 
Internal  memoranda  were  prepared,  of  which  the  last  five  (V  thru  IX)  are 
appended  to  this  report  (Appendix  B) .  The  first  four  tracts  were  concerned 
with  measurement  techniques  and  data  reduction  methods,  including  such 
topics  as  frequency  corrections  for  finite  transducer  size. 

Tract  V  titled  "Estimate  of  the  Low  Frequency  Power  Spectrum  of  Wall 
Pressure  Fluctuations"  details  some  of  the  basic  ideas  in  what  is  herein 
referred  to  as  a  convective  waveguide  model.  It  starts  with  a  review  of 
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2.5.  Normalized  power  spectrum  showing  acoustic  contamination 
below  ui*  =  0.3  and  finite  transducer  averaging  effects  a 
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Experimental  plot  of  the  location  of  the  convective  peaks  versus  the  normalized  frequency 
£}  for  the  buoyant  body  facility.  The  solid  line  is  a  fit  of  Eqs.  5  and  7  for  the  para¬ 
meters  given  in  the  legend. 
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some  Ideas  embodied  in  Che  Kolmogoroff  Theory  of  Turbulence  as  interpreted 
and  explained  by  G.  K.  Batchelor  in  his  monogram.  The  waveguide  model  was 
developed  on  the  hypothesis  that  in  a:,  experiment  of  finite  dimensions 
(such  as  in  a  water  tunnel  or  buoyant  body  facility)  the  wavenumber 
spectrum  would  actually  cut-off  below  some  finite  value  of  wavenumber,  k  , 
greater  than  zero.  This  was  based  on  the  intuitive  idea  that  eddies  above 
a  certain  size  (or  below  a  certain  wavenumber)  were  too  large  to  be  con- 
vected  by  an  incompressible  fluid.  Since  water  has  some  compressibility, 
however  small,  there  probably  would  not  be  a  sharp  cut-off  or  cut-on  as  the 
case  may  be,  but  rather  the  phenomena  would  occur  gradually  as  the  observed 
frequency  is  reduced. 

As  outlined  in  tract  V,this  development  led  to  a  rather  simple  rela¬ 
tionship  between  frequency  ta  and  wavenumber  k  of  the  form 

a2  -  V  2 (k2 -k  2)  k  >  k  (1) 

coo 

where  V is  called  an  •  convection  velocity  and  the  cut-on  wavenumber 


c 


where  2nA  is  the  perimeter  of  the  largest  eddy  in  the  guide.  In  the 
absence  of  better  information  it  was  decided  to  set  A  *  d  where  6  is  a 
boundary  layer  thickness  defined  in  the  jsual  way.  The  parameter  V 
8 till  is  interpreted  as  an  eddy  convectior  velocity  but  not  a  phase  velo¬ 
city  and  hence  is  unobservable  directly  from  phase  measurements  (except 
perhaps  at  very  high  frequencies) . 

However,  the  observed  phase  velocity 


1 s  well-defined. 


In  terms  of  normalized  variables,  Eq.  3  can  be  rewritten  as 
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a  for®  that  requires  knowledge  of  how  V  /V*  varies  with  frequency  and  a 
numerical  value  for  6*/S  where  5*  isc  the  displacement  boundary  layer 
thickness  and  where  Q  ■  a6*/Vm  .  If  Vg/V*  ■  constant,  then  Eq.  4  predicts 
that  VD/V»  will  increase  from  zero  and  gradually  approach  the  value 


Vc  as 


is  increased. 


and  Plunkett  reveal  a  peak  in  the  Vp/V, 
Si  which  requires  a 
that  a 


'C/V< 


choice  for  Vc/Va 


However,  Blake's  data  and  the  data  of  Narayan 
characteristic  at  some  low  value  of 
also  to  vary  with  Q  .  By  trial  it  was  determined 
of  the  form 
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(5) 


Si  *  0 


would  provide  a  satisfactory  fit  of  the  Vp/Va.  characteristic,  Eq.  4,  to 
observed  phase  velocity  data.  Here  (Vp/V*,^  ■  constant  is  the  asymptotic 
value  of  Vp/Vo,  observed  at  high  frequencies*.  Figure  V2.7  shows  a  fit  of 
Blake's  data  for  the  choice  (Vp/Va,)^  ■  0.60.  Note  Eqs.  4  and  5  also  pre¬ 
dict  a  peak  in  Vp/V„  of  about  0.83  at  Hi  =  0.30  is  almost  in  exact  agreement 
with  Blake' 8  measurements.  This  was  the  first  indication  that  the  convec¬ 
tive  waveguide  model  could  be  validated. 

The  data  of  Narayan  and  Plunkett  also  reveal  a  similar  peak  of  about 
the  same  value  at  certain  measurement  distances.  However,  they  did  not 
normalize  their  frequency  scale  so  a  direct  comparison  was  not  made. 

Further  properties  of  the  convective  waveguide  model  were  then  deve¬ 
loped  and  of  special  note  is  a  suggested  experiment  in  Appendix  B  of  tract 
VIII.  As  noted  there,  the  convective  waveguide  model  predicted  that  the 
longitudinal  wavenumber/frequency  power  spectrum  P(k^ ,  0,  w)  would  exhibit 
a  peak  at  the  wavenumber  value  kp  given  by 


or  in  terms  of  normalized  variables 


Note  Eq.  7  predicts  that  kp  ♦  1/6  as  0  0  and  kp  >  u>/Vc  as  Q  ♦  •  by 
employing  Eq.  5.  Thus,  the  convective  waveguide  model  exhibits  a  form  of 
frequency  dispersion  not  found  in  other  theoretical  formulations. 

The  suggested  experiment  was  to  isolate  the  peaks  in  the  P(ki,  0,  w) 
spectrum  characteristic  as  accurately  as  possible  and  then  plot  kp  versus 
(t)  in  order  to  see  if  the  resulting  dispersion  characteristic  did  indeed 


follow  Eq.  7  and  cutoff  as  predicted .  It  was  felt  that  information  about 
the  location  of  peaks  In  the  convective  ridge  of  the  longitudinal 
wavenumber/frequency  spectrum  would  be  the  best  information  about  the  model 
that  could  be  obtained  even  before  the  experiments  were  conducted.  It  was 
also  felt  that  other  critical  information  about  the  P(ki,  0,  m)  charac¬ 
teristic  would  be  masked  by  acoustic  background  noise. 

The  suggested  experiment  was  performed  in  the  buoyant  body  facility 
and  the  peaks  plotted  as  a  function  of  frequency,  see  Fig.  V2.8.  It  was 
decided  to  use  a  nominal  value  of  S*/6  *  0.125  in  the  absence  of  experimen¬ 
ts  data  on  the  actual  boundary  layer  profile.  The  data  of  Fig.  V2.8 
approach  a  cutoff  wavenumber  Vcq  «  40  m-*  as  ai  tends  to  zero  via  empirical 
fits  of  either  Eq.  6  or  Eq.  7.  This  result  then  yields  a  value  of  5  « 
0.025  cm  for  the  boundary  layer  flow  over  the  buoyant  body  in  the  region 
where  the  spectral  measurements  were  made.  This  value  is  about  2.5  times 
the  thickness  expected  on  the  basis  of  discussions  heretofore. 

The  high  frequency  data  in  Fig.  V2.8  also  reveals  a  value  of  Vc/V„  - 
0.80  using  a  measured  value  of  Vw  =*  8.55  m/s  for  the  speed  of  the  buoyant 
body. 


A  predicted  dispersion  diagram  plotting  o>  versus  x  is  shown  1’ 
Fig.  V2.9.  Two  curves  are  plotted  for  comparison  purposes,  one  showing  > 
plot  of  u  ■  Vck  based  on  Eq.  5  (dotted  line)  and  the  other  showing  u  ver¬ 
sus  k  based  on  Eqs.  1  and  7  revealing  a  sharp  cut-on  at  k  ■  kQ .  The 
characteristics  are  based  on  empirical  data  for  the  buoyant  body. 

Finally,  a  predicted  Vp/V*,  characteristic  based  on  this  empirical 
information  for  the  buoyant  body  is  plotted  in  Fig.  V2.10.  There  exists  a 
predicted  peak  of  about  0.9  that  occurs  at  S  ■  0.30  vhich  is  in  fair 
agreement  with  the  work  of  others.  Note,  however,  the  predicted  high  fre¬ 
quency  asymptote  (Vp/V,,,,)^  *  0.80  which  is  somewhat  higher  than  usually 
observed  in  wind  tunnels  and  in  the  water  tunnel  data  of  Narayan  and 
Plunkett. 

Various  other  theoretical  predictions  based  on  a  convective  waveguide 
model  were  made  in  tracts  V-IX.  Of  note  are  predictions  of  power  spectrum 
in  tract  V  and  an  estimate  of  the  low  wavenumber/frequency  spectrum  in 
tract  VI.  All  of  the  numerical  work  in  tracts  V-IX  was  based  on  empirical 
information  obtained  from  wind  tunnel  data.  The  various  characteristics 
would  have  to  be  recomputed  in  order  to  make  comparisons  with  either  water 
tunnel  or  buoyant  body  spectrum  measurements  especially  in  view  of  the 
thicker  boundary  layer.  Values  of  6  and  6*  influence  most  of  the 
numerical  predictions. 
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Predicted  wave  dispersion  characteristics  shown  by  solid  line  based  on  Eq.  1 
for  the  buoyant  body  and  the  equation  w  =»  V  k  based  on  Eq.  5  as  shown  by  the 
dotted  line.  Note  the  dotted  line  is  not  a  straight  line. 


3*  Structural  Response  to  Pressure  Fluctuation 


R.  Plunkett,  Principal  Investigator 


Introduction 


The  major  purpose  of  this  project  was  to  investigate  the  validity  of 
the  assumptions  which  underlie  the  weighted  integral  approach  which  is  now 
used  to  calculate  the  vibration  response  of  water  loaded  plates  subjected 
to  turbulent  boundary  layer  excitation.  The  most  important  assumption  is 
that  the  plate  vibration  has  so  little  effect  on  the  turbulent  boundary 
layer  that  the  power  spectral  density  (spectrum)  of  the  pressure  field  on 
the  wall  is  affected  only  to  a  negligible  extent.  The  second  major  assump¬ 
tion  is  that  the  plate  response  is  linear  so  that  a  superposition  can  be 
used. 


These  two  assumptions  are  sufficient  to  establish  the  integral 
equation  for  the  cross  spectrum  of  the  vibration  response  of  the  plate 
(Lin,  1967). 

**.yy(*^ ****  ^  ^ K  A  ^  *P  2*® 2 ^ (P i >9  2 j)d^Cp 2^  (^) 

where  is  the  cross  spectrum  of  the  normal  pressure 

exerted  by  the  boundary  layer  on  the  non¬ 
vibrating  plate. 

H(x,p;u>)  is  the  response  of  the  plate  of  x  due  to  a 
unit  sinusoidal  force  at  p  . 

and  Syy(x^  ,*2  ;w)  *s  t*ie  cross  spectrum  of  the  response. 

Tack  and  Lambert  (1962)  have  demonstrated  that  this  approach  accu¬ 
rately  predicts  the  vibration  level  for  a  narrow  beam  in  air  exposed  to  a 
turbulent  boundary  layer  in  a  wind  tunnel.  Most  of  the  papers  in  the  open 
literature  are  concerned  with  ways  of  approximating  H  and  r  so  as  to  make 
the  Integration  relatively  simple  (Powell,  1958;  Chandirlmani,  1977).  Our 
primary  objective  was  to  compare  the  measured  response  spectrum  of  a  flat 
plate  with  that  predicted  from  the  measured  cross  spectra  of  the  pressure 
field  and  the  measured  response  functions  of  the  plate.  In  this  way,  for 
at  least  one  case,  we  will  be  able  to  determine  what  the  controlling 
variables  are  and  what  accuracy  is  necessary  in  the  two  cross  functions  to 
obtain  acceptable  accuracy  for  the  vibration  level. 

Work  on  this  project  has  been  reported  in  previous  progress  reports 
(Arndt,  1985).  Up  to  November  1,  1985,  the  pressure  cross  spectrum  was 


measured  and  preliminary  results  were  obtained  for  the  vibration 
autospectrum.  Since  then,  vibration  cross  spectra  have  been  measured  for 
the  vibrating  plate,  the  pressure  cross  spectrum  has  been  measured  upstream 
and  downstream  of  the  vibrating  plate,  and  the  boundary  layer  velocity  pro¬ 
file  has  been  measured.  All  of  this  work  is  now  being  written  for  a  thesis 
in  partial  fulfillment  of  a  Ph.D.  for  N.  Narayan.  The  work  remaining  to  be 
done  is  the  actual  calculation  of  the  predicted  PSD  for  the  vibration,  its 
comparison  with  the  measured  values,  and  the  preparation  of  the  thesis  in 
its  final  form.  When  that  is  finished,  the  major  results  will  be  embodied 
in  papers  to  be  submitted  to  serial  journals  in  the  open  literature. 


All  of  the  measurements  were  made  in  the  7.5  inch  square  low  noise 
water  tunnel  of  the  St.  Anthony  Falls  Hydraulic  Laboratory.  Pressure 
spectrum  measurements  were  first  made  on  the  surface  of  the  thick 
plexiglass  celling  of  the  tunnel.  The  measured  vibration  levels  in  this 
wall  were  so  low  that  the  voltage  induced  in  the  pressure  transducers  by 
the  vibration  contributed  less  than  1%  (-  40  dB)  to  the  measured  pressure 
levels.  The  pressure  transducers  were  held  in  rotatable  holders  so  that 
the  cross  spectra  could  be  measured  for  a  grid  of  pairs  of  points  covering 
180  mm  axial  separation  and  50  mm  transverse.  The  results  were  reported  in 
an  ASHE  symposium  volume  (Narayan  and  Plunkett,  1985). 


The  roof  of  the  tunnel  was  then  modified  to  take  a  stainless  steel 
plate  200  mm  long  by  40  mm  wide  and  1.2  mm  thick.  The  plate  was  backed  by 
a  chamber  filled  with  air  and  pressurized  to  the  static  pressure  of  the 
tunnel.  Five  accelerometers  were  mounted  on  the  back  of  the  plate  to 
measure  the  vibration  cross  spectra  on  the  centerline  for  axial  separa¬ 
tions.  Pressure  transducers  were  mounted  just  upstream  and  just  downstream 
of  the  vibrating  plate  to  determine  if  the  plate  vibration  affected  the 
pressure  cross  spectra.  Measurements  were  made  at  flow  velocities  of  7,  9, 
12,  and  15  m/sec. 


The  dynamic  response  functions  from  50  to  2000  Hz  were  measured  for 
the  plate  mounted  in  the  ceiling  separator  in  air,  in  a  static  water  tank, 
and  in  the  water  tunnel  by  exciting  each  of  the  five  accelerometer  mounting 
points  sinusoidally  through  a  mechanical  vibration  impedance  head.  The 
mode  shapes  and  natural  frequencies  were  deduced  from  the  response  func¬ 
tions.  The  first  five  mode  shapes  had  no  longitudinal  nodes  and  from  zero 
to  four  transverse  nodes.  The  water  loaded  frequencies  ranged  from  about 
460  Hz  to  1500  Hz  which  agreed  reasonably  well  with  the  predicted  values 
when  the  edge  fixities  were  suitably  adjusted.  The  entrained  water  mass 
loading  had  about  one  third  the  effect  of  the  accelerometer  mass  so  that 
there  was  not  much  difference  between  the  air  and  water  loaded  frequencies. 
The  amplification  factor  at  resonance  (a  measure  of  damping)  was  about  10 
to  20  for  the  air  loaded  system  and  about  3  to  5  for  the  installation  in 
the  water  tunnel. 


Results 


i 
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The  velocity  profile  in  the  boundary  layer  was  measured  with  a  dynamic 
pressure  probe  which  had  a  rectangular  opening  0.2  mm  high  by  4  mm  thick 
and  a  ground  down  lip  thickness  of  0.02  mm.  The  probe  was  traversed  from 
just  touching  the  wall  to  39  mm  away.  A  plot  of  the  measured  points  at  a 
free  stream  velocity  of  3.8  m/sec  is  shown  in  Figure  V3.1  along  with  the 
predicted  profile  based  on  the  l/7th  power  law. 

There  was  negligible  difference  (usually  less  than  1  dB)  for  the  auto 
and  cross  spectra  measured  just  upstream  and  just  downstream  of  the 

vibrating  plate.  The  plate  vibration  level  was  about  0.01  of  the  boundary 

layer  thickness. 

The  general  characteristics  of  the  plate  vibration  acceleration 
response  autospectrum  (Fig.  V3.2)  is  a  rather  smooth  curve  falling  about 
20  dB  from  100  Hz  to  2500  Hz  with  two  peaks  at  750  Hz  (3rd  mode)  and  1500 

Hz  (5th  mode).  The  peaks  rise  about  10  dB  above  the  general  trend. 

Subject  to  verification  in  detail,  the  general  response  seems  to  be  caused 
by  the  reverberant  acoustic  field  and  the  two  peaks  are  caused  by  the  con- 
vected  turbulent  field. 

The  transverse  cross  correlation  pressure  field  among  the  three 
pressure  pickups  in  each  of  three  positions  (two  upstream  and  one 
downstream)  showed  a  much  greater  correlation  length  with  the  fixed  pickups 
than  had  been  measured  before  with  the  transducers  in  the  rotatable 
holders.  With  no  axial  separation  the  correlation  length  for  y*  *  a>y/Uc 
was  about  0.6  for  the  measurements  made  with  the  rotatable  holders.  It  is 
about  2  for  the  fixed  transducers.  These  new  results  make  it  possible  to 
fit  the  dimensionless  cross  spectrum  with  a  function  of  the  form: 

T (x,y ,  ;w)  -  r1(0,0;u)r2(x*,0)r3(0,y*) 

jl2  jlO 

where  r2  ■  exp(-ux  )  and  -  exp(-8y  ) 
o  -  0.0070;  B  -  0.183; 
and  x*  ■  o)x/Uc,  y*  *  wy/Uc. 


Corcos  (1964)  had  suggested  a  similar  expression  with  |x*|  and  |y* |  instead 
of  the  quadratic  terms.  The  difference  is  most  pronounced  for  small  values 
of  x*  and  y*  ;  a  careful  examination  of  our  short  distance  data  does  not 
show  the  linear  falloff  that  would  be  necessary  for  the  absolute  value  for¬ 
mulation.  His  formulation  has  the  added  disadvantage  of  being  difficult  to 
deal  with  analytically.  While  it  would  be  possible  to  improve  the  fit 
slightly  by  having  some  interaction  between  the  axial  (x*)  and  transverse 
(y*)  dimensionless  coordinates,  the  results  are  within  the  margin  of 
experimental  error.  This  is  in  contrast  to  the  conclusions  drawn  in  the 
earlier  report  (Narayan  and  Plunkett,  1985)  where  a  substantial  interaction 
was  found  because  of  the  very  short  transverse  correlation  length  at  zero 
axial  distance.  We  do  not  yet  have  a  reasonable  explanation  for  this 
discrepancy  but  assume  it  is  caused  by  the  slight  steps  on  the  edges  of  the 
rotating  holders. 
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Figure  V3.2.  Plate  vibration  response  autospectrum. 
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Current  Status 


All  of  the  experimental  work  has  been  completed.  All  of  the  auto  and 
cross  spectra  for  both  pressure  and  acceleration  have  been  reduced.  The 
predicted  vibration  cross  spectra  are  being  calculated  from  the  dimen¬ 
sionless  pressure  spectra  and  the  measured  plate  response  functions.  When 
this  is  completed,  the  results  will  be  compared  with  the  measured  vibration 
cross  spectra. 

The  thesis  Is  being  written  at  the  same  time.  When  that  Is  completed, 
the  results  will  be  put  Into  the  form  of  published  papers  in  the  open 
literature. 
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4.  Tip  Vortex  Cavitation  and  lfolee 

Roger  E.A.  Arndt  &  H.  Higuchi  Principal  Investigator 


Introduction 


The  tip  vortex  research  program  was  broken  down  into  three  basic  com¬ 
ponents:  a)  measurement  of  cavitation  Inception  and  fully  developed  cavi¬ 
tation  characteristics,  b)  measurement  of  noise  due  to  cavitation,  and 
c)  LDV  measurement  of  the  velocity  field,  construction  of  a  seml-empirlcal 
model  of  a  trailing  vortex  and  estimation  of  the  pressure  within  the  core. 
Four  different  facilities  were  used  in  the  experimental  program. 
Cavitation  and  noise  were  studied  In  a  low  noise  recirculating  water 
tunnel,  velocity  measurements  were  made  in  a  second  water  tunnel,  dye 
studies  were  carried  out  In  a  flume,  and  oil  film  studies  were  carried  out 
in  a  wind  tunnel.  The  objectives  of  the  study  were:  to  determine  the 
vortex  roll-up  process  and  the  minimum  pressure  in  the  vortex,  to  relate 
the  measured  minimum  pressure  to  cavitation  observations  and,  finally,  to 
determine  the  amplitude  and  spectral  characteristics  of  the  noise  signal 
and  relate  it  to  the  fluid  mechanical  aspects  of  tip  vortices. 

Experimental  Procedures 

Cavitation  inception  and  noise  were  studied  in  a  high  speed  water  tun¬ 
nel  having  a  test  section  with  a  190  mm  square  cross  section,  1250  mm  in 
length.  The  experimental  set  up  is  shown  schematically  In  Fig.  V4.1.  The 
top  wall  of  the  test  section  was  adjusted  to  produce  zero  axial  pressure 
gradient. 

Hydrofoils  having  an  elliptic  planform  and  an  aspect  ratio  of  3.0  were 
the  test  bodies.  Two  different  cross  sections  were  studied.  The  bulk  of 
the  investigation  utilized  a  hydrofoil  with  a  NACA  662~415  cross  section. 
This  shape  was  selected  to  produce  tip  vortex  cavitation  in  the  absence  of 
surface  cavitation.  A  second  hydrofoil  having  a  66-012  cross  section  was 
utilized  to  study  surface  cavitation  in  the  absence  of  tip  vortex  cavita¬ 
tion.  A  comparison  of  the  cavitation  characteristics  of  the  two  hydrofoils 
is  shown  In  Figure  V4.2. 

The  acoustic  signature  was  monitored  by  four  Bruel  &  Kjaer  (B  &  K) 
model  8103  hydrophones  (nominal  sensitivity  -211dB  re  1  V/Pa)  positioned  in 
a  stagnant  water  chamber  separated  from  the  test  section  by  an  almost 
acoustically  transparent  wall.  The  positioning  of  the  hydrophones  is  also 
shown  in  Fig.  V4.1.  The  signals  from  the  individual  hydrophones  were 
simultaneously  digitized  at  2  MHz  per  channel  by  a  Nicolet  model  4094  four- 
channel  digital  oscilloscope  with  a  disk  drive.  The  data  were  stored  on 
disks  and  were  analyzed  in  an  on-lim  HP9F36A  micro  computer.  The  digital 
oscilloscope  was  set  to  capture  the  single  noise  events  corresponding  to 
either  cavitation  Inception  or  continuous  cavitation  noise.  An  HP  digital 
wave  form  analysis  program  which  includes  a  standard  fast  Fourier  transform 


37 


Side  View 
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Top  View 


Figure  V4 . I .  Experimental  set  up 


routine  and  a  Levlngson  autoregressive  technique  was  used  with  minor  modi¬ 
fications. 

Cavitation  Inception  was  detected  both  visually  and  acoustically.  The 
Inception  was  marked  by  a  distinct  burst  of  noise  In  the  flow  field.  By 
pre-trlggerlng  the  A/D  converter  with  the  hydrophone  signal  itself,  Indivi¬ 
dual  traces  of  instantaneous  pressure  time  traces  at  four  different  loca¬ 
tions  could  be  recorded.  By  measuring  the  time  delay  of  arrival  among  the 
four  hydrophones,  the  sound  location  can  be  estimated  by  a  triangulatlon 
method.  Time  delays  among  the  Initial  pulses  on  the  traces  were  measured 
with  a  cursor  on  the  digital  oscilloscope  to  an  accuracy  of  0.5  micro 
second.  A  cross-correlation  technique  applied  to  the  pressure  traces  of 
the  initial  regions  produced  the  same  time  delay  as  that  obtained  directly 
from  the  traces.  Given  three  time  delays  among  four  hydrophones  as  inputs, 
implicit  simultaneous  equations  for  the  three  coordinates  of  the  sound 
source  were  solved  iteratively. 

The  pressure  traces  indicate  that  a  significant  amount  of  rever¬ 
beration  is  present.  By  limiting  the  source  location  analysis  to  the 
direct  signals  from  single  events  prior  to  the  reflection  from  the  tunnel 
walls,  the  effects  of  reverberation  can  be  treated.  The  present  method  of 
sound  source  identification,  unfortunately,  is  not  applicable  to  the  noise 
signal  from  more  developed  cavitation. 

The  four  hydrophones  were  placed  so  the  separation  from  each  other  was 
maximized  with  the  constraint  that  direct  radiation  from  the  noise  sources 
arrives  at  each  hydrophone  in  advance  of  the  arrival  of  any  reflected 
pulses. 

Acoustic  data  collection  was  concurrent  with  video  observation. 
Thus,  various  types  of  noise  signals  could  be  correlated  with  different 
types  of  cavitation  phenomena. 

A  very  wide  set  of  test  cases  were  studied,  as  shown  in  Fig.  V4.3. 
This  was  necessary  in  order  to  sort  out  the  very  complex  interrelationships 
between  viscous  effects,  dissolved  gas  content  and  free  gas  content  (small 
bubbles  or  nuclei). 

Velocity  measurements  were  made  in  water  using  laser  Doppler  velocl- 
metry.  Because  of  the  very  small  core  radius  at  the  higher  Reynolds  number 
tested,  the  optics  had  to  be  modified  to  improve  the  spatial  resolution. 
The  bulk  of  these  studies  were  made  in  a  second  recirculating  water  tunnel, 
having  the  same  size  test  section,  but  specially  modified  with  optical 
quality  glass  walls. 

Extensive  flow  visualization  was  also  carried  out  in  both  air  and 
water.  Surface  oil-flow  visualizations  were  conducted  in  air  to  study  the 
behavior  of  the  foil  boundary  layer.  These  visualizations  were  carried  out 
in  a  low  speed  open-return  wind  tunnel  with  a  43.2  cm  by  30.5  ,cm  test  sec¬ 
tion,  and  free  stream  velocities  which  could  be  varied  from  18  m/s  to  70 
m/s.  Two  Reynolds  numbers  were  tested:  Re  *  5.3xl05  and  27xl05 ,  based  on 
the  root  chord  (129.4  mm).  A  mixture  of  a  light  mineral  oil  and  titanium 
dioxide  was  used  as  a  flow  tracer. 


The  experiments  in  water  were  carried  out  in  a  30.5  cm-wide  water 
channel.  The  dye  visualizations  were  conducted  using  multiple  dye  injec¬ 
tion  holes  on  the  foil  and  observing  the  streakline  development  downstream. 

Cavitation  Studies 

The  variation  of  inception  index  with  angle  of  attack  for  various 
Reynolds  numbers  is  shown  in  Fig.  V4.4  for  a  gas  content  value  of  4  ppm. 
The  angle  of  attack  dependency  was  found  to  be  consistent  with  other 
investigators.  With  respect  to  the  Reynolds  number  dependence,  this  figure 
indicates  two  distinct  zones,  one  at  low  angle  of  attack  and  another  at 
high  angle  of  attack.  The  Reynolds  number  effect  is  seen  to  be  more  promi¬ 
nent  at  high  angle  of  attack.  This  point  is  more  clearly  shown  in  Fig. 

V4.5  where  the  inception  data  are  presented  as  a  function  of  Reynolds 
number  for  fixed  values  of  angle  of  attack  with  the  same  gas  content.  The 
figure  shows  that  for  angles  of  attack  less  than  7.5  degrees,  very  little 
variation  of  a.  with  Reynolds  number  is  observed,  while  for  higher  angles 
of  attack,  the  effect  of  the  Reynolds  number  is  more  prominent.  The  slope 
of  the  lines  appears  to  increase  monotonically  with  Increasing  angle  of 
attack.  It  may  be  pointed  out  here  that  details  of  the  flow  visualization 
study  showed  various  viscous  flow  regimes  at  different  angles  of  attack,  In 
particular,  separated  flow  existed  slightly  downstream  of  the  mid-chord 
portion  at  low  angles  of  attack.  However,  at  higher  angles  of  attack, 
separation  with  reattachment  shifted  to  the  leading  edge  resulting  in 
complex  flow  patterns.  The  latter  studies  were  limited  to  Reynolds  number 
of  5.7x10  ,  but  similar  effects  are  likely  to  exist  at  higher  Reynolds 
number  at  least  up  to  1  million.  It  is  very  likely  that  these  observations 
on  viscous  effects  are  very  relevant  to  the  observed  dependency  of 
values  on  Reynolds  number  and  angle  of  attack. 

Next,  additional  data  similar  to  those  presented  in  Fig.  V4.4  are 
shown  in  Figs.  V4.6  and  V4.7  at  Cg  values  of  8  ppm  and  12.5  ppm,  respec¬ 
tively.  The  results  for  8  ppm  appear  very  similr  to  those  for  4  ppm; 
however,  there  is  a  distinct  difference  In  the  trei~  it  Cg  ■  12.5  ppm.  The 
observed  strong  Reynolds  number  dependence  at  hi  .  angle  of  attack  for  4 
and  8  ppm  are  greatly  reduced  at  12.5  ppm.  It  i~  speculated  that  this  may 
be  related  to  the  nuclei  supply  conditions  which  are  Influenced  by  many 
factors,  some  of  these  being  facility  dependent  and  others  associated  with 
viscous  flow  characteristics  of  the  foil  itself.  However,  the  present 
results  bring  out  the  complex  interactions  between  various  physical  factors 
which  Influence  the  inception  process.  Some  of  these  questions,  to  the 
author’s  best  knowledge,  have  not  been  addressed  before.  In  order  to 
further  clarify  the  gas  content  effect,  a  gas  content  index  was  defined  as 
follows. 

G  -  Cg8  /(0.  SpU2  ) 


where  6  is  Henry's  constant.  Figure  V4.8  presents  a  set  of  results  using 
this  index  rather  than  the  absolute  value  of  gas  content  in  ppm.  Firstly, 
the  figure  shows  that  the  equilibrium  theory  for  correcting  gas  osntent 
effect  is  not  valid  since  the  slopes  do  not  remain  constant,  which  would  be 
required  for  the  equilibrium  theory.  The  slope  apparently  increases  with 


Figure  V4  .4 .  Cavitation  inception  index  as  a  function  of 
angle  of  attack  gas  content  4  ppm. 


Figure  V4.5.  Cavitation  inception  index  as  a  function  of 
Reynolds  number. 
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increasing  Reynolds  number,  which  in  this  case  was  varied  by  increasing 
velocity*  This  is  curious  since  the  residence  time  decreases  with 
increasing  velocity  which  would  have  indicated  the  reverse  trend.  Any 
further  conclusion  should  be  based  only  on  the  observations  carried  out 
with  changing  of  Reynolds  number  by  altering  the  size  of  the  foil. 

It  is  well  known  that  the  nuclei  characteristics  of  the  facility  can 
be  altered  by  subjecting  the  tunnel  water  to  various  initial  pressures 
prior  to  inception  testing.  The  effect  of  this  procedure  was  examined  in 
Fig.  V4.9.  It  is  clear  that  there  is  a  cut-off  pressure  below  which  the 
inception  index  is  strongly  influenced  by  pre-pressurization  history.  In 
other  cases,  very  similar  results  were  observed,  though  the  cut-off 
pressure  was  found  to  be  dependent  on  angles  of  attack  and  Reynolds  num¬ 
bers.  It  was  concluded  from  this  study  that  the  trend  of  the  dependence  of 
on  other  parameters  was  consistent  irrespective  of  the  test  procedure 
used  for  inception  observations.  However,  the  absolute  values  of  were 
found  to  be  dependent  on  this  procedure  as  is  evident  from  Fig.  V4.9. 

A  detailed  analysis  of  the  effects  of  nuclei  supply  was  made  and 
recently  reported  (Arakeri  et  al.,  1986).  Nuclei  can  be  supplied  to  the 
tip  vortex  from  the  free  stream  or  from  "body"  nuclei.  These  are  nuclei 
that  are  generated  within  separated,  transitional  or  turbulent  boundary 
layer  regions  on  the  hydrofoil  itself  (Katz,  1984).  It  was  found  that  the 
conditions  for  the  supply  of  body  nuclei  are  strongly  dependent  on  the 
viscous  flow  characteristics  of  the  body.  Boundary  layer  computations  were 
made  on  both  the  suction  and  pressure  side  of  the  particular  hydrofoil 
under  consideration.  The  theoretical  pressure  distribution  at  a  given 
angle  of  attack  was  computed  using  the  method  described  in  Abbott  and  Von 
Doenhoff  (1959).  These  do  not  Include  any  three-dimensional  effects  which 
are  likely  to  be  important  in  the  tip  region.  Next,  the  boundary  layer 
growth  calculations  were  made  using  the  Thwaites  method,  and  the  separation 
bubble,  if  existing,  was  decided  to  be  either  short  or  long  based  on  the 
results  of  Gaster  (1967).  One  set  of  computed  results  on  the  suction  side 
at  at  angle  of  attack  of  a-a  =  5.5°  is  shown  in  Fig.  V4.10.  A  laminar 

separation  bubble  is  predicted  as  indicated.  Since  on  an  elliptic  planform 

the  local  Reynolds  number  decreases  sharply  towards  the  tip,  a  long  separa¬ 
tion  bubble  is  predicted  near  the  region.  The  point  of  changeover  from 
short  to  long  separation  bubble  shifts  towards  the  tip  as  Reynolds  number 
is  increased,  Indicated  by  a  small  arrow  mark.  Only  at  a  Reynolds  number 
of  10  ,  a  short  separation  bubble  is  predicted  to  exist  almost  all  the  way 
to  the  tip.  Approximate  stability  estimates  were  also  made  to  predict  the 
transition  location  within  the  attached  boundary  layer.  Separation  is  pre¬ 
dicted  just  downstream  of  this  at  about  63  percent  chord.  Using  the  maxi¬ 
mum  spatial  amplification  factor,  it  was  predicted  that  laminar  separation 

would  prevail  up  to  a  Reynolds  number  of  2x10° . 

On  this  foil,  separation  is  also  predicted  on  the  pressure  side  at  the 
angles  of  attack  of  present  interest.  In  addition,  the  location  of  separa¬ 
tion  and  pressure  coefficient  at  separation  were  not  found  to  be  a  strong 
function  of  angle  of  attack  as  indicated  below  for  the  suction  side. 


Figure  V4.I0.  Dependence  of  a.  on  the  initial  cavitation 
number  at  the  stari  of  the  inception  tank. 


Figure  V4.ll.  Conditions  used  for  estimating  dominant 
nuclei  size  and  predicted  behavior  of 

a .  •  y  *  4  ppm . 
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Assuming  that  separated  regions  become  sources  of  nuclei  when  the  flow 
locally  becomes  saturated  with  gas,  the  critical  values  of  a  for  this  to 
occur  on  the  suction  side  or  the  pressure  side  can  be  calculated  as 
°crl  an<*  CTcr2*  acr3  will  be  used  to  denote  supersaturation  in  the  test 
section.  Close  examination  of  the  inception  data  shows  certain  interesting 
general  trends  when  compared  with  the  behavior  of  critical  a's  for  nuclei 
supply.  One  of  the  clear  trends  is  that  if  a  >  ocr ^  ,  then  03  is  a  strong 
function  of  Re  with  a 3  increasing  with  an  increase  in  Re.  The  second  trend 
which  is  also  apparent  is  that  in  the  Reynolds  number  range  for  which 
°cr3  <  ai<CTcr2»  CTi  ®lmo«t  independent  of  Re  and  finally  we  note  that  if 
Of  <  acr^,  then  03  increases  with  a  decrease  in  Re.  It  is  not  obvious  here 
that  this  trend  is  due  to  gaseous  cavitation  normally  observed  since  we  do 
not  know  the  precise  magnitude  of  -Cp^n*  Using  these  ideas,  along  with  an 
estimation  of  dominant  nuclei  size,  as  discussed  in  Arakeri  et  al.  (1986), 
the  salient  features  of  the  cavitation  inception  process  can  be  predicted. 
An  example  is  shown  in  Fig.  V4.ll.  In  this  case,.  R*,  -  6.3xl0~&  (1.9  pm) 
and  -Cpn^n  -  3  values  were  used.  For  Re  >  6.4x10  ,  the  observed  03  depen¬ 
dence  with  Re  can  be  explained  solely  on  the  basis  of  tension  effects.  For 
this  Reynolds  number  range  also,  there  is  no  obvious  source  of  nuclei  based 
on  static  pressure  considerations,  and  hence  the  assumption  that  Rq  remains 
constant  may  be  a  valid  one. 

For  Re  <  6.4  x  105 ,  the  observed  03  follows  a  different  trend  from  the 
predicted  one.  In  this  region,  at  inception,  the  extent  of  supersaturation 
being  proportional  to  ocri  -  03  Increases  with  decrease  in  Re  or  velocity. 
Thus  we  may  not  expect  Rq  to  remain  constant  here  but  to  increase  roughly 
in  proportion  to  (ocr^  -  03)  ~  1/1/2  pU2  .  In  the  same  Reynolds  number 
range  A  a  required  for  inception  increases  like  l/2plr  and  the  two  effects 
may  cancel  giving  a  constant  03  value.  This  is  roughly  in  agreement  with 
the  observations. 

In  summary,  there  is  definite  evidence  that  the  observed  complex 
behavior  of  03  for  tip  vortex  cavitation  inception  with  Reynolds  number 
(actually  velocity  since  body  size  was  not  varied),  angle  of  attack  and 
dissolved  gas  content  are  related  to  the  critical  conditions  for  nuclei 
supply .  It  is  necessary  to  examine  the  latter  for  both  the  free  stream  and 
separated  regions  that  may  exist  on  the  foil.  In  the  regions  of  Reynolds 
numbers  (velocities)  where  there  is  no  obvious  source  for  nuclei  supply, 
the  03  value  is  a  strong  function  of  Reynolds  number  and  is  almost  indepen¬ 
dent  of  the  dissolved  gas  content  value.  From  the  measured  03  values,  it 


was  possible  to  estimate  the  dominant  nuclei  size  responsible  for  inception 
and  these  are  found  to  be  surprisingly  small  being  of  the  order  of  2  pm. 
Thus,  some  of  the  Reynolds  number  dependence  on  oj  can  be  ascribed  to  the 
tension  effects.  It  appears  that  in  dealing  with  inception  characteristics 
of  flows  with  high  -Cproin  values,  it  is  necessary  to  carefully  examine  the 
nuclei  supply  conditions.  For  flows  with  lower  -Cp,,,^  values  like  axisym- 
metric  bodies,  the  nuclei  supply  conditions  may  not  be  as  critical.  The 
above  conclusions  are  valid  for  inception  conditions. 


Vortex  Trajector  as  Visualized  by  Cavitation 

The  tip  vortex  roll-up  process  and  the  vortex  trajectory  are  of  prac¬ 
tical  Importance  as  evidenced  by  a  substantial  number  of  recent  numerical 
and  experimental  studies.  In  the  present  investigation,  cavitation  pro¬ 
vided  an  excellent  method  for  visualizing  the  mechanism  of  vortex  roll-up. 

Stroboscopic  lighting  with  a  5  micro  second  duration  was  used.  First 
the  velocity  of  the  incoming  flow  was  set,  and  the  static  pressure  of  the 
test  section  was  gradually  lowered.  Cavitation  first  occurs  in  the  low 
pressure  core  region,  and  the  inception  of  tip  vortex  cavitation  is 
observed  to  be  the  Intermittent  appearance  of  a  continuous  cavity  starting 
within  one  or  two  cord  lengths  from  the  tip  of  the  foil.  As  the  pressure 
is  further  decreased,  a  developed  tip  vortex  cavity  formed.  While  the 
cavity  core  radius  is  small,  the  flow  field  about  the  cavity  can  be  con¬ 
sidered  to  be  representative  of  the  single  phase  flow.  Numerous  pho¬ 
tographs  of  the  tip  vortex  cavity  were  taken  at  different  angles  of  attack, 
Reynolds  number,  and  gas  content.  These  photographs  were  analyzed  and  the 
curvature  and  the  spanwlse  position  of  the  cavity  were  studied.  The  vortex 
trajectory  up  to  5.8  chord  lengths  was  measured  and  at  4  chord  lengths,  the 
spanwlse  position  reaches  within  2%  of  the  theoretical  asymptotic  position 
for  an  elliptically  loaded  wing  (y/s  “  0.785)  The  trajectory  agrees  well 
with  the  experimental  data  of  Spreiter  and  Sachs  (1951)  for  a  larger  aspect 
ratio  (AR»5)  than  the  present  foil.  Figure  V4.12  summarizes  the  data  In 
the  vicinity  of  the  foil  taken  at  various  test  conditions.  Lack  of  corre¬ 
lation  with  gas  content  is  as  expected,  but  it  is  to  be  noted  that  the 
general  be  lavior  of  the  trajectory  is  almost  independent  of  the  Reynolds 
number  and  the  angle  of  attack.  Spreiter  and  Sachs  (1951)  applied  Kaden's 
theory  in  the  near  field  to  estimate  the  rate  at  which  the  vortex  trajec¬ 
tory  approaches  its  asymptotic  value.  According  to  this  analysis,  the 
downstream  distance  at  which  the  vortex  trajectory  crosses  the  asymptotic 
value  is  expected  to  be  inversely  proportional  to  the  lift  coefficient, 
however,  the  variation  among  the  present  experimental  data  was  found  to  be 
significantly  smaller  than  predicted.  For  comparison,  the  vortex  trajec¬ 
tory  was  also  calculated  using  the  Betz'  (1913)  roll-up  model  and  Moore's 
(1974)  numerical  solution.  The  computed  approach  of  the  trajectory  to  the 
asymptote  lagged  experimental  data  considerably,  and  the  agreement  was 
poorer  than  Kaden's  model.  Naturally,  these  classical  models  are  too 
simplistic  since  the  flow  over  the  wing  or  the  three-dimensional  roll-up 
process  is  not  taken  into  account.  Recently,  Krasny  (1987)  carried  out  a 
computation  of  vortex  sheet  roll-up  using  a  large  scale  discrete  vortex 
method.  The  results  for  an  elliptically  loaded  wing  agreed  with  the 
Kaden's  analysis  in  the  near  field,  and  also  agreed  fairly  well  with  the 
experiment  on  the  asymptotic  vortex  position.  However,  the  transformation 
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from  the  Trefftz  plane  to  the  physical  plane  is  linearly  related  to  the 
wing  circulation,  and  the  result  becomes  more  sensitive  to  the  angle  of 
attack  than  experimentally  observed.  Numerical  simulation  with  the  full 
Navier-Stokes  equations  of  the  vortex  roll-up  from  a  rectangular  wing  has 
been  reported  by  Srinlvasan  et  al.  (1986),  but  no  numerical  data  exist  on 
the  low  aspect  ratio  elliptic  wing. 


The  Structure  of  Trailing  Vortices 

In  order  to  relate  the  observed  cavitation  characteristics  to  the 
minimum  pressure  within  the  tip  vortex,  a  semi-empirical  model  of  the  near 
wake  structure  was  developed.  This  model  was  based  on  LDV  measurements 
made  in  a  water  tunnel  and  a  flume. 

The  structure  of  vortices  trailing  from  wings  of  finite  span  has  been 
treated  extensively  by  previous  investigators.  The  impetus  for  this  pre¬ 
vious  work  was  aircraft  safety.  The  results  of  this  research  do  not  pro¬ 
vide  sufficient  information  for  the  tip  vortex  cavitation  problem 
(McCormick,  1962;  Platzer  and  Souders,  1979).  The  minimum  pressure  in  the 
vortex  is  not  of  much  interest  in  aerodynamics,  since  the  aircraft  hazard 
problem  is  concerned  with  the  structure  of  the  far  wake  and  not  the  near 
wake  where  cavitation  inception  occurs.  In  much  of  the  previous  work  on 
vortex  cavitation,  the  data  are  analyzed  on  the  basis  of  an  assumed 
Rankine  vortex  model  (Billet  and  Holl,  1979;  Latorre  et  al.,  1984).  The 
core  radius  in  these  studies  was  assumed  to  scale  with  the  boundary  layer 
thickness  on  the  pressure  side  of  the  hydrofoil.  Recently,  it  has  been 
pointed  out  that  there  are  inconsistencies  in  this  approach  (Arndt  et  al., 
1985). 

In  order  to  deal  with  the  trailing  vortex  analytically,  many  inviscld 
techniques  (Bilanin  and  Donaldson,  1975)  have  been  developed.  These  models 
have  been  expanded  to  include  laminar  and  turbulent  diffusion  effects 
(Lamb,  1945;  Moore  and  Saffman,  1973;  Phillips,  1981;  Squire,  1965) 
However,  none  of  the  available  theoretical  models  are  successful  in  accu¬ 
rately  describing  the  near  wake  characteristics  of  a  trailing  vortex.  It 
was  found,  for  example,  that  the  functional  dependence  of  core  radius  ok 
both  Reynolds  number  and  the  angle  of  attack,  which  has  been  confirmed  by 
the  recent  LDV  measurements  (Orloff,  1974;  Quadrelli,  1985),  is  not  pre¬ 
dicted  in  any  of  the  existing  analyses.  Moreover,  the  predicted  core 

pressure  diverges  near  the  tip. 


The  evolution  of  a  trailing  vortex  with  downstream  distance,  x,  may  be 
broken  down  into  four  basic  regions  as  illustrated  in  Fig.  V4.13.  In  the 
fully  rolled-up  region,  the  roll-up  process  is  complete  and  observed  self¬ 
similar  characteristics  are  consistent  with  the  previously  mentioned  analy¬ 
tical  work.  In  the  initial  roll-up  region  of  interest,  more  complicated 
near  wake  characteristics  are  evident.  For  example,  the  vortex  circulation 
I”*  is  not  necessarily  equal  to  the  mid-span  circulation  r0«  In  addition, 
the  tangential  velocity  profiles  are  asymmetric.  The  asymmetry  of  the  flow 
is  difficult  to  address  and  has  been  ignored  in  this  simple  analysis.  An 
important  feature  of  the  model  is  the  consideration  of  turbulent  diffusion 
in  the  core  region.  The  modelled  "inner"  and  "outer”  core  regions  are 
similar  in  structure  to  Squire's  (1965)  turbulent  vortex.  With  the 
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Figure  V4.13.  Evolution  of  trailing  vortex 


modelled  velocity  field,  the  pressure  drop  In  the  vortex  core  can  be  esti¬ 
mated  by  integrating  the  momentum  equation  for  an  axisymmetrlc  swirling 
flow 

-  »vVr 

The  model  depends  heavily  on  measured  velocity  data.  In  addition  to 
the  measurements  reported  in  Quadrelli  (1985),  extensive  water  tunnel 
measurements  were  made  at  a  Reynolds  number  of  5.6x10  .  A  comparison  of 
the  model  and  measurements  is  shown  in  Figures  V4.14  and  V4.15.  Details 
are  found  in  Ikohagi  et  al.  (1986).  Figure  V4.14  compares  the  measured 
circulation  as  a  function  of  distance  from  the  centerline  of  the  vortex. 
The  agreement  with  the  model  is  much  better  on  the  suction  side  of  the  vor¬ 
tex  (as  shown)  than  on  the  pressure  side.  In  fact,  an  axisymmetrlc  model 
is  only  a  fair  approximation  to  the  measured  velocity  data  as  shown  in 
Figs.  V4.l5a  and  V4.15b. 

Comparisons  of  tangential  velocity  are  made  for  various  a  and  x/C0 
to  examine  the  validity  of  the  proposed  semi-empirical  modelling.  Although 
there  is  some  asymmetry  of  the  velocity  profiles,  with  higher  velocities 
noted  on  the  suction  side  of  foil,  fair  agreement  is  obtained  except  for 
the  case  a  *  15°.  The  data  also  indicate  that  the  core  radius  increases 
with  increasing  a.  This  suggests  that  the  incipient  cavitation  number  for 
a  tip  vortex  scales  with  a  smaller  power  of  a  than  a2 .  The  axial  gradient 
of  pressure  is  almost  zero  within  a  few  chord  lengths  downstream  after  a 
8 harp  pressure  drop  near  the  tip.  The  semi-empirical  approach  used  here 
circumvents  the  problem  of  pressure  divergence  near  the  tip,  found  in 
existing  theoretical  models. 

Although  the  near  wake  characteristics  are  quite  different  from  the 
far  wake  •  haracteristics,  and  the  near  wake  can  be  adequately  modelled 
seml-empir  cally,  the  issue  of  what  physical  mechanisms  are  involved  still 
remains  tc  be  solved.  A  much  more  sophisticated  analytic  approach  is 
necessary  to  adequately  describe  the  near  wake  region  of  a  trailing  vor¬ 
tex,  at  least  for  the  purpose  of  further  elucidating  the  cavitation 
problem. 


Cavitation  Noise 

Different  types  of  cavitation  phenomena  were  observed  and  recorded  on 
still  photos  and  on  video  tapes.  Reynolds  number  (free  stream  velocity), 
angle  of  attack,  dissolved  gas  content  and  the  test  section  pressure  were 
varied. 

Surface  cavitation  on  the  symmetric  hydrofoil  was  as  expected.  The 
first  appearance  of  cavitation  was  along  the  leading  edge.  The  appearance 
of  the  cavitation  was  in  qualitative  agreement  with  the  assumption  that  the 
hydrofoil  was  elliptically  loaded,  though  a  more  careful  study  indicated 
three-dimensional  behavior  of  the  foil  boundary  layer  with  separation. 

The  Inception  of  tip  vortex  cavitation  occurred  either  as  an  unsteady 
appearance  of  a  single  bubble  within  one  or  two  chord  lengths  from  the  tip 


i£ 

3r 


53 


(rc/Ii  =  0.28) 


Figure  V4.14.  Measured  circulation  as  a  function  of  distance 
from  the  centerline  of  the  vortex. 


Figire  V4.I6.  Various  types  of  tip  vortex  cavitation. 
(;)  Inception  ( 0-1.4,  XiO=2.4  ,  Re=6.1xl0  )  Q 

(b)  Continuous  tip  vortex  cavity  (j=I.15,  :<-jO=9.5 

(c)  Tip  vortex  cavitation  with  surface  cavitation  ( 
Re=5.3xl05) . 


2=5.3x10  ) 

.43,  i-'tO-  11.5 
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or  as  a  continuous  cavity  suddenly  appearing  from  about  the  same  position 
and  attaching  itself  to  the  tip.  Typical  cavitation  phenomena  with  the 
cambered  foil  are  shown  in  Fig.  V4.16.  In  particular,  Figs.  V4.16b  and 
V4.l6c  were  taken  with  a  strobe  with  a  5  micro  second  duration.  Figure 
V4.16a  is  an  example  of  intermittent  cavitation  inception  where  a  single 
bubble  is  captured  on  the  film  in  the  vortex  roll-up  region.  In  other 
cases,  video  tapes  showed  that  the  inception  occurred  as  a  continuous 
cavity  travelling  upstream  and  attaching  itself  at  the  foil  tip.  An 
example  of  this  type  of  a  developed  tip  vortex  cavitation  is  shown  in 
Fig.  V4.16b.  In  either  case  of  inception,  distinct  sound  pulses  are 
radiated.  As  the  test  section  pressure  is  further  decreased,  the  cavity 
extends  over  the  entire  tip  region,  and  subsequently  the  vortex  develops  a 
twisted  ribbon  shape.  The  noise  level  decreases  with  the  appearance  of 
fully  developed  cavitation  until  a  cavity  starts  to  cover  a  portion  of  the 
foil  tip  with  further  lowering  of  the  cavitation  index.  As  the  pressure  is 
further  decreased,  surface  cavitation  develops.  As  shown  in  Fig.  V4.16c, 
cavitation  is  observed  in  the  region  where  trailing  edge  boundary  layer 
separation  was  predicted  and  observed  in  oil  film  visualization  studies. 

Three  types  of  pressure  traces  were  noted  during  tip  vortex  cavita¬ 
tion  in  inception.  Typical  inception  data  shown  in  Fig.  V4.17a  have  a 
sharp  pressure  peak  at  the  onset  followed  by  a  series  of  rebounds  due  to 
bubble  oscillation  and  tunnel  reverberation.  Similar  sharp  pulses  were 
observed  either  when  the  intermittent  short  cavitation  bubble  becomes 
visible  at  one  or  two  chord  lengths  downstream  in  the  tip  vortex,  or  when  a 
continuous  cavity  appears  at  a  similar  axial  position  and  rapidly  moves 
upstream  and  attaches  Itself  to  the  foil  tip.  In  other  cases,  very  regular 
low  frequency  oscillations  occurred  (see  Fig.  V4.17b).  These  low  frequency 
oscillations  were  observed  at  different  gas  contents,  Reynolds  numbers  and 
angles  of  attack.  The  frequency  of  oscillation  varied  from  2.5  KHz  to  6 
KHz  depending  on  various  test  conditions.  In  many  cases,  particularly  at 
low  angles  of  attack,  cavitation  inception  had  the  appearance  of  elongated 
bubbles,  presumably  the  result  of  discrete  portions  of  the  vortex  core 
intermittently  cavitating  as  shown  in  Fig.  V4.16a.  Although  the  natural 
frequency  of  such  elongated  bubbles  is  unknown,  the  observed  noise  fre¬ 
quency  range  is  comparable  to  the  natural  frequencies  of  spherical  bubbles 
of  equivalent  size.  One  problem  with  this  hypothesis  as  a  cause  of  the 
noise  is  that  this  type  of  noise  was  observed  over  a  wide  range  of  gas 
content.  Unfortunately  there  is  no  discernible  relationship  between  fre¬ 
quency  and  gas  content.  Further  analysis  of  a  cylindrical  cavity  is  needed 
to  more  accurately  describe  the  source  of  the  low  frequency  noise.  The 
third  type  of  noise  source  occurs  when  the  gas  content  is  low,  or  when  the 
cavitation  inde  c  is  lowered  rather  rapidly.  Under  such  circumstances  a 
steady  cavity  core  often  appears  suddenly.  Figure  V4.17  corresponds  to 
such  a  case,  vhich  displays  a  negative  transient  pressure  due  to  this 
cavity  formation. 

Once  continuous  tip  vortex  cavitation  occurs,  the  noise  level 
decreases  markedly,  as  shown  in  Fig.  V4.17c.  As  the  pressure  is  further 
decreased,  the  cavity  core  radius  increases  and  the  foil  tip  starts  to 
cavitate.  The  sound  level  significantly  increases  with  the  onset  of  a  sur¬ 
face  cavity. 

Within  a  very  narrow  band  of  cavitation  number  and  at  a  surface  cavity 
angle  of  attack  (»  10±2°),  a  discrete  tone  is  generated  in  the  frequency 


1  Volt 


Time  (m  sec) 


Figure  V4.17.  Time  traces  of  tip  vortex  cavitation  (  *12.5  ppm) 

(a)  Inception  (a=2.76,  a-a0=13.5°,  Re=6.1xloj?) 

(b)  Inception  (a=3.97,  a-aO=!2.5°,  Re=9.5xl(r) 

(c)  Continuous  tip  vortex  cavity 

(a=3.29,  a-aO= 15.5°,  Re=9.Ixl05) 

(d)  Tip  vortex  cavitation  with  surface  cavitation 

(cr=0.52,  a-a0=9.5°,  Re=7.5xl05) 


frequency  (KHz) 

Figure  V4.I8.  Power  spectral  density 
(a  &  b)  Inceptions 

(c)  Continuous  tip  vortex  cavity 

(d)  Tip  vortex  cavitation  with  surface  cavitation 

(same  as  in  Fig.  V4.I7) 


range  340-900  Hz.  Holding  free  stream  velocity  constant,  it  is  observed 
that  the  frequency  of  the  tone  varies  continuously  within  a  certain  range 
of  cavitation  number,  a  higher  pitch  corresponding  to  higher  cavitation 
number  and  smaller  cavity  core  radius.  Limited  data  indicate  that  the  fre¬ 
quency  varies  almost  linearly  with  velocity  at  constant  cavitation  number. 
This  phenomenon  is  apparently  associated  with  cavitation  on  the  suction 
side  of  the  tip  and  the  braided  structure  of  the  vortex  cavity.  Under 
strobscopic  Illumination,  it  is  observed  that  the  tip  cavity  pulsates  and 
the  braided  structure  rotates  in  synchronism  with  the  emitted  tone.  At 
first  it  was  conjectured  that  the  phenomenon  was  due  to  a  flow  Induced 
vibration  of  the  foil  since  the  resonance  frequency  of  the  foil  was  deter¬ 
mined  to  be  about  400  Hz.  TLe  fact  that  the  tonal  frequency  varied  con¬ 
tinuously  with  free  stream  velocity  and  cavitation  number,  discounts  this 
theory.  The  actual  mechanism  for  noise  generation  in  this  regime  remains 
to  be  investigated. 

As  the  cavitation  index  is  further  lowered,  surface  cavitation  appears 
in  the  vicinity  of  the  region  where  a  laminar  separation  and  a  turbulent 
reattachement  occurs.  This  region  has  been  observed  to  be  in  agreement 
with  theoretical  prediction  and  the  flow  visualization  studies.  The  noise 
level  increases,  mainly  in  the  high  frequency  range,  as  shown  in  Fig. 
V4.17d.  No  marked  effect  of  different  gas  content  (between  4.5  ppm  and 
12.5  ppm)  was  observed. 

The  effect  of  tunnel  reverberation  was  studied  for  the  cases  of  sur¬ 
face  cavitation  inception  as  reported  by  Higuchi  and  Arndt  (1985a).  By 
ensemble  averaging  the  cavitation  in< eption  data,  a  wave  form  of  a  single 
bubble  collapse  in  the  free  field  was  obtained.  A  simple  acoustic  model  of 
the  reverberant  water  channel  was  uti  ized  to  quantify  the  reverberation  in 
a  computer  simulation.  The  power  s>ectrum  of  a  simulated  single  bubble 
collapse  was  compared  to  the  power  spectrum  of  a  simulated  pulse  with 
reverberations  and  the  Fourier  transform  of  the  measured  pressure  trace. 
It  was  found  that  the  effects  of  reverberation  on  the  measured  spectrum  was 
only  on  the  relative  magnitude.  The  jeneral  shape  of  the  spectrum  remained 
unchanged  and  was  in  general  agreement  with  Fitzpatrick  and  Strasberg's 
(1956)  theoretical  spectrum  for  a  single  bubble. 

Arakeri  and  Iyer  (1986b)  conducted  a  study  of  the  noise  due  to  bubble 
induced  by  nuclei  seeding  by  an  electrolysis  technique.  They  found  that 
the  spectrum  of  the  noise  from  multiple  cavitation  collapses  retains  the 
shape  cf  the  single  bubble  spectrum.  The  amplitude  scales  with  the  rate  of 
cavitation  events.  They  also  noted  that  the  mutual  interaction  among 
bubbles  is  negligible.  This  is  consistent  with  our  present  computer  simu¬ 
lation  of  the  single  bubble  collapse  accompanied  by  reflected  waves.  The 
Important  point  is  that  cavitation  noise  can  be  studied  in  a  highly  rever¬ 
berant  water  tunnel. 

Acoustic  spectra  of  tip  vortex  cavitation  at  different  stages  of  deve¬ 
lopment  are  present  in  Fig.  V4.18.  As  mentioned  earlier,  two  types  of 
acoustic  signals  were  noted  during  tip  vortex  cavitation  inception.  The 
power  spectrum  corresponding  to  an  inception  with  a  sharp  sound  pulse  (Fig. 
V4.17a)  exhibited  a  gradual  decrease  of  spectrum  level  at  high  frequency. 
Unlike  the  spectra  for  surface  cavitation  inception,  no  high  frequency 
spectral  peak  is  noted.  The  power  spectrum  shown  in  Fig.  V4.18b 
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corresponds  to  inception  when  there  is  a  low  frequency  regular  oscillation 
(Pig.  V4.17b).  A  dominant  peak  is  evident  at  low  frequency,  and  the 
spectrum  drops  off  rather  rapidly.  The  noise  level  decreases  abruptly  as 
steady  tip  vortex  cavitation  develops  to  form  a  well  defined  cavity  core 
trailing  from  the  tip  of  the  foil  (e.g.  Fig.  V4.16b).  Figure  V4.18  shows 
the  corresponding  spectrum.  Within  a  narrow  range  of  the  cavitation  number 
where  the  ribbon  type  cavity  appears,  the  distinct  tone  is  generated  with  a 
sharp  spectral  peak  at  340-900  Hz.  The  pitch  decreases  and  then  the  noise 
level  decreases  again  at  lower  cavitation  numbers.  The  noise  level 
Increases  dramatically,  especially  in  the  higher  frequency  bends  when  sur¬ 
face  cavitation  appears,  as  shown  in  Fig.  V4.18d.  The  power  spectrum  of 
the  combined  tip  vortex  and  surface  cavitation  is  very  similar  to  the 
spectrum  of  continuous  surface  cavitation  alone  on  the  symmetric  foil. 
Thus  it  can  be  assumed  that,  when  it  occurs,  surface  cavitation  is  the 
dominant  noise  source. 

The  present  inception  data  under  various  conditions  were  examined  for 
their  source  locations  and  the  results  are  tabulated  in  Rogers’  thesis 
(1987).  Substantial  scatter  in  the  resulting  positions  (as  much  as  30  mm 
under  the  similar  conditions)  were  observed.  Due  to  the  small  vertical 
separation  between  the  hydrophone  No.  4  and  the  other  hydrophones,  the 
numerical  convergence  of  the  vertical  coordinate  was  often  poor.  In  these 
cases  the  vertical  position  of  the  visualized  vortex  was  given  as  alter¬ 
native  input.  The  latter  was  the  method  reported  earlier  when  only  three 
hydrophones  were  available  (Higuchi  et  al.,  1985b).  In  spite  of  the 
scatter  in  the  data,  the  present  results  clearly  showed  that  the  noise 
sources  are  concentrated  very  close  to  the  hydrofoil  tip.  They  were  either 
within  a  couple  of  chord  lengths  downstream  from  the  tip  corresponding  to 
the  vortex  roll-up  region,  or  in  the  Immediate  neighborhood  of  the  foil  tip 
corresponding  to  the  upstream  end  of  the  continuous  cavity.  These  obser¬ 
vations  were  in  good  agreement  with  the  visual  observation  recorded  on 
video  tape. 


The  findings  of  this  investigation  are  summarized  below: 

Several  different  types  of  cavitation  inception  noise  signatures  were 
noted.  Cavitation  inception  was  observed  to  be  either  the  sudden 
appearance  of  a  continuous  cavity  or  the  intermittent  appearance  of  an 
elongated  bubble  extending  axially  over  a  relatively  small  portion  of  the 
tip  vortex.  The  transient  noise  signal  due  to  the  former  type  of  inception 
consisting  initially  of  a  negative  pressure  pulse  implies  that  sudden 
expansion  of  a  cylindrical  cavity  was  important  in  the  noise  generation 
process.  In  the  latter  case,  either  high  frequency,  intermittent  pulses, 
or  a  relatively  low  frequency  damped  oscillation  was  observed  in  the  noise 
signature.  It  was  not  possible  to  relate  the  occurrence  of  either  type  of 
signal  with  Reynolds  number  or  gas  content,  nor  was  it  possible  to  relate 
the  noise  signals  to  gas  bubble  oscillation  (pseudo  cavitation)  or  to  the 
growth  and  collapse  of  true  vapor  cavities. 

A  hydrophone  array  was  utilized  to  identify  the  noise  source  location. 
In  spite  of  limitations  on  spatial  resolution  and  repeatability,  the 
results  of  this  study  clearly  indicate  that  cavitation  noise  sources  are 
concentrated  in  the  vortex  rollup  region,  less  than  one  chord  length  from 
the  tip  of  the  hydrofoil. 


The  noise  intensity  associated  with  fully  developed  tip  vortex  cavita¬ 
tion  is  substantially  lower  than  inception  noise,  or  surface  cavitation 
noise,  when  it  occurs.  At  high  angle  of  attack,  and  under  special  cir¬ 
cumstances,  a  very  strong,  low  pitch  tone  occurs.  The  frequency  of  the 
tone  varied  with  cavity  core  radius  and  free  stream  velocity,  discounting 
the  thought  that  this  might  be  a  flow  Induced  vibration  of  the  hydrofoil 
itself.  Although  no  conclusive  statement  could  be  made  from  the  present 
set  of  data,  it  was  speculated  that  the  oscillation  of  the  twisted  ribbon- 
type  cavity  is  responsible  for  this  tone. 

If  surface  cavitation  occurs,  it  is  the  dominant  noise  source, 
overwhelming  any  type  of  vortex  cavitation  noise.  Surface  cavitation  in 
the  absence  of  tip  vortex  cavitation  was  examined  with  a  hydrofoil  having  a 
symmetric  foil  section.  Well  defined  noise  pulses  were  observed,  which 
correspond  to  individual  bubble  implosions  on  the  hydrofoil  surface.  The 
noise  sources  were  concentrated  in  a  narrow  region  on  the  foil. 

Although  these  tests  were  carried  out  in  a  highly  reverberant  water 
tunnel,  experiments  and  computer  simulations  demonstrated  that  the  shape  of 
the  measured  spectra  was  unaffected. 
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The  Covariance- Constrained  Maximum  Likelihood  Method  (CCMLM) 
defines  a  class  of  power  spectrum  estimators  which  adjust  Capon's  Maximum 
Likelihood  Method  to  best  fit  the  observed  covariance  of  the  random  process.  In 
this  paper,  a  particular  form  of  the  CCMLM  is  introduced,  several  of  its  proper¬ 
ties  are  discussed,  and  examples  of  its  application  are  provided. 


1.  Introduction 

In  recent  years,  there  has  been  renewed  interest 
in  spectrum  estimation  based  on  non- parametric  pro¬ 
cedures  such  as  Capon's  Maximum  Likelihood  Method 
(MLM)  or  the  Blackman-Tukey  (BT)  method.  This 
interest  stems,  in  part,  from  the  observation  that  the 
performance  of  parametric  procedures  may  be  highly 
sensitive  to  errors  in  the  assumed  model  [1|.  Since 
MLM  and  BT  spectra  are  not  formulated  with  respect 
to  a  parametric  representation  of  the  spectrum,  these 

firocedures  have  proven  to  be  more  robust  to  model¬ 
ing  errors  than  parametric  methods.  However,  the 
loss  in  resolution,  or  fidelity,  that  is  associated  with 
these  methods  has  remained  one  of  their  primary 
disadvantages.  Significant  research  effort,  therefore, 
has  been  put  into  the  development  of  procedures  for 
improving  the  performance  of  these  non-parametric 
spectrum  estimators. 

Improvements  have  been  achieved  through 
"post-estimator’’  procedures  in  which  the  original 
estimated  spectrum  is  modified  to  satisfy  additional 
constraints.  These  constraints  are  most  often  stated 
with  respect  to  the  observed  data  or  their  statistics. 
Most  research  has  focussed  on  the  BT  spectrum  [2,3]. 
With  respect  to  the  MLM,  D’Assumpcao  |4]  has  pro¬ 
posed  an  "inverse-fit”  procedure  that  transforms  the 
MLM  spectrum  into  one  with  a  covariance  matrix  R' 
such  that  R‘ R'x  as  /,  where  R  is  the  observed  covari¬ 
ance  matrix. 

In  the  following  we  propose  a  new  method  for 
modifying  the  MLM  spectrum.  The  Covariance- 
Constrained  MLM  ICCMLM)  modifies  the  MLM  spec¬ 
trum  to  minimize  the  difference  between  the  observed 
covariance  matrix  and  that  of  the  spectrum  estimate. 
The  paper  is  divided  into  three  sections.  In  the  first 
section,  the  basic  approach  of  the  CCMLM  is 
presented  and  the  unity-gain  iterative  solution  to  the 
covariance-matching  problem  is  proposed.  The  second 
section  considers  briefly  issues  of  convergence  of  the 
unity-gain  procedure.  Finally,  the  third  section  com¬ 
pares  the  performance  of  the  unity-gain  CCMLM  with 
that  of  MLM  for  several  simulated  observation 


sequences. 

1.  Derivation  of  the  CCMLM 

1.1.  Definitions 

Let  P(u)  be  the  original  spectrum  with  covariance 
matrix  R  of  dimension  NxN.  We  shall  refer  to  PM 
as  the  MLM  spectrum  that  is  obtained  from  R.  Let 
the  covariance  matrix  of  PM  be  A,  and  the  error 
covariance  matrix  £,  be  defined  as 

*»“*»-*  (1) 

Since  the  MLM  spectrum  does  not  match  the  observed 
autocorrelation  sequence,  £,  will  be  nonzero  and  is,  in 
fact,  guaranteed  to  be  positive  definite  since 

/»,(«)  >  F(w)  [51.  Thus,  Ei  is  a  valid  covariance  matrix 
and  corresponds  to  the  error  spectrum  £(u)  defined  as: 

£(u)-P,(u)-P(w)  (2) 

3.2.  Transformation  of  Pju (w) 

We  are  interested  in  transforming  PA w)  (»-*-. 
f’xuufw))  into  P  (u)  such  that  R  ‘  -  R  —  0.  The  obvious 
transformation  is  to  choose 

P*(«)-P,(w)-£(w)  (3) 

which,  of  course,  cannot  be  achieved  since  we  observe 

the  truncated  error  covariance  matrix  £,.  Therefore, 
it  is  necessary  to  estimate  E[u)  from  £,  and  use  £,(u), 
the  estimated  error  spectrum,  to  modify  the  MLM 
spectrum. 

The  CCMLM  solves  the  problem  of  estimating 
the  error  spectrum  by  using  the  MLM.  Thus 


£,M  - 


x  *  M£  l'1*  (w) 


where  x(u>)  is  the  vector  of  Fourier  coefficients.  The 
original  spectrum  is  then  updated  by  computing  the 
difference 

PM  -  PM  -  EM  (5) 


A- 2 


Since  the  covariance  matrix  for  £,(w)  ia  not  necessarily 
equal  to  £ t, 

KtfiR  («) 

and  a  new,  non-xero  error  matrix  results: 

(7) 

The  results  above  suggsst  the  following  three-step 
procedure  for  iteratively  updating  the  CcMLM  spec¬ 
trum.  On  the  i-th  stage,  we  are  given  £,(w)  from 
which  we 

1.  construct  the  error  covariance  matrix 

£,-£,-£  (8) 


2.  compute  the  MLM  spectrum  for  the  matrix  E, 
E‘  ^  £*  (j)£,'‘x  (w) 


(9) 


3.  update  the  current  CCMLM  spectrum 

P.+M  -£,(«»-£(«)  (10) 

The  procedure  is  terminated  when  £,  is  sufficiently 
close  to  the  zero  matrix. 


the  following,  we  propose  a  suboptimal  solution  to 
this  problem  ia  which  weights  are  selected  indepen¬ 
dently  at  each  stage  of  the  iteration  so  that  the 
unity-gain  constraint  of  the  original  MLM  spectrum  is 
preserved.  Under  these  simplifying  assumptions,  it 
nas  been  shown  that  the  update  equations  depend  on 
a  single  scalar,  t,  which  must  be  selected  to  satisfy  the 
positivity  constraints  |6|.  Therefore,  the  update  equa¬ 
tions  for  the  unity-gain  CCMLM  are: 


(»*) 


and 

£.*-£,  (13) 

where  the  matrix  D,  is  the  covariance  matrix  for  the 
estimated  error  spectrum  at  the  i-th  stage,  £.(w).  The 
sequence  is  initialised  by 


f.twI-.-f  ..  .i 

nM*  i(w) 


(14) 


and 


£,-£,-£  (15) 

The  fi.'t  are  chosen  at  each  stage  such  that  £.(u>)  >  o 
and  £,*,  is  positive  definite.  In  the  examples  that  fol¬ 
low,  g,  was  selected  by  an  iterative  search  procedure. 


Although  other  spectrum  estimation  methods 
could  be  used  to  estimate  the  error  spectrum  at  each 
step,  the  MLM  is  chosen  for  several  reasons.  Firstly, 
the  MLM  is  often  used  in  those  cases  for  which 
parametric  methods  are  judged  inappropriate  because 
of  their  poor  approximation  or  statistical  properties. 
Therefore,  having  failed  to  estimate  accurately  the 
original  spectrum,  it  is  not  clear  that  a  parametric 
model  should  estimate  the  error  spectrum  with  any 
greater  fidelity.  Secondly,  among  nonparametric 
methods,  the  MLM  is  better  than  BT  since  it  gives  a* 
unbiased  minimum-power  estimate  at  each  frequency. 
Finally,  the  upper-bound  [5]  and  local  properties  [0]  of 
the  MLM  suggest  that  poor  fits  to  the  error  spectrum 
over  certain  regions  of  the  speetrum  will  be  confined 
primarily  to  neighborhoods  around  those  regions. 

The  upper-bound  property  of  the  MLM,  however, 
may  also  cause  problems  for  the  iterative  procedure. 
Let  us  consider  the  equation  for  the  first  update: 

PM  -P.M -£,(«>  (11) 

Since  £,(u)  >£H,  PM  may  become  negative  in 
regions  of  the  spectrum  where  the  difference  between 
EM  and  £(uf)  are  large  relative  to  PM-  In  such 
cases,  PM  “  no  longer  a  valid  power  spectrum  and 
the  updated  error  matrix  E,  is  not  necessarily  a 
covariance  matrix.  Since  the  MLM  is  defined  for 
covariance  matrices,  the  non-positive  matrix  E, 
presents  a  problem  at  the  next  stage  of  the  iteration. 
One  way  to  correct  this  problem  is  to  introduce 
weights  in  the  update  formulas  for  P,  (u)  and  t  so 
that  (i)  the  spectrum  is  always  non-negative  and  (si) 
the  error  matrix  is  always  positive  definite. 

2.3.  Unity-gain  procedure 

There  are  several  methods  for  incorporating  the 
positivity  constraints  into  the  update  equations,  la 


3.  Convergence  Properties. 

Several  commenta  with  respect  to  the  conver¬ 
gence  of  the  unity-gain  procedure  are  necessary  before 
proceeding  with  examples  of  the  CCMLM. 

Tr  £,*.  It  can  be  shown  that  there  will  always 
exist  a  non-tero  5,  at  each  stage  of  the  iteration  that 
satisfies  the  positive-definiteness  constraint.  For  such 
a  P- sequence,  the  trace  of  the  square  of  the  updated 
sequence.  Furthermore,  this  result  can  be  used  to 
show  that  the  updated  error  speetrum  is  stable,  in  the 
sense  that  is  does  not  become  unbounded  as  the 
number  of  iterations  becomes  large.  In  proving  these 
assertions,  however,  it  is  dear  that  the  unity-gain 
CCMLM  will  not,  in  general,  converge  to  a  solution 
for  which  the  error  covariance  matrix  goes  to  zero. 
Therefore,  the  unitv-gain  CCMLM  reduces  the  dispar¬ 
ity  between  the  observed  and  computed  covariance 
matrices,  but  does  not  necessarily  converge  to  an 
exact  match. 

Slopping  rriferie.  Stopping  criteria  have  been 
studied  with  respect  to  £  and  an  alternative,  £..,, 
which  is  formed  by  taking  the  difference  of  the  two 
covariance  matrices  after  each  is  normalised  with 
respect  to  variance  along  the  main  diagonal.  Our 
results  indicate  that  the  trace  of  £,*,  is  a  better  meas¬ 
ure  of  performance.  For  many  test  cases,  this  alterna¬ 
tive  measure  of  performance  exhibits  a  local  minimum 
which  c  jrrespond*  to  the  best  fitting  estimate  in  the 
spectral  domain.  This  result  is  consistent  with  the 
observe. ion  that  further  reductions  in  Tr  £'  reflect 
adjustments  ia  the  overall  displacement  of  the 
estimated  spectrum,  rather  than  refinements  in  its 
shape.  Therefort ,  we  suggest  that  a  stopping  criterion 
based  os  this  rormalisM  error  covariance  matrix  m 
preferable  to  one  based  os  the  non- normalised  form 


Seiuitivitg  to  tkt  /-itoutnet.  Based  oa  simula¬ 
tions.  the  unity-gain  CCMLM  does  not  appear  to  be 
sensitive  to  the  choice  of  ^-sequence;  different  9- 
sequences  appear  to  converge  to  the  same  spectrum. 
The  choice  of  0-sequence  appears  to  determine  only 
the  rate  of  convergence, 


4.  Examples 

Figures  1-4  compare  the  CCMLM  and  MLM 
methods  for  several  test  cases.  The  First  two  figures 
show  results  from  a  simulation  based  on  a  proposed 
model  of  the  wall-pressure  spectrum  beneath  a  tur¬ 
bulent  boundary  layer  using  a  linear  array  of  sensors. 
The  remaining  figures  show  results  from  simulations 
based  on  ARMA  spectra. 

4.1.  Wall-pressure  Wavenumber  S  pee  tram 

Considerable  interest  exists  among  fluid- 
dyn&micists  with  respect  to  the  wavenumber  spectrum 
of  the  wall-pressure  field  beneath  a  turbulent  boun¬ 
dary  layer.  Measurements  have  obtained  accurate 
estimates  of  this  spectrum  only  in  the  neighborhood  of 
the  convective  ridge  which  appears  as  a  broad  peak  in 
the  wavenumber  domain.  There  is,  however,  consider¬ 
able  controversy  over  the  nature  of  the  wavenumber 
spectrum  at  very  low  wavenumbers.  The  dashed  line 
in  Figures  1  and  2  represent  one  such  model  of  this 
wavenumber  spectrum  in  which  the  spectrum  is 
approximately  SO  dB  below  the  convective  peak  at 
wavenumbers  of  interest  [7]. 

The  non  homogeneous  nature  of  the  turbulent 
boundary  layer  restricts  the  aperture  of  the  sensor 
array  which  is  used  to  measure  the  wall  pressure  along 
the  boundary.  Thus,  this  wavenumber  spectrum  esti¬ 
mation  problem  is  similar  to  that  found  in  other  array 
signal  processing  problems  such  as  bearing  estimation 
in  which  the  number  of  sensors  is  typically  limited. 
However,  the  performance  criterion  differ*  from  that 
most  often  assumed  in  these  other  applications  since 
the  estimation  problem  is  focused  primarily  on  the 
general  shape  of  the  spectrum,  rather  than  on  the 
resolution  of  spectral  peaks.  We  are,  therefore, 
interested  in  a  high-fidelity  rather  than  higb- 
resolution  procedure. 

When  high-resolution  techniques  are  applied  to 
this  problem,  it  is  found  that  thev  perform  poorly  in 
tbe  regions  of  the  noise  floor  and,  therefore,  are  not 
appropriate  [8|.  Both  the  MLM  and  the  CCMLM,  on 
tbe  other  hand,  exhibit  considerably  better  Fidelity. 
The  results  of  4  simulations  based  on  12  sensors  and 
2S  snapshots  are  shown  by  the  solid  curves  for  the 
MLM  (Figure  1)  and  the  CCMLM  (Figure  2).  In  gen¬ 
eral,  the  bias  of  the  CCMLM  in  the  low  wavenumber 
region  of  tbe  spectrum  is  reduced  by  more  than  10  dB 
from  that  exhibited  by  the  MLM.  This  decrease  in 
bias  is  offset  somewhat  by  an  increase  in  tbe  variance 
of  the  estimate,  although  the  effecte  in  this  simulation 
are  not  great. 
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Fi*«r«  1 

SioelstioBS  of  tbe  MLM  (solid  hen)  for  the  Chase  model  of  the 
wall-pressure  wavesumber  spectrum  (dashed  lie*)  beaeath  a  tur¬ 
bulent  boaadsry  layer.  Simulations  arc  based  oa  25  saapsbota 
of  a  12-eeasor,  aaiform  tiaear  array.  Waveaamber  is  aormalised 
with  respect  to  tbe  eoavective  peak. 


Figare  2 

CCM.M  results  Sue  tbs  sum  conditions  ia  Fifaw  One. 

4.2.  ARMA  apnetraa  natlmation. 

One  of  the  intereating  properties  of  the  CCMLM 
is  that  it  is  a  pole/iero  form  of  estimator.  This  raises 
the  possibility  that  the  CCMLM  may  be  useful  in 
estimating  ARMA  spectra.  The  examples  shown  in 
Figures  3  and  4  are  based  on  25  simulations  of  the 
MLM  and  CCMLM.  The  dashed  line  represents  the 
actual  ARMA(2,2)  spectrum  (see  [0]  for  details  on  the 
spectra).  Each  simulation  was  based  on  40  samples  of 
toe  ARMA  process  and  estimates  of  tbe  First  12  terms 
in  tbe  sutocorrelation  sequence.  The  estimated  mean 
and  variability  (within  1  standard  deviation  of  the 
mean)  of  the  CCMLM  spectra  are  shown  by  the  solid 
and  dotted  lines,  respectively. 

The  results  support  the  assumption  that  the 
CCMLM  is  a  useful  procedure  for  estimating  ARMA 
spectra.  Several  trends  in  the  data  are  worth  noting, 
however.  Although  the  Fidelity  of  tbe  CCMLM  is 
superior  to  that  of  MLM,  the  variance  of  the  estima¬ 
tor  is  also  larger.  When  compared  with  that  for  other 
ARMA  methods,  however,  this  variance  is  not  unrea¬ 
sonable.  Secondly,  the  CCMLM,  unlike  other  A'  MA 
methods,  does  not  involve  an  axpUeit  computation  of 
model  order.  An  initial  ordtr  ia,  of  course,  implicit  in 
the  choice  of  dimension  for  the  covariance  matrix. 
This  selection  trades  bias  against  variance  and 
requires  further  study. 
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Figure  3 

MLM  (top)  »d  CCMLM  (bottos)  retires ire  for  ARMA(2.2) 
time  series  (dashed  liar)  breed  oa  40  suspire  red  estimates  of 
tbc  Tint  12  aotoeorrelatioe  lags.  The  (olid  bar  ia  racb  paaoi 
■bowa  tbc  estimated  meat  aad  the  dotted  bare  repreaeat  a 
range  of  1  ataadaid  deviatioa  about  tbc  mea*. 


Figure  4 

MLM  (top)  acd  CCMLM  (bottore)  for  a  differeut  ARMA(2,2) 
tine  aerica.  Coamtioa*  are  tbc  same  re  ia  Figure  Three. 


S  grantee  that  this  difference  goat  to  sero.  Neverthe- 
a,  the  fidelity  of  the  unity-gain  CCMLM  i>  superior 
to  that  of  the  MLM  and  other  parametric  procedure* 
when  evaluated  for  spectra  which  are  not  dominated 
by  multiple,  narrowly  spaced  peaks. 

Acknowledgments  This  work  was  supported  by  a 
Selective  Research  Opportunity  in  Hydroacoustics  by 
the  Office  of  Naval  Research. 


References 

1.  Tong,  P.  S.  aad  Rees,  F.  L.,  *TBL  Flow  Noise 
Measurements  Using  a  Maximum  Entropy  Method,* 
Binary  Systems  Inc.,  ONR  Contract  No.  N00014-72- 
C-0318,  April,  1879. 

2.  Konvalinka,  1.  S.,  'Iterative  nonparametric  spec* 
trum  estimation,"  IEEE  Transactions  on  Acoustics, 
Speech  and  Signal  Processing,  vol.  ASSP-32,  (1884), 
59-59. 

3.  Steinbardt,  A.  O.,  Goodrich,  R.  K-,  and  Roberts, 
R.,  "Spectral  estimation  via  minimum  energy  correla¬ 
tion  extension,”  IEEE  Transactions  on  Acouitics, 
Speech  and  Signal  Processing,  vol.  ASSP-33,  (1985), 
1509-1515. 

4.  d'Assumpcao,  R  A,  "Some  new  signal  processors 
for  arrays  of  sensors,*  IEEE  Transactions  on  Informa¬ 
tion  Theory,  vol.  IT-28,  (1980),  441-430. 

5.  Marietta,  T.  L.,  *A  new  interpretation  for  Capon’s 
maximum  likelihood  method  of  frequency- 
wavenumber  spectral  estimation,”  IEEE  Transactions 
on  Acoustics,  Speech,  and  Signal  Processing,  vol. 
ASSP-31,  (1983),  445-449. 

6.  Wakefield,  Gregory  H.,  High-fidelity  Estimation  of 
Continuous  Spectra,  Ph.D.  dissertation,  Department 
of  Electrical  Engineering,  University  of  Minnesota, 
Minneapolis,  MN,  1985. 

7.  Chase,  D.  M.,  ’Modeling  the  wavevector-frequency 
spectrum  of  turbulent  boundary  layer  wall  pressure,* 
Journal  of  Sound  and  Vibration,  v.  70,  (1980),  29-87. 

8.  Wakefield,  G.  H.  aad  Kaveb,  M.,  "Frequency- 
wavenumber  spectral  estimation  of  non-planar  ran¬ 
dom  fields,*  Proceedings  of  ICASSP-8t5,  Tampa, 
Florida  (1985),  806-811. 

9.  Kaveh,  M.  and  Bruisone,  S.  B.,  ’A  Comparative 
overview  of  ARMA  spectral  estimation,”  Proceedings 
of  the  1st  ASSP  Workshop  on  Spectra  Estimation, 
Hamilton,  Ont.,  Canada  (1981). 


6.  Conclusions 

The  CCMLM  is  a  general  procedure  which  modi¬ 
fies  the  MLM  spectrum  to  obtain  an  estimate  that  is 
more  consistent  with  the  observed  covarianee  matrix. 
An  iterative  suboptimal  implementation  of  this  pro¬ 
cedure  reduces  the  differences  between  the  observed 
aad  computed  covarianee  matrices,  but  does  not 
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ABSTRACT 

We  have  investigated  the  application  of  two  non¬ 
linear  data-adaptive  technique*  to  the  problaa  of 
estimating  the  frequency-wavenuatber  spectrua  of  the 
preeaure  field  beneath  a  Turbulent  Boundary  Layer. 
Although  the  Maximus  Likelihood  Method  (MLM)  exhibit* 
greater  bia*  at  wavenuabere  below  the  convective 
ridge,  Honte  Carlo  aiaulatione  demonstrate  that  it* 
variance  at  each  wavenuaber  i*  aaaller  than  that  of 
the  Haxlaua  Entropy  Method  and  i*,  therefore,  pre¬ 
ferred  in  caaee  of  Hatted  nuaber*  of  sensors  and 
•napshots  of  data.  Subapace  extension  of  the  HLH  1* 
shown  to  laprova  it*  perforaanc*  at  low-wavenuabara. 

INTRODUCTION 

The  problaa  of  estlaating  the  frequency- 
wavenuaber  (w-k)  spectrua  of  a  randoa  field  i*  a  topic 
of  great  lntereet  in  the  physical  science*.  While  a 
"large-sample"  theory  for  this  estlaation  problaa 
exists,  in  practice  aany  of  the  aasuaptlons  essential 
to  the  theory  are  violated.  Frequently,  a  randoa 
field  1*  stationary  and  hoaogeneoua  only  over  rela¬ 
tively  saall  spatial  separation*  and  teaporal  delay*. 
According  to  the  definition,  a  ee-fc  apactrua  for  such 
a  field  does  not  exist.  Nevertheless,  it  is  often 
useful  to  characterize  the  local  properties  of  the 
field  by  a  w~k  spectrua.  In  this  case,  on*  could 
apply  an  eatiaatlon  procedure  that  assustes  long  data 
records  end  aany  spatial  saaple*  to  a  problaa  with 
very  few  teaporal  and  spatial  saaplee  and  hope  that 
the  results  will  be  satisfactory.  Unfortunately,  they 
rarely  are. 

The  problaa  of  Halted  teaporal  and  spatial 
records  has  genarated  a  nuaber  of  alternative  spectral 
estimators.  Huch  of  the  original  work  in  this  area 
was  in  tlae-series  analysis  but  has  been  extended, 
over  the  last  twenty  years,  to  the  analysis  of 
tlaa/space  series,  particularly  with  respect  to  the 
bearing  estlaation  problea  in  radar  and  sonar. 


This  paper  la  as  auch  an  sffort  to  proaote 
cross-fertilisation  of  ideas  among  investigator*  i.i 


fluid  aschanics  and  spectral  eatiaatlon  as  it  is  an 
investigation  and  extension  of  spectral  aatiaatlon 
techniques  to  wall-pressure  spectra.  After  reviewing 
definitions  and  issues  concerned  with  sampling,  we 
present  a  brief  tutorial  on  issues  and  a* t hods  of 
Importance  in  spectral  estimation  theory.  He  avoid  a 
review  of  all  known  spectral  estimators  and  instead 
focus  attention  on  three  eetlaators  that  have  dom¬ 
inated  research  for  aany  years.  Thus,  this  section  is 
not  Intended  as  a  survey  of  the  available  technique* 
but  rather  as  an  introduction  to  the  general  problem 
of  spectral  estlaation  based  on  Halted  saaple*  of  a 
random  field.  We  than  dlscues  an  extension  of  the 
Haxlaua  Likelihood  Method  and  investigate  its  applica¬ 
tion  for  estlaating  wall-pressure  spectra. 

DEFINITIONS 

The  autocorrelation  function  (acf)  of  a  random 
field  is  def ined  aai 

■  El  pNtj.fcj)  P(t2.S2>  1  <1> 

where  p(t,$j  is  the  pressure  of  the  field  at  position 
x  and  at  ties  t,  superscript  **’  denotes  complex  con¬ 
jugal!  on  and  E[.)  is  the  statistical  expectation 
operator. 

For  a  stationary,  homogeneous  randoa  field,  the 
acf  is  invariant  with  respect  to  teaporal  lag  and  spa¬ 
tial  separation  and  can  be  written  asi 

♦  (',£)  »  E[  p\t,2>  p(t+t,*+£>  ]  (2) 


The  frequency-wavenumber  (*-&)  epectrua  of  e  eta- 
tionary,  hoeogeneoue  randoa  fiald  is  the  Fourier 
trenafore  of  tha  fiald'*  acf  with  raspect  to  temporal 
lag  and  spatial  saparation: 

P(-,fc)  ■=  /  /  ♦(!,!)  dtdC  O) 
whar#  "t"  denote*  th*  transpose  of  tha  vector. 

Although  tha  acf  and  w-fc  spectrum  specify  com¬ 
pletely  the  eecond-order  statistic*  of  a  stationary, 
homogenaous  randoa  field,  on*  or  th*  other  is  often 
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Fig.  1.  Theoretical  wavenuaber  apactra  for  two  models 
(Corcoa  and  Chaaa)  of  tha  wall  praaaura  banaath  a  TBL 
ara  shown.  Tha  abaciaaa  la  plottad  In  unlta  of  k  nor¬ 
malized  by  tha  ratio  of  radian  fraquancy  to  convactlva 
velocity.  The  ordinate  la  plottad  In  dB  ralativa  to 
tha  maximum  value  of  each  individual  epectrua. 


aora  appropriate  whan  It  la  neceaaary  to  characterize 
particular  propartlaa  of  tha  field.  For  example,  tha 
acf  la  useful  In  determining  tha  taaporal  lacs  and 
apatlal  separations  over  which  events  In  tha  field 
tend  to  be  highly  correlated.  Tha  a>-k  spectrua 
character  lias  directly  tha  average  distribution  of 
power  In  tha  field  across  fraquancy  and  wavenuaber. 
Tha  latter  characterisation  la  important  In  tha  physi¬ 
cal  description  of  tha  pressure  field  Induced  by  a 
Turbulent  Boundary  Layer  (TBL)  along  tha  wall  boun¬ 
dary. 

MODELS  OF  THE  HALL  PRESSURE  SPECTRUM 

For  aany  years,  fluid  dynaalciats  have  been 
Interested  In  tha  fraquancy-wavanuabar  characteristics 
of  fluctuations  In  wall  pressure  that  are  induced  by 
tha  Turbulent  Boundary  Layer.  This  Interest  has  lad 
to  several  nodale  of  tha  wall  pressure  epectrua,  e.g., 
Corcos  (1),  Chase  (2),  end  Ffocz-Hllllaaa  (J[). 

There  tends  to  be  general  agreeaent  aaong  theore- 
tlctane  ae  to  the  properties  of  the  wall  pressure 
epectrua  *i  *uy»tlH»btC*  cmMPPMtni  l£  convection 
tMCtolini  pressures.  At  a  given  frequency,  the 

wavenuaber  spectrua  In  the  direction  of  flow  Is  a 
aaooth,  continuous  function  of  wavenuaber  which 
reaches  a  aaxlaua  at  the  eo-called  convective  peak. 
The  location  of  thle  peak  is,  to  a  first  approxima¬ 
tion,  proportional  to  the  ratio  of  frequency  w  and 
convective  velocity  u  .  The  greatest  dlsagreeaent 
aaong  theoreticians  cSncerns  properties  of  the 
wavenuaber  epectrua  below  the  convective  peek.  This 
disagreement  centers  around  two  leeueei  the  general 
roll-off  In  spectrua  level  with  k  and  the  presence  of 
spatially  narrowband  componenta  at  vary  low 

wavenuabers. 

An  exaaple  of  the  differences  In  spectral  roll¬ 
off  lj  shown  In  Figure  One  where  the  Chase  and  Corcos 
models*  of  the  wavenuaber  spectrum  ^re  coppered  for  a 
normalised  frequency  w  of  0.3  (  w  «  ud  /tL,  where  u. 
■  10  m/s,  le  the  free-etreaa  velocity  and  4^  ■  0,0012b 
a  is  the  boundary  layer  displacement  thlckneee). 


MEASUREMENT  OF  THE  HALL  PRESSURE  FIELD 

Measurement  of  the  well  pressure  field  Is,  of 
course,  necessary  ae  e  first  step  towards  eetlaatlng 
the  w-fc  spectrua.  Various  approaches  have  been  sug¬ 
gested  for  obtaining  euch  measurements.  The  approach 
taken  In  this  paper  measures  the  wall  pressure  field 
by  an  array  of  pressure  transducers  that  are  mounted 
flush  along  the  surface  of  the  wall.  Ha  consider  only 
the  one  dimensional  w-k  spectrum  estimation  problems 
therefore,  the  gtometry  of  the  sensor  array  le  res¬ 
tricted  to  the  1  ne  and  a  uniform  separation  d  between 
sensors  Is  ass>ueed.  Nevertheless,  the  techniques 
developed  for  the  one-dl mens tonal  linear  uniform  array 
can  be  extended  directly  to  the  estimation  problem  for 
the  two-dimensional  w-l  spectrum  based  on  a  rectangu¬ 
lar  uniform  array. 

To  avoid  aliasing,  the  wall  pressure  field  must 
be  sampled  at  frequencies  and  wavenumbers  greater  than 
twice  the  highest  corresponding  frequency  and 
wavanumbere  of  the  field.  Uhen  the  spectrum  is  not 
bandllmited,  pre-filtering  le  neceesary  to  restrict 
the  bandwidth  and  thereby  minimize  the  effects  of 
aliasing.  In  the  time  domain,  such  pre-f Uterlng  Is 
accomplished  by  lowpass  filtering  tha  sensor  output 
prior  to  time  sampling.  Pre-fllterlng  In  the  apatlal 
domain  le  much  more  difficult  since  one  cannot 
manipulate  the  field  prior  to  sampling  at  the  array. 
Therefore,  the  recovery  of  the  wavenumber  spectrum  at 
a  given  frequency  Is  limited  directly  by  £  ,  tha 
separation  between  the  centers  of  adjacent  eeneore; 
components^  at  wavenumbers  above  the  Nyqulst 
wavenumber  k„  will  be  allaeed  Into  the  sampled  spec¬ 
trum.  " 


The  models  of  the  wall  pressure  w-k  spectrum  sug¬ 
gest  that  for  each  frequency  the  wavenumber  epectrum 
has  a  bandpass  characteristic  and,  except  for  the  Cor¬ 
cos  model.  Is  predominantly  one-sided.  The  sampling 
wavenuaber  for  one-sided  bandpass  processes  can  be 
reduced  by  a  factor  of  two  elnce  components  f  on  k„  to 
2k„  are  aliased  into  the  zero  region  of  the  original 
spectrum.  If  we  assume  a  one-sided  epectrua  and 
furthermore  assusw  that  the  spectrua  is  highly 
attenuated  above  the  convective  peak,  then  ths  aaxlaua 
frequency  for  which  the  w-k  spectrua  can  be  recovered 
after  spetlel  tamp ling  lei 


where  the  convective  peak  la  located  at  k  ■  w/u  . 

c 

The  effect  of  finite  eeneor  size  Is  to  lowpass 
filter  the  random  field  In  the  wavenumber  domain.  If 
the  face  of  the  sensor  Is  assumed  to  be  circular  with 


1.  For  the  Chase  model,  the  proportionality  constants 
were:  convective  velocity/f ree-streaa  velocity  ■  0.8, 
friction  veloclty/eonvectlve  velocity  ■  0.056,  and 
boundary  layer  displacement /boundary  layer  displace¬ 
ment  thickness  •  8.0.  The  model  parameters  are  from 
(2)  and  were  obtained  originally  by  fitting  the  aodel 
to  Bull's  data  (4):  v  >  0.176,  a  ■  0.766,  b  -  0.756, 
b.  *  0.378,  and  r  ■  0.389.  The  paraaetef  8  of  the 
Corcos  model  was  set  to  0.11. 


2.  The  Nyqulst  wavenumber  k  Is  defined  si  k.  •  1/2  k 
where  k  »  2t/d  Is  the  sampling  wavenumber  (also  knowX 
as  either  the  grating  lobe  wavenumber  or  the  aliasing 
lobe  wavenuaber). 
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Fit-  2.  Chasa  spectre  tn  shown  by  ths  solid  line*  For 
various  frequencies  (100  to  1000  Hz  In  100  Hz  steps) 
and  tha  spatial  responas  of  an  Ideal  transducer 
(radius  «  0.005a)  Is  shown  by  the  dashed  line.  The 
solid  vertical  Una  denotes  the  highest  wavenumber 
sampling  frequency  that  can  be  obtained  In  theory  (  d 
-  2  x  radius  );  in  practice,  the  sampling  wavenuabar 
will  be  significantly  lower  due  to  physical  con¬ 
straints  in  sensor  placement. 


radius  r  and  uniformly  sensitive,  than  the  wavenumber 
response  characteristic  of  the  sensor  can  be  written: 

7  ♦jftkt) 

|h(£)i  -  — 1 -  (5) 

Ur)2 

-  2  2(1/2) 

where  k  *  (  k.  ♦  k*  )  and  J.(.)  Is  a  first  order 

Bessel  function. 

The  width  of  the  main  lobe  of  this  lowpass  filter 
is  inversely  proportional  to  the  diameter  of  the  sen¬ 
sor.  Therefore,  greater  lowpass  filtering  of  the 
wavenumber  spectrum  is  achieved  only  by  Increasing  the 
size  of  the  sensor.  Such  an  increase  In  size,  how¬ 
ever,  Increases  the  lower  limit  of  the  sensor  spacing 
and  therefore  decreases  f  Since  the  lower  bound 

on  d  Is  the  diameter  of  tReX»ensor ,  finite  sensor  size 
cannot  act  as  a  "pre-filter"  to  attenuate  the  convec¬ 
tive  region  of  the  spectrum.  This  is  demonstrated  in 
Figure  2  where  the  wavtnjmber  response  of  the  ideal 
circular  sensor  of  radius  0.005  m  Is  shown  along  with 
a  family  of  wsvenumbsr  spectra  generated  by  the  Chase 
model  with  frequency  aa  tha  parameter.  Wavenumber  le 
normalize^  to  the  boundary  layer  displacement  thick¬ 
ness,  4  *  0.00125  m.  The  maximum  unaliaeed 

wavenumber  of  inversely  proportions!  *,  the  eps^ing  o( 
the  eeneore,  which  spacing  car  be  as  small  as  the  sen¬ 
sor  diameter  0.01  m;  this  maximum  wavenumber  le  indi 
catad  in  non-dlmtnelonal  coordinates  by  ths  solid 
vertical  line  In  Figure  2.  Finite  sensor  elze  help*, 
to  some  extent,  by  further  attenuating  components  in 
the  wavenumber  apectrum  that  are  well  abcvs  the  con 
vectlve  region  of  the  apectrum  but  it  cannot  prevant 
severe  allaalng  if  spectra  are  eetimeted  at  frequen¬ 
cies  above  f  »  u  /d  . 

max  c 

ESTIMATION  OF  THE  FREQUENCT -WAVENUMBER  SPECTRUM 

The  preceding  section  on  measurement  formulated 
tha  abaoluta  limits  for  the  recovery  of  a  continuous 

random  field  from  lta  samples  In  time  and  apses 
These  ltmite  eleo  restrict  any  aatiasta  of  tha  w-k 
spectrum  to  the  unaliaeed  region  of  the  spectrua. 
Besides  these  sampling  limits,  tha  fact  that  ths  field 
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is  measured  over  a  finite  period  of  time  by  e  finite 
aperture  array  Implies  that  even  within  this  unaliaeed 
region  of  the  spectrum,  errorless  recovery  le  Impossi¬ 
ble;  statistical  errors  In  tha  estimator  and  res¬ 
tricted  samples  of  the  acf  limit  tha  performance  of 
any  procedure. 

Tha  spectral  estimation  problem  Is  posed  with 
respect  to  limitations  In  tha  measurements:  given  a 
finite  sample  at  discrete  points  In  time  and  apace  of 
a  random  field,  m  are  Interested  In  processing  the 
data  to  provide  tha  st*  estimate  of  tha  w-k  spec¬ 
trum. 

Statistical  eb.lmatora  g(jt)  are  functions  that 
operate  on  data  x  to  eatimate  soma  characteristic  or 
parameter  (8)  of  the  data.  Since  tha  data  are  random, 
each  g(x>  "will  have  a  probability  distribution  func¬ 
tion  (pdf)  that  depends  directly  on  tha  stochastic 
representation  of  the  data.  Statistical  estimators 
have  many  properties  by  which  they  era  Judged  to  be 
good.  The  bias  b  (8)  and  variance  var  (g)  of  an  esti¬ 
mator  (  8=  g(x))  Ire  the  flret  and  iecond  central 
momenta  of  pdf,  i.a., 

bg(g)  -  E[  g(g)  -  8  ]  (6) 

var|(8)  -  Et  (g(g)  -  g)2  1  O) 

There  are  numerous  ways  to  estimate  a  given 
characteristic  of  tha  data.  In  general,  bias  and 
variance  are  criteria  that  era  used  to  determine  which 
estimator,  from  a  given  sat,  la  best.  In  the  Ideal 
caee,  the  beet  estimator  1s  unbiased  and  exhibits  the 
least  variance  among  those  estimators  In  the  eat;  in 
practice,  this  is  not  alwaye  true.  Optimum  estimation 
techniques,  such  as  maximum  likelihood  estimation,  are 
often  useful  in  deriving  a  "beat*  estimator  for  a 
given  problem.  For  this  reason,  the  development  of 
spectrum  estimation  theory  has  been  baaed  on  "empiri¬ 
cal"  research  in  the  sense  that  the  theory  consists  of 
numerous  forms  of  estimators  each  of  which  performs 
well  for  a  given  class  of  problems  but  none  of  which 
is  necessarily  optimum  among  the  clasa  of  ell  possible 
estimator,,. 

This  section  reviews  three  general  forma  of  spec¬ 
tral  estimator  that  have  been  proposed  In  tha  litera¬ 
ture  and  diecuetet  tha  performance  of  each  with 
respect  to  blae  and  variance  for  the  Corcos  end  Chase 
models  of  the  well  pressure  w-k  spectrum.  We  conaider 
ths  problem  of  blae  and  variance  separately.  Ini¬ 
tially,  the  effects  of  truncating  tha  acf  on  tha 
deterministic  component  of  the  estimator  blae  are  stu¬ 
died  by  computing  the  estimated  apectrum  based  on 
exact  values  of  the  truncated  ecf.  Later,  we  consider 
bias  and  variance  in  greater  detail  by  uae  of  Bonte 
Carlo  simulation. 

Without  lose  of  generality,  each  estimator  la 
introduced  by  considering  the  estimation  of  the  power 
spectrum  of  a  stochastic  process:  the  results  ere  the 
same  for  the  wavenumber  spectrum  of  a  stochastic  spa¬ 
tial  field. 

Let  r(r).  be  the  autocorrelation  function  of  the 
etocb, astir  p-ocese.  and  let  r(:.)  be  its  sampled  ver¬ 
sion.  In  practice,  rii^i  must  be  eetimeted  from  the 
sampled  prcceee  by  some  approprlsts  technique.  Let 
'  ,  )  denote  such  an  aatimats.  Finally,  assus*  that 

p.wj,  tha  power  apectruu,  is  a  lowpass  process  and 
that  (he  sampling  f-squenev  is  sufficiently  high  to 
prevent  aliasing. 


Fourier  Transform  al  lA*  Truncated  *£[.!  Blackaan-Tuker 
IM  Parlodograa  r-'tbodl 

Tha  Blackaan-Tukey  method  estimates  P(w)  by  per¬ 
forming  a  dtaerata  Fourlar  tranafcra  (OFT)  on  a  wln- 
dowad  version  of  length  L  of  tha  aatlaatad  acf: 

PK(a)  -  DFT[  w<tj)  ?(tt)  ]  (6) 

Ora  coaputational  disadvantage  of  tha  Blackaan- 
Tukey  aathod  la  that  tha  saapla  acf  auat  ba  aatlaatad 
flrat  lroa  tha  data.  Tha  parlodograa  la  an  altarna- 
tlva  procadura  that  aatlaataa  tha  power  apactrua 
directly  froa  tha  data  x(t  )  by  segmenting  tha  data 
Into  !'  sections  of  length  M,  performing  a  OFT  on  aach 
aactlon,  and  a  (raging  ovar  aactlon  tranaforaa.  Tha 
aatlaata  of  thi  powar  apactrua  la  tafcan  aa  tha  average 
of  tha  squared  aagnltuda  of  tha  tranaforaa.  Thuai 

Pp(w)  =  (1/S)  l  1  DFT[  B.fTj)  ]  |2  (9) 

TT a  Black san-Tukey  and  parlodograa  aathoda  ylald 
different  aat  aatad  apactra  but,  on  tha  avaraga,  tha 
period! graa  aa  hod  la  tha  aaaa  aa  tha  Blackaan-Tukey 
aathod  with  a  Bartlatt  (triangular)  window  appllad  to 
tha  flrat  N-l  lags. 

Blaa.  The  affact  of  truncating  tha  acf  at  lag  L 
on  oatlaator  blaa  can  ba  atudiad  by  aubatltutlng  tha 
axact  axpraaalon  of  r(-t  )  and  coaputlng  tha  P.-lu). 
For  an  arbitrary  window,  tha  P„(u)  can  ba  re-written 
In  tha  frequency  doaialn  aa  the  convolution  of  tha  true 
apactrua  with  '.ha  window  function! 

PBT(u)  -  H(u)  »  P(w)  (10) 


tha  expected  aatlaata  of  tha  powar  aubjact  to  tha  con- 
atralnt  that  tha  apactral  component  at  a  la  paaeed 
with  unity  galm 

aln  E[  I jjh(w)j| 2  ]  with  *h(w)*(w)  -  1  (16) 


Thla  alnlalsatlon  problaa  can  ba  aolvad  ualng  Lagrange 
aultlpllera  to  ylald: 


H<“>  “  ^ - 1] - 

«.(->«  »(a) 


(17) 


and 


W-) 


1 

lh(a)t_1t(a) 


(18) 


Capon* a  aathod  haa  alao  bean  termed  tha  Haxlaum 
Likelihood  Method  (HLH)  although  It  la  not  tha  maximum 
likelihood  aatlaata  of  powar  apactral  danalty. 
Mevertheleaa,  tha  term  HLH  contlnuea  to  ba  uaed  whan 
referring  to  Capon*a  method  and  wa,  too,  ahall  adhere 
to  thla  (alaleading)  convention. 


Blaa.  Truncating  tha  acf  In  tha  Fourier-based 
aathoda  la  equivalent,  for  tha  deteralnletlc  caaa  to 
convolving  tha  power  apactrua  with  that  of  tha  lag 
window.  Tha  effecta  of  truncation  on  tha  HLH  apactrua 
are  more  difficult  to  charactarlxa  alnca  tha  "window* 
change*  aa  a  function  of  frequency  and  depend*  on  tha 
covariance  matrix  K.  In  general,  blaa  can  only  ba 
evaluated  by  comparing  tha  HLH  apactrua:  with  example* 
froa  the  particular  claa*  of  apactra  that  are  under 
conaldaratlon. 


where  dano  .a*  tha  convolution  operator.  Thua,  the 
resulting  spetral  estimate  la  aaoothed  dua  to  (Bear¬ 
ing  of  adjacent  components  In  the  apactrua  by  tha  aaln 
lobe  and  to  leakage  of  more  distant  components  through 
tha  sldalobas  if  W(u).  Bias  can  be  reduced  to  arbi¬ 
trarily  aaall  values  but  only  by  taking  more  sample* 
of  tha  acf. 


Cmn’i  Bgaliuq  UllLLAaad  Hat  hod. 


Tha  ganai al  fora  of  tha  averaged  parlodograa  i*i 

(11) 


P  (w)  -  (1/S)  £  I  hNu)^  I2 


where  superscript  "h"  denotes  tha  conjugate  transpose, 
and,  for  the  specific  case  of  tha  parlodograa, 
g(w)  *  £(w),  tha  discrete  Fourier  operator! 


*(«)  »  (  1,  a 


Ju(H-1)t 


V 


(12) 


and  g(  la  tha  data  Xj  wrlttan  In  vector  fora: 

*.  ■  <  Vi . Vh  >*  <l3> 

Taking  t ha  expectation  of  both  aldaa  of  Eq.  (11) 
yield*! 

P(w)  •  ih(w)  R  f(»)  (14) 

P 

which  1*  a  quadratic  fora  Involving  the  covariance 
matrix  R  of  the  randoa  vector  gi 

R  -  Et  i  ih  ]  (15) 


Maximum  Entropy  Method 


Given  the  eaaple  acf,  r(t  ),  truncated  after  lag 
L,  an  Infinite  nuaber  of  P(a)  exist  that  have  r(t.)  as 
the  first  L+l  term*  In  thalr  acf  sequence  (j).  1  The 
Haxlaua  Entropy  Method  selects  that  P'w)  from  this 
Infinite  set  which  has  tha  maxlaum  entropy.  This 
optlalsatlon  problaa  can  ba  written  asi 

max  /  In  P(w)  dm  (19) 


such  that  P(u)  aattsflee 

r(tj)  «  /  P(w)  e 


JV 


•,  1  -  0,  L 


(20) 


For  the  one-dlaenelonal  case  with  uniformly 
spaced  aaaples,  the  HEH  solution  reduces  to  another 
fora  of  spectral  estimator,  tha  all-pole  or  auto¬ 
regressive  (AR)  estiMtor,  In  which  the  wavenumber 
apactrua  Is  approximated  by  an  all-pola  spectrua. 
Thus, 


PHEH<“> 


I  1 


»*!<•> 1 2 


(21) 


where  the  coefficient  vector  g  Is  determined  from  the 
covariance  matrix  R  ($)i 


B 


(22) 


with 


r<vL))t 


In  the  parlodograa  aathod,  g(w)  has  the  same  fora 
regardless  of  the  particular  w.  In  Capon’ a  formula¬ 
tion,  g(w)  is  chosen  for  each  frequency  to  alntalse 


t  •  (r(tj),  r<tj> 


(23) 


CSD  using 


gj||.  By  maximizing  the  entropy,  the  HEM  eolu- 
tlon  forcaa  tha  laaat  amount  of  "structure"  onto  tha 
correlation-matching  P(w).  In  tha  caaa  of  one- 
dlaenslonal  unlfora  sampling,  thla  aolution  iapllaa 
that  an  all-pole  nodal  of  ordar  L  la  the  aoit  antroplc 
choice.  It  la  not  clear  how  blaa  and  antropy  are 
related.  If,  In  fact,  tha  apectrua  la  all-pole  of 
ordar  L,  then  tha  MEM  aolution  la  exact  and  there  la 
no  blaa.  Uhen  tha  apectrua  la  unknown  and  few  aseump- 
tlone  can  be  aada  g  priori  aa  to  tha  appropriataneaa 
of  an  all-pole  aodal  of  (Ivan  ordar,  tha  question  of 
blaa  raaalna  unanswered.  Considerable  attantlon  has 
bean  paid  to  thla  lasua  and  It  raaalna  an  actlva  area 
of  research. 

gp.allnp.tr  Pata-adaptlva  procedures. 

Both  HLH  and  AR  ars  examples  of  data-sdaptiva 
procedures  In  which  tha  weighting  vector  (w)  or  tha  AR 
coefficient  vector  (4)  depend  nonlinearly  on  tha 
aaapla  covariance  R.  Thla  contrast*  with  tha  llnaar 
dapandsr.c*  of  tha  two  Fourlar  aathods  on  R. 

Intsrast  In  nonlinaar  dats-adapt ive  procedural 
davslopad  quickly  whan  It  was  shown  in  tha  ald-1960'a 
that  such  aathods  ar*  batter  able  to  resolve  cloaaly 
spaced  peaks  In  tha  apectrua  than  Fourlar  aethoda  for 
equal  nuabere  of  lags  and  high  elgnal-to-nolae  ratio. 
Tha  tara  "high-reeolutlon"  has  often  been  applied  to 
eetlaatore  such  aa  AR  and  HLH  with  rafaranca  to  thalr 
superior  resolution  of  spectra  that  conslat  of  several 
line  coaponants  In  a  whits  background  nolsa.  Tha 
study  of  spectral  estimation  has  lines  bean  doalnstad 
by  tha  problaa  of  resolving  peaks  In  tha  apectrua 
under  conditions  of  vary  few  noisy  saaplaa  of  tha  acf. 
Therefore,  aany  of  tha  high-raeolutlon  spectral  estl- 
matora  that  have  bean  developed  over  tha  peat  decade 
should  bo  aors  appropriately  labeled  paranetar  estima¬ 
tors  since  tha  aaasura  of  psrforaanca  reflects  the 
accuracy  with  which  tha  eatlMtor  locataa  tha  spectral 
peaks. 

In  general,  the  apectrua  predicted  by  any  of  tha 
aodels  of  tha  wall-praaaura  field  is  not  doalnatsd  by 
narrowband  peaks  nor  does  the  Interest  11a  particu¬ 
larly  In  aaaeurlng  tha  apectrua  within  the  region  of 
the  broad  convact’ve  peak.  Inataad,  tha  raaaarch 
focus  la  on  eetlastlng  accurately  the  wavenuabar  spec- 
trua  over  a  broad  range  of  k.  Thus,  tha  estimator 
auat  exhibit  alnlaua  bias  or  high  fidelity  over  t  tar- 
gettad  region  of  tha  apectrua  rather  than  raaolva 
peaks  that  are  located  throughout  tha  apectrua.  Since 
both  HLH  and  AR  aodal ing  are  nonlinaar  data-adaptlve 
procedures  that  have  baan  designed  for  resolution 
(parameter  estimation)  rather  than  spectral  fidelity 
(spectral  estimation),  it  la  not  claar  whether  resolu¬ 
tion  for  line  spectra  extends  to  fidelity  for  broad, 
saooth  spectra. 

Spectral  asJlaatlop;  taaporal  component  a£  the  w-k 
■mactrua. 

Tha  w-k  apectrua  can  ba  sstlaatad  in  one  step  by 
applying  a  multidimensional  version  of  one  of  tha 
three  aathods  discussed  above,  or,  eltsrnat lvely ,  In 
two  stapa  by  first  estimating  the  cross-spectral  den¬ 
sity  (CSD)  as  a  function  of  frequency  whara 

r(w,E)  -  /  ♦(t.t)  e'J“t  dt  (24) 

and  than  estimating  tha  wavenuabar  apectrua  froa  tha 


P(w,k)  -  f  r(ta»,  C)  dC  (25) 

Because  of  Ite  computational  afficiancy  and 
wall-developed  theory  of  blaa  and  variance,  tha 
Fourlar  approach  la  generally  preferred  whan  there  1* 
aufficlent  data  available  for  processing.  In  tha  case 
of  TBL  studies,  the  amount  of  data  sampled  over  time 
la  Halted  primarily  by  tha  availability  of  storage 
specs  and  by  tha  tlaa  span  of  tha  experiment.  Tha 
nuabar  of  spatial  saaplaa,  however,  la  usually  vary 
restricted  both  because  of  tha  difficulties  associated 
with  calibrating  and  maintaining  large  arrays  and 
bacauaa  of  tha  relatively  email  range  over  which  tha 
pretsure  field  can  ba  assumed  to  ba  homogeneous. 
Therefore,  a  two-staga  approach  to  tha  estimation 
problem  la  appropriate  In  which  tha  frequency  com¬ 
ponent  Is  estimated  froa  DFTa  of  the  aaapled  tlaa 
series  of  each  of  tha  aanaora. 

Spectral  estimation:  apetiai  component  &i  ItU  ur-k 
spectrum. 

Assume  that  tha  pressure  (laid  is  aansad  by  an 
M-elaaant  llnaar  unlfora  array.  Let  s#(Ja»)  be  tha  e- 
th  anapshot  of  tha  llnaar  array  at  frequency  w,  i.e., 
2  (Ju)  la  a  complex  M-d leans lonal  vector  with  com¬ 
ponents  x  (Ju)  equal  to  1-th  tensor  output  that  la 
aval j*ted*'*t  complex  frequency  Ju  over  the  a-th  tlaa 
segment  of  the  data.  Aa  noted  above,  x  (Ju)  can  ba 
obtained  by  performing  a  DFT  on  tha  t lae-aagaented 
data.  Tha  sasiple  cross-apectral  density  matrix  (CSDH) 
at  frequency  u,  R(u),  la  foraad  by  summing  acroaa  tha 
outer  product  of  tha  snapshots: 

R(u>  *  (1/S)  l  ^(Ju)  s^(Ju)  (26) 

For  notatlonal  convenience,  the  dependence  of  R  and  x 
on  u  la  assumed  Implicitly,  ao  tha  argument  u  la 
dropped. 

A  perlodogram  aat lasts  of  tha  wavenuabar  apectrua 
Is  equivalent  to  appropriately  steering  tha  beam  of 
the  arrey  for  a  desired  wavenuabar  k  and  taking  tha 
squared  magnitude  of  tba  atasred  array  output  as  tha 
estlaata  of  tha  wavenuabar  apectrua: 

P  (k)  »  £h(k)  R  i(k)  (27) 

P 

whara  £00  la  the  discrete  Fourier  operator  with 
respect  to  the  spatially  aaapled  data. 

Recall  that  tha  primary  concern  with  tba  Fourier 
approach  la  tha  potential  bias  that  can  exist  due  to 
smearing  of  tha  true  spectrum  by  tha  main  lobe  of  the 
lag  window  and  leakage  of  tha  apectrua  through  the 
side  lobes.  For  a  aaall  nuabar  of  aanaora,  tha 
degrees  of  fraadoa  available  for  designing  a  suitable 
lag  window  are  restricted  ao  that  an  array  response 
with  sufficient  selectivity  at  tha  k  of  tntaraat  and 
nulla  away  froa  that  k  la  difficult,  if  not  impossi¬ 
ble  ,  to  design  for  arbitrary  spectra  and  for  all  k. 

Spatial  analogs  for  AR  and  for  HLH  follow 
directly  froa  tha  taaporal  results  by  denoting  tha 
output  of  the  spatial  aanaora  at  a  given  u  as  tha  data 
and  substituting  the  aaapled  CSD  and  CSDH  for  £(tj) 
and  tha  covariance  matrix  R,  respect lvely . 

Tha  performance  of  tha  perlodogrea  la  compared 
with  HLH  and  AR  aodal ing  in  tha  two  panala  of  Figure  3 
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Fi*.  3.  Th*  performance  of  three  spectral  estimators, 
Fouriar  tranaf oraation  (FT),  Haxlaua  likelihood  Hat hod 
(HLN),  and  Auto-Regreaal ve  Nodal In,  (AX)  ara  coaparad 
for  tha  Corcoa  and  Chaaa  aodala  of  the  wall-praaaura 
wavanuabar  apactrua.  Tha  coaputad  apactra  ara  baaad  on 
axact  valuaa  of  tha  CSDN  for  a  12-aanaor  llnaar  unt- 
fora  array. 

for  tha  Chaaa  and  Corcoa  aodala  (dot tad  llnaa)  of  tha 
wall  praaaura  wavanuabar  apactrua.  For  coaparlaon 
purpoaaa,  tha  aathoda  ara  evaluated  baaad  on  a  llnaar 
unifora  array  of  12  aanaora  .  Tha  raaulta  ara  baaad 
on  tha  axact  CSDN  for  aach  aodal.  Each  raault  la 
rafarancad  with  raapact  to  ita  paak. 

Tha  parforaanca  of  tha  parlodograa  la  raaaonably 
accurata  for  tha  Corcoa  aodal  but  ia  conaidarably 
biaaad  for  tha  Chaaa  aodal.  Tharafora,  tha  Corcoa 
aodal  la  aufflctantly  broad  and  aaooth  that  convolving 
it  with  tha  wavanuabar  apactrua  of  tha  lag  window, 
a.g.,  tha  baaapattarn  of  tha  12-aanaor  array,  haa  lit— 
tla  affact  on  tha  final  aatiaata.  On  tha  othar  hand, 
tha  parlodograa  ia  unabla  to  follow  tha  apactral 
roll-off  of  tha  Chaaa  aodal  dua  prlaarily  to  tha 
affacta  of  aidaloba  leakage.  Tha  diffarancaa  batwaan 
tha  two  coaputad  apactra  ara  vary  ali«ht.  If  tha 
affacta  of  atatlatlcal  arror  in  eetiaatlng  tha  CSDN 
ara  includad,  tha  dlatrlbutiona  of  tha  aatlaator  for 
tha  two  apactra  ovarlap  conaidarably.  Tharafora,  tha 


3.  Aa  tha  final  goal  of  tha  eatiaetion  daaign  ia  to 
procaaa  data  froa  tha  SAFHL  (Saint  Anthony  Falla  Hy- 
draulica  Laboratory)  taat  facility  at  tha  Univeraity 
of  Ninnaaota,  all  aathoda  wara  avaluatad  within  tha 
llatta  of  tha  nuabar  of  aanaora  aval labia  for  tha  ax- 
par  iaant. 
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Fig.  4.  Tha  parforaanca  of  two  apactral  aatiaatora, 
HLN  and  AR,  ara  coaparad  whan  tha  CSDN  la  aatiaatad 
froa  25  anapahota  of  tha  12-aanaor  llnaar  unifora 
array.  Tha  Chaaa  apactrua  la  ahown  by  tha  daahad 
llnaa  and  tha  four  alaulatlona  for  each  aatlaator  ara 
ahown  by  tha  aolld  llnaa. 


parlodograa  (beaaforalng)  aathod  la  Inadequate  aa  an 
aatiaata  of  tha  wavanuabar  apactrua,  particularly  in 
tha  vary  low  wavanuabar  raglaa. 

Tha  AX  aodal  and  tha  HLN  follow  tha  general  trend 
of  tha  two  aodala  qulta  wall.  Khan  rafarancad  to  the 
paak  level ,  HLN  providea  an  uppar  bound  on  tha  under¬ 
lying  apactrua  in  both  caaaa.  AX,  on  tha  other  hand, 
axhibite  laaa  biaa,  on  tha  aver ago,  particularly  for 
tha  Chaaa  aodal,  but  alao  introducaa  rippling  In  both 
of  the  coaputad  apactra.  Thle  rippling  raflacta  tha 
poor  fit  batwaan  thaaa  apactra  and  onaa  generated  by 
all-pole  aodala  of  order  11.  Tha  two  aathoda  differ 
conaidarably  with  raapact  to  blaa  In  tha  vary  low 
wavanuabar  region  of  tha  apactrua.  HLN  a how*  .ora 
than  10  d>  of  blaa  for  wavanuabera  that  are  laaa  than 
1  decada  below  tha  convective  paak.  AX,  on  tha  othar 
hand,  exhiblta  little  biaa  over  thia  region. 

gitlmtt  Si  1M  aatiaataa.  It  la  difficult,  in 
general,  to  obtain  a  cloaed-fora  expreaalon  for  the 
variance  of  tha  AR  apactral  aatiaata  for  the  aaall- 
aaapla  caaa.  Therefore,  Nonte  Carlo  alaulatlona  wara 
par f oread  to  etudy  tha  eenelttvity  of  AR  and  of  HLH  to 
at  itletlcal  errora  In  the  aatiaata  of  tha  CSDN,  X, 
f r  '■  th«  data.  Xneteed  of  eiaulatlng  tha  data  vector 
x  jndar  <»  directly,  R  waa  alaulatad  uatng  a  coaplax 
V! ahart  diatributlon  (7).  Figure  4  coaparea  tha 
raaulta  of  four  aiaulatione  with  tha  actual  apactrua 
whan  R  la  aatiaatad  froa  25  anapahota  of  tha  data.  Aa 
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Fit-  S.  95X  confidence  intervale  baeed  on  200  simula- 
tiona  of  the  AR  and  HLH  estimators  are  shown  by  the 
solid  lines  and  the  siean  is  shown  by  the  dashed  line. 
Monte  Carlo  simulations  ware  based  on  25  snapshots 
from  a  12-sensor  linear  uniform  array- 


can  be  seen,  the  effects  of  statistical  error  on  the 
performance  of  AR  are  substantially  greater  than  those 
on  HLM.  Figure  5  shows  95%  confidence  intervals  that 
are  based  on  200  simulations  which  support  these  con¬ 
clusions. 

Evaluation  the  three  methods 

The  effect  of  truncating  the  cross-spectrum  on 
the  bias  can  be  reduced,  for  each  estimator,  only  by 
increasing  the  total  number  of  sensors.  Since  this 
value  is  determined  by  the  constraints  of  the  experi¬ 
ment,  a  perlodogram  (beamforming)  approach  cannot 
recover  the  underlying  spectrum  with  sufficient  fidel¬ 
ity  in  the  low-wavenumber  region.  The  truncation 
affects  are  less  severe  for  both  HLH  and  ARi  HLH 
overestimates  the  underlying  spectrum  whereas  AR  exhi¬ 
bits  fidelity  which  is  superior  to  that  of  HLH,  par¬ 
ticularly  at  low  wavenumbers  for  the  Chase  spectrum, 
but  introduces  false  "peaks’  or  "ripples"  throughout 
the  computed  spectrum.  Nonetheless,  AR  is  consider¬ 
ably  more  sensitive  to  statistical  srrors  than  HLH. 
Thus,  the  superior  fidelity  of  AR  can  be  realized,  in 
practice,  only  by  increasing  substantially  the  number 
of  snapshots  used  to  estimate  the  CSDH.  HLH,  on  the 
other  hand,  is  less  sensitive  to  these  effects  but  is 
also  more  biased  than  AR  in  the  region  of  interest. 

In  addition  to  statistical  srror,  the  performance 
of  each  method  can  be  degraded  by  mismatches  in  phase 
and  amplitude  across  the  sensors  and  by  errors  in  sen¬ 
sor  placement.  Studies  of  these  effects  have  shown 
that  AR  is  also  more  sensitive  than  HLH  to  these 
errors  in  the  design  of  the  linear  array  (8,  9). 


SUBSPACE  EXTENSION  OF  CAPON'S  METHOD;  se-HLH 

The  relatively  stable  behavior  of  HLH  suggests 
that  it  is  a  useful  estimate  of  the  upper  bound  of  the 
true  we  pressure  wavenumber  spectrum  at  a  given  u 
In  the  ollowlng,  we  develop  a  new  method  whereby  the 
performance  of  HLH  can  be  improved  over  a  targeted 
region  of  the  spectrum. 

Consider  the  eigen-decomposition  of  the  CSDH  R: 

*  •  I  »,  Ij  l}  (28) 

where  k  and  |,  are  the  i-th  eigenvalue  and  eigenvec¬ 
tor  oflR,  respectively.  Substituting  this  expression 


Fig.  6.  Eigsnbeams  for  the  Chase  model  (dashed  line) 
based  on  a  12-sensor  linear  uniform  array  are  shown  by 
the  solid  lines  for  the  first  three  largest  eigen¬ 
values  . 


into  Eq.  (18), 
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find 


!h(k)<  l  1/kj  ij  >*(k> 


(29) 


PHLH(1°  = 


Z  1/A, 


Zj  (k) 


(30) 


where  Zj(k)  =  z  f(k).  Thus,  the  HLH  estimate  at  aach 
k  is  inversely  proportional  to  the  weighted  sum  of  the 
power  in  each  Fourier  transformed  eigenvector  at  that 
k. 


Let  us  consider  the  spectral  properties  of  the 
eigenvectors  of  the  Chesj  CSDH.  In  Figure  6,  the 
eigenbeams,  e.g.,  |  Z.(k)  I  ,  of  the  eigenvectors  that 
correspond  to  the  three  largest  eigenvalues  of  the 
Chase  spectrum  are  shown.  Thsss  eigenbeams  vary  with 
eigenvalue  dominance.  The  main  lobe  of  the  largest 
eigenvector  (the  eigenvector  with  the  largest  eigen¬ 
value)  is  centered  on  the  convective  peak  of  the  spec¬ 
trum.  For  the  second  and  third  eigenvectors,  the 
malnlobe  splits  in  two  and  a  slight  null  is  introduced 
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•t  the  convective  peak.  For  smaller  eigenvectors  (not 
ahown  In  tha  Figure)  a  single  main  lobe  structure  re- 
appaara  ft rat  at  low  wavenumbers  and  finally  in  tha 
negative  wavenuaber  region  of  tha  apactrua.  Along 
with  thia  algratlon  in  tha  location  of  tha  aainlobe, 
tha  null  in  tha  elgenbaaa  bacoaaa  atrongar  and  broadar 
at  tha  lowar  aiganvaluaa  and  la  local lsad  in  tha 
nalghborhood  of  tha  convactiva  paak. 

Reduction  ai  LLit;  SMtom 

Ao  arguad  abova,  tha  parfornanca  of  tha  HLH  ovar 
a  apacific  raglon  of  tha  wavanuabar  apactrua  la  detar- 
ainad  by  tha  ralativa  doalnanca  of  tha  waightad  powara 
of  tha  aiganvactora  in  that  raglon.  Tbua,  an  laprove- 
want  in  the  fidelity  aay  be  achieved  by  perforaing  tha 
HLH  on  an  approprlata  aubapaca  of  R,  rather  than  on 
tha  full  apace. 

To  dataralna  an  ouch  appropriate  aubapaca,  tha 
parforaanca  of  tha  HLH  ovar  all  poaalble  aubapacea  of 
R  auat  be  dateralnad  for  the  target  tad  region  of  tha 
apactrua.  Given  a  aultabla  choice  of  error  aatric 
d( . , . ) ,  the  problaa  becoaaa: 

aln  d(  p(k),  p(k)>  (31) 


Fig.  ?.  Tha  parforaanca  of  tha  HLH  ualng  the  subapar' 
of  t' a  CSDH  apanned  by  tha  10  aaallaat  aiganvactora, 
i.e.,  tha  aiganvactora  correapondlng  to  th<  10  eaal- 
laat  aiganvaluaa,  la  coaparad  with  tha  Cha  e  apactrua. 
Each  solid  line  denotaa  tha  slaulated  parforaanca  of  a 
12-aenaor  linear  unlfora  array  ualng  25  snapshots  of 
the  data. 


whore  the  alnialzatlon  la  perforaed  ovar  all  possible 
aubapacea  R  of  R  and  p(k)  la  taraad  tha  design  spec- 
trua.  Tha  sSbspace  astlaate  takas  the  forai 


Pss-MLH(I° 


l  1/A j  |  *j(k)  |‘ 


where  the  euaaatlon  extends  over  tha  indices  of  tha 
choaan  aubapaca  R(. 

Tha  search  procedure  for  tha  optlaal  aubapaca  is 
ovar  all  posaible  coabinetlona  of  tha  CSDH  oigonspacaa 
and  thus  can  bacoae  coaputationally  intanslva  ovan  for 
a  aoderate  nuaber  of  sensors.  An  alternative  search 
procedure  is  used  for  designing  tha  aubapacea  for  tha 
Chase  and  Corcos  spectra.  This  procadura  searches 
through  an  eabedded  set  of  aubapacea  by  reaoving,  at 
each  stags,  the  elgenapace  of  the  largest  eigenvalue 
of  the  subspace  and  selecting  as  tha  "optlaal"  sub¬ 
space  that  which  alnlalzes  the  naan-square  error 
between  the  estimate  and  the  design  apactrua  over 
wevenuabers  between  0  and  the  convective  peak. 

Figure  7  shows  results  froa  four  Honte  Carlo 
elaulatlona  in  which  the  HLH  waa  perforaed  on  tha 
subspace  spanned  by  the  10  aaalleat  eigenvectors  of 
tha  ostlaatod  CSDH.  The  procedure  reduces  the  bias 
froa  10  tc  4  dB  within  one  decade  below  the  convective 
paak.  Furthernore,  tha  procedure  la  r  t  overly  sensi¬ 
tive  to  statistical  error  as  waa  the  case  for  AR. 
Slailar  results  were  obtained  for  the  Corcos  model. 
For  this  apactrua,  however,  the  alnialzlng  aubapaca 
was  spanned  by  tha  11  aaalleat  eigenvectors. 


Since  tha  weighting  acbeaa  depends  on  tha  exact  eigen¬ 
values,  statistical  errors  In  the  aatlaated  eigen¬ 
values  are  equivalent  to  "de-tuning"  the  reflnaaent 
procedure.  Therefore,  tha  success  of  fractional 
weighting  depends  to  a  large  extant  on  the  sensitivity 
of  tha  spectral  error  aatric  to  perturbations  in  the 
weights  which,  in  turn,  is  a  function  of  the  design 
spectrua. 

CONCLUSIONS 

Estiaatlon  of  the  frequency-wavenuabar  apactrua 
of  a  randoa  field  is  a  difficult  problea  when  rela¬ 
tively  little  data  is  available.  Sources  of  error  in 
tha  estlaatad  spectrua  include  truncation  of  the  scf 
and  statistical  errors  in  estlaatas  of  the  acf.  Non¬ 
linear  spectral  estiaators  exist  and  are  known  to  par- 
fora  wall  for  certain  classes  of  fraquancy-wavenuabar 
spectra  under  a  resolution  criterion.  These  estima¬ 
tors,  however,  are  not  necessarily  tha  bast  choice  for 
all  classes  of  spectra.  In  particular,  the  class  of 
saooth,  broad  spectra  that  are  thought  to  be  charac¬ 
teristic  of  the  wall-pressure  fluctuations  beneath  a 
Turbulent  Boundary  Layer  are  estlaatad  with  greater 
fidelity  by  the  Haxtaua  Likelihood  Hethod  rather  than 
by  Auto-Regressive  modeling,  i.e.,  tha  Haxlaua  Entropy 
Hethod  for  the  case  of  a  linear  unlfora  array. 
Further  improvement  in  the  fidelity  of  tha  estimate 
can  be  realized  by  performing  the  Haxlaua  Likelihood 
Hethod  on  an  appropriately  selected  subspace  of  the 
cross-spectral  density  matrix  R. 


Further  Kfflntunt 

The  subspace  technique  can  be  further  refined  by 
Including  fractional  weighting  of  the  sigenspaces. 
Under  such  an  approach,  each  elgenapace  la  assigned  a 
value  be* ween  0  and  1  rather  than  tha  binary  value  0 
or  1.  Thf  advantage  of  fractional  weighting  Is  that  a 
better  fit  to  the  design  spectrua  over  the  targetted 
region  can  be  achieved.  This  refineaent  of  the  pro¬ 
cedure  la  far  more  sensitive  to  statistical  errors. 
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ABSTRACT 

Spectral  estimation  for  a  class  of 
smooth,  broad  spectra  generated  by  the 
Turbulent  Boundary  Layer  is  discussed  in 
terms  of  spectral  fidelity  and  sensitivity 
to  statistical  error.  When  the  number  of 
sensors  is  limited,  which  is  often  the 
case  in  array  problems,  MLM-based  estima¬ 
tors  exhibit  better  performance  than  all¬ 
pole  models.  These  performance  differ¬ 
ences  underscore  the  fact  that  for  certain 
classes  of  spectra  high-resolution  and 
h i gh-f I de 1  I ty  are  competing  criteria. 

1 NTRODUCT ION 

The  past  twenty  years  have  witnessed 
a  growth  in  the  number  of  spectral  estima¬ 
tors  used  in  array  processing,  particu¬ 
larly  in  the  area  of  source-angle  estima¬ 
tion.  Although  this  problem  and  Its 
time-series  equiva'ent,  the  estimation  of 
frequencies  of  sinusoids,  have  often  been 
approached  withlr,  the  context  of  spectral 
estimation,  both  are  more  appropriately 
characterized  as  paramet er -est i mat  I  on 
problems  for  which  source-wavenumber  (or 
frequency)  serves  as  the  unknown  parame¬ 
ter.  New  estimators  of  source  angle  are 
often  eva'uated  on  the  extent  to  which  the 
spectrum  is  resolved  into  peaks  that 
correspond  to  the  correct  bearing  of  each 
cf  the  sources.  In  the  more  general  case 
of  spectral,  as  opposed  to  parameter, 
estimation,  the  'typical'  spectrum  may  not 
be  dominated  by  strong  peaks  so  that  esti¬ 
mators  which  exhibit  high-resolution  may 
exhibit  worse  performance  than  others 
because  of  poorer  fidelity  in  regions  away 
from  the  spectral  peaks.  Our  work  *ocuses 
on  the  problem  of  spectral  estimation  as 
opposed  to  parameter  estimation,  and  is 
motivated  specifically  by  a  measurement 
problem  In  fluid  dynamics. 

WAVENUMBER  SPECTRUM  OF  THE  TBL 

For  well  over  a  century,  turbulence 
has  been  a  dom  Inant  rocu«,  0f  research  in 
fluid  dynamics,  both  f'or  a  theoretical  as 
well  as  measurement  standpoint.  Due  to 
its  stochastic  nature,  most  efforts  have 
been  directed  to  characterizing  the  first 
and  s< roi  u  order  statistics  of  the  tur¬ 


bulent  field.  Models  based  largely  on 
theory  predict  a  broad,  unlmodal 
wavenumber  spectrum  In  the  direction  of 
streamwlse  flow  for  wall  pressure  beneath 
a  turbulent  boundary  layer  (TBL).  For 
wavenumbers  near  the  convective  peak,  the 
models  are  consistent  with  physical 
measurement,  but  for  low  wavenumbers, 
where  reliable  data  is  almost  nonexistent, 
these  models  differ  significantly  with 
respect  to  level  and  smoothness  of  the 
spectrum.  Although  the  spectrum  level  at 
low  wavenumbers  is  believed  to  be  small, 
e.g.,  30-50  dB  below  that  of  the  convec¬ 
tive  peak,  accurate  estimates  are  essen¬ 
tial  to  an  understanding  of  turbulence  and 
Its  effect  on  the  structural  stability  of 
non-rigid  boundaries.  Thus,  an  estimator 
that  performs  well  primarily  In  the  peak 
region  of  the  TBL  spectrum  Is  not  sat i s- 
factory:  only  one  that  performs  uniformly 
well  over  a  (possibly)  large  dynamic  range 
will  adequately  solve  the  TBL  measurement 
problem. 

The  measurement  problem  is  all  the 
more  difficult  because,  in  this  case, 
estimation  of  the  spectrum  cannot  be 
translated  Into  a  parameter-estimation 
problem.  The  theoretical  models  are  pri¬ 
marily  qualitative  descriptions  of  the 
tru *  wavenumber  spectrum  and  thus  could 
bio.  estimates  based  on  inaccurate 
par  imetrlc  forms.  In  addition,  the  specu- 
i  at  ve  nature  of  the  models  prevents  their 
use  as  vehicles  for  hypothesis-testing 
formulations  of  the  spectra  I -est I mat i on 
problem.  At  best,  the  models  can  serve  as 
test  cases  for  evaluating  the  performance 
of  a  given  estimator;  at  worst,  they  can 
lead  to  erroneous  conclusions  about  an 
estimator’s  performance  over  a  more  gen¬ 
eral  class  of  TBL  spectra. 

Several  other  aspects  of  the  TBL 
measurement  problem  Increase  the  diffi¬ 
culty  of  determining  an  appropriate  spec¬ 
tral  estimator.  As  is  often  the  case  In 
array  processing,  physical  factors  prevent 
the  use  of  large  arrays  so  that  the  class 
of  possible  estimators  is  restricted  to 
those  that  perform  we  1  based  on  very  *>w 
samples  of  the  cross- spectrum.  Finally, 
the  uncertainty  with  respect  to  spectral 
smoothness  In  the  1 ow-wavenumber  region 
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Source-bearing  estimators  based  on 
the  projection  of  the  SCM  onto  the  signal 
subspace  have  proven  to  yield  excellent 
performance  in  many  signal /noise  environ¬ 
ments.  The  success  of  these  estimators 
depends,  to  a  large  extent,  on  the  eigen¬ 
values  of  the  SCM:  only  when  these  break 
down  into  two  distinct  sets  can  the 
covariance  space  be  decomposed  accurately 
Into  signal  and  noise  subspaces.  The  SCM 
for  the  Chase  model,  however,  has  no  such 
clear-cut  distinction  between  different 
portions  of  the  elgenspace;  the  eigen¬ 
values  are  distributed  almost  uniformly  In 
dB  over  a  50  dB  range.  Nevertheless,  a 
subspace  formulation  of  the  MLM  can 
Improve  performance  over  a  restricted 
region  of  the  spectrum. 


The  amplitude  spectra  of  the  SCM 
eigenvectors  reveals  that  each  eigenvector 
is  dominant  over  a  relatively  restricted 
portion  of  the  wavenumber  spectrum.  Thus, 
one  method  of  improving  performance  at  low 
wavenumbers  Is  to  project  the  SCM  onto  the 
space  spanned  by  those  eigenvectors  with 
beam  patterns  having  main  lobes  away  from 
the  convective  peak.  Since  the  spectral 
density  is  greatest  In  the  vicinity  of  the 
convective  peak,  this  projection  is 
equivalent  to  projecting  the  SCM  onto  the 
a  smallest  eigenvectors.  Let  z.  denote 
the  l-th  eigenvector,  ranked  in  order  of 
descending  eigenvalues  (X  ).  Then  the  MLM 
based  on  the  subspace  1  spanned  by  the  a 
smallest  eigenvectors  can  be  written  as*. 

P(k)  =  (  £h(k> 

where  £(k)  is  the  direction  vector  for 
wavenumber  k.  This  formulation,  in 
effect,  throws  out  what  would  be  con¬ 
sidered  the  "signal  space"  In  the 
bearing-estimation  problem  and  estimates 
the  spectrum  from  the  "noise  space". 


The  simulations  shown  In  Fig.  4  are 
the  result  of  projecting  the  SCM  onto  the 
subspace  spanned  by  the  8  smallest  eigen¬ 
vectors  for  a  12-sensor  array  based  on  the 
Chase  spectrum.  These  results  show  that 
projecting  the  SCM  onto  the  appropriate 
subspace  Improves  the  performance  of  the 
MLM  over  particular  portions  of  the 
wavenumber  spectrum.  The  demonstrated 
success  of  this  "crude"  subspace  version 
of  the  MLM  suggests  the  use  of  more  gen¬ 
eral  weighting  schemes  of  the  elgenspace 
to  optimize  performance  over  restricted 
regions  of  the  wavenumber  spectrum. 


CONCLUSIONS 


The  MLM  is  superior  to  AR  for  TBL 
spectra  with  respect  to  spectral  fidelity, 
smoothness,  and  statistical  stability.  An 
added  advantage  of  MLM  over  other  classes 
of  estimators,  particularly  when  little  Is 


known  about  the  true  spectrum.  Is  that  it 
always  forms  an  upper  bound  of  the  true 
spectrum  [7].  Nevertheless,  the  optimality 
of  the  MLM  or  of  its  subspace  version  has 
certainly  not  been  demonstrated.  Our 
efforts  on  the  TBL  measurement  problem 
continue  to  be  dominated  by  a  search  for 
estimators  which  are  robust  statistically 
and  are  not  overly  sensitive  to  modeling 
errors. 

The  primary  purpose  of  this  paper  was 
to  draw  attention  to  the  growing  need  for 
spectral  estimators  designed  for 

wavenumber  spectra  that  are  typically 
smooth  rather  than  peaked.  This  problem 
is  particularly  acute  in  the  array  case 
where  the  number  of  spatial  samples  is 
very  small.  The  design  of  such  estimators 
must  be  considered  with  some  care.  Wren 
spectral  peaks  are  not  desired,  estimators 
such  as  low-order  AR  may  Introduce  inap¬ 
propriate  peaks  in  the  spectrum  and  lead 
to  false  conclusions  about  the  true  spec¬ 
trum.  Frequently,  the  processes  under  I y- 
I ng  the  random  field  are  poorly  under¬ 
stood,  as  is  the  case  In  turbulent  boun¬ 
dary  layer  research.  These  cases  pose, 
perhaps,  the  most  difficult  proD'erm  In 
spectral  estimation  in  as  much  as  the 
estimator  must  be  designed  around  a  broad 
class  of  possible  spectra;  performance  for 
one  or  two  exemplars  of  such  a  class  may 
not  be  indicative  of  that  for  the  entire 
c I  ass . 
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Fig.  1.  C8,  RLR,  and  AR  estimate*  of  it*  tbl  wevenue«>er  spectrjs 
#re  shown  for  two  theoretical  models  o(  the  spectrus:  Chase  (left) 
and  Cor cos  (r  Ight). 
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Fig.  2.  951  confidence  intervals  bated  on  Rente  Carlo  » i me  I  at i on* 
are  shown  f or  AR  (left)  and  for  HR  (right).  The  simulation*  were 
based  on  24  snapshot*  from  a  ! 2-sensor  array  and  used  the  Che** 
model  of  the  TBL  spectra.  The  so' id  !  in#  denotes  tne  average 
spectra. 
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Fig.  3.  Performance  of  the  RLR  (solid)  is 
compared  with  the  Chase  spectrie*  (dotted) 
for  mlnlim**  redm dancy  (RR).  sensor-equi¬ 
valent  uniform  (U5),  and  aperture-equl va¬ 
lent  uniform  (UIO)  arrays.  All  spectra  are 
based  on  the  exact  SCR. 


Fig.  4.  Three  simulation*  (solid)  of  a 
tub space  formulation  of  the  RLR  are  shown 
for  the  Chase  model  (dotted)  of  the  TBl 
spectres. 
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penalizes  estimators  that  introduce  dis¬ 
tortion  i  rt  the  form  of  peaks  or  ripples  In 
the  spectrum. 

EVALUATION  Of  THREE  SPECTRAL  ESTIMATORS 

In  an  attempt  to  Isolate  properties 
that  are  essential  to  the  TBL  measurement 
problem,  the  performance  of  three,  we  11  - 
known  estimators  [I]  was  studied.  Evalua¬ 
tions  of  conventional  Oeamforming  (CB) 
with  a  triangular  taper,  Auto-Regressive 
Modeling  ( AR )  and  the  Maximum  Likelihood 
Method  ( ML M )  were  based  on  spectral  fidel¬ 
ity  and  sensitivity  to  statistical  error 
using  the  Chase  [2]  and  Cor cos  [3]  models 
of  the  TBL  spectrum.  These  models  were 
chosen  for  their  extremes  In  dynamic  range 
and  therefore  provide  a  means  to  assess 
per  f  o'-rr.anc  e  across  the  class  of  spectra 
that  .ary  In  rolloff  below  the  convective 
peak  . 

Spec  t  r  a i  f i de i ity 

The  spectral  fidelity  of  the  estima 
t'.rs  car-  be  determined  by  c  jmparing  the 
<mputed  ( so  1  I d  lines)  against  the 
theoretics'  (dotted  lines)  sce'tra  in  the 
two  pane's  of  fly.  i.  T  Or  Purposes  of 
i  i  us*  >  at  i  on ,  a  i  '  spec  ha  ait  based  on  the 
e-av*  spa* 'a'  cOvaria'ce  matrix  (5CM)  at  a 
y  '  ..  wp  t  ■  eo  jjnc  y  for  a  t  we  I  ve  -  sensor  , 

1  >near  tr  (form  array. 

The  p.  »or  f  I  de  l  1 1  y  of  v.  B  at  1  ox 

wa  v  e  ..umber  s  s  nut  Su'  pr  s  '  ny  g  •  ven  the 
eaxjye  of  power  f 'on  * ne  convective 
'ey  -  w*  '  h  is  expected  tr  .rr.  a  twelve- 
s  w  n'if  t  p  irnf./iier  ;  it  is  ■  lea'  f r  om  the 
s  '  jn t  •  f  f  ei  en<  es  het  xwer  the  t  wo  c  om  - 

Put  e--<  pe  *  r  a  t*  at  e-  eryv  -ejr  the  spec 

»  -  a  i  pear  1  Km  •  r  at  es  t  he  '  b  spe'  *  '  om  at  a  i  i 

wa  -  -■ '  .-it  e  r  s 

xp.p  w  3  w  t  *  t  *  e  b  fe  its,  h.  M  a*  d 
AR  p.r  vide  T.jd  better  *  M  to  the 
i  hno'  et  a  I  spe  t'a.  Thr  super  'or  dynamo 


linear  uniform  array. 

flfl.  2  shows  95%  confidence  Intervals 
(solid  II  nes  )  a  I  ong  with  t he  average  es t,  i  - 
mat  or  performance  based  on  a  12x12  SCM  for 
AR  (left  panel)  and  for  ML M  (right  panel). 
The  results  are  consistent  with  known  sta¬ 
tistical  properties  of  the  estimators 
[4,5];  the  confidence  intervals  for  AR  are 
substantially  greater  than  those  for  ML M . 
This  greater  sensitivity  of  AR  to  statist¬ 
ical  error  is  often  compensated  by  Its 
higher  resolution  in  those  cases  where  the 
spectrum  is  dominated  by  narrowband  peaks 
1 5  ]  .  The  present  results  indicate  that  AR 
loses  Its  advantage  over  the  statistically 
more  stable  ML M  for  smooth  spectra  such  as 
the  TBL  spectrum. 

EXTENSIONS  Of  THE  MLM 

The  good  spectral  fide1 ity  and  rela¬ 
tive  insensitivity  to  statistical  error 
exhibited  by  the  MLM  for  TBL  spectra  is 
encouraging  and  raises  the  possibility 
that  better  performance  may  be  achieved 
through  extensions  of  the  MLM.  We  report 
on  two  attempts  at  improving  the  perfor¬ 
mance  of  the  MLM,  particularly  in  the 
l ow -wavenumber  region. 


Non  uniform  an 


geometr 1 1 


Non-uniform  array  geometries  have 
produced  Improved  source-bearing  estimates 
by  increasing  the  array  aperture  for  a 
fixed  number  of  sensors.  A  non-uniform 
formulation  of  the  MLM  based  on  a  Minimum 
Redundancy  array  [6]  of  five  sensors  pro¬ 
duced,  somewhat  unexpectedly,  a  degrada¬ 
tion  in  performance  from  that  achieved  by 
a  uniform  array  of  five  sensors  for  the 
Cr.ase  and  Corcos  spectra.  These  results 
ar e  shown  in  f ig.  3  for  the  Chase  spectrum 
a 1 ong  „lth  the  omputed  HpM  spectrum  for  a 
10  sensor .  linear  array  that  has  the  same 
aperture  as  the  MR,  five  sensor  array. 


»  -jr  ye  of  AR,  *-  :,»vpr  ,  *  s  a  s.  j  near  '  ed 

by  sift'll  r  •  p  p  ’  I  r.y  »'•!)  t  *  .,s  umpr  urn i  ses 
1'  v  -.rr..  S  C'  i*er  '  in.  ’O  a.n'eve 

a  r  j r  ji  f  ,  smrx.a  ’  tits  tr  s  r  *>  '  f)L  spe  t  r  a  . 

'  s.e  u  )e'  <  *  r  hr  AP  *.■  )e  i  siust  be 

•nr  'eased  ~e  I  be,'_.nd  t'e  '  iist.X'  of  sensor 
e  enie*  t  s  ;  1  owe  r  order  rt.1  ■■  :e  ■  s  p  ‘  d'CC  ever 

gr  ea*  r"  d  s  t  or  *  '  cm  i  t  r-.e  *  ■  ’  •  jf  r  i  pi  - 

p  1  •  *  g  ■ 

5*  it.  stjuj  1  er  (o' 

T  r  ix  '  leva'  '  s  o'i  s  i  ri  -hr  p  'ey  s  u  S  '  jt 

s t  w >r i  were  Eased  •  »n  t. ' .e  e ,  u '  '  .  a  r ns  * 

t  ne  v  M  .  T ....  determ  ,  'e  t'e  e»fe  r  ,  c  f  s '  a 

t  i  st  '  i.  a  i  error.  t'e  per  f .  ■>  ma  re  of  ea  .  *• 

est  ina'or  was  st  orbed  by  model  i  r>g  *  ne 

snapshot.  -  are' aged  es'  ■  rrv*t  e  o*  tne  V.M  as  a 
c  i jmp  1  ex  -  W  I  Shar  t  r  andom  mat  rix  and  per  form- 
i  ng  Monte  Car  I  o  s  a*  onr.  Results  for 

the  Chase  spec  t '  urn  a’ e  r  epx  >r  t  ed  and  are 
based  on  ,4  snapshots  from  a  I*.’  sensor. 


The  apparent  discrepancy  between 
sour ce  Dear i ng  and  the  TBL  results  can  be 
arrounted  for  by  differences  in  the  beam- 
lotterns  of  the  MR  and  uniform  arrays. 
a’ though  MLM  does  not  depend  directly  on 
*  he  array  beampattern,  as  is  the  case  for 
■  ir.pnt i ona I  beamforming.  It  nevertheless 
mo  i ds  the  optimum  steering  vector  out  of 
the  array  beampattern.  E -ami  nation  of  the 
t'eampat  t  er  ns  for  the  two  arrays  reveals 
narrower  ma r n I  one  and  higher  s i de 1 obe  lev¬ 
y's  in  the  MR  a-ray  than  for  the  uniform 
array  which  would  enhance  the  abl I ity  of 
the  HR  array  to  suppress  components  near 
•re  target  wavenumber  but  reduce  Its  abl'- 
•ly  to  suppress  components  far  away.  This 
unde’ scores  the  fact  t  fiat  resolution  and 
spectral  fidelity  may  t-  competing  cri¬ 
teria,  depending  on  the  nature  of  the 
underlying  spectrum. 
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APPENDIX  B 

INTERNAL  MEMORANDA  ON  A 
HAVE  GUIDE  MODEL  OP 
BOUNDARY  LAYER  PRESSURE  FLUCTUATIONS 
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V.  Estimate  of  the  Low  Frequency  Power  Spectrum 

of 

Wall  Pressure  Fluctuations 
by 

Robert  F.  Lambert 
Institute  of  Technology 
University  of  Minnesota 


Introduction 

Up  to  the  present  most  descriptions  of  wall  pressure  fluctuations  in  wind 
tunnels,  water  tunnels,  and  the  like  have  been  based  on  two  point  space-time 
correlation  measurements.  Regardless  of  experimental  techniques  whatever  facts 
we  know  about  these  fluctuations,  say  power  spectra,  must  be  gleaned  from 
experimental  observation.  The  best  information  to  date  which  comes  from  wind 
tunnel  data  employing  flush  mounted  pressure  sensitive  transducers  of  various 
sizes  lies  in  a  restricted  normalized  frequency  range 

0.02  <  rf^-  =  Z  <  10  (1) 

*  do 

where  is  the  free  stream  velocity  of  the  fluid  and  6*  is  the  displacement 
boundary  layer  thickness.  Power  spectra  measurements  at  lower  frequencies  have 
been  contaminated  by  acoustic  noise  which  can  dominate  the  wall  pressures  and 
hence  are  not  deemed  valid. 

Nevertheless,  we  would  like  reasonable  estimates  of  the  wall  pressure  power 
spectra  due  to  turbulence  at  lower  frequencies,  i.e.  <3  <  0.02  in  order  to  design 
experiments,  predict  the  performance  of  new  electronic  instrumentation  designed 
to  make  pressure  measurements  in  low  noise  facilities,  and  make  other  acoustic 
power  predictions. 

The  work  here  presented  describes  one  approach  to  that  low  frequency 
estimate.  It  is  a  semi-empi ncal  approach  based  on  experimental  information 
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about  the  space-time  correlations  and  power  spectrum-gathered  at  higher 
frequencies  together  with  an  estimate  of  the  low  frequency  characteristics  in 
the  spectra*  obtained  by  viewing  the  turbulent  boundary  layer  as  a  sort  of 
convective  waveguide  of  finite  width.  It  Is  called  a  "convective  waveguide 
model".  This  information  in  turn  can  be  used  to  estimate  a  half-power  bandwidth 
and  center  frequency  for  the  spectrum.  In  order  to  estimate  the  spectral  level 
at  and  below  the  low  frequency  cut-on  for  boundary  layers  of  finite  thickness  a 
wave  number  translation  is  hypothesized.  The  model  also  is  useful  in  estimating 
the  wavenumber/ frequency  pressure  spectrum  at  low  frequencies  and  low 
wavenumbers  at  low  mach  numbers. 

Background 

In  order  to  make  the  extension  to  lower  frequencies  a  normalized  space-time 
correlation  function  RnU.n.T)  and  a  normalized  power  spectrum  ♦n(o»)  will  be 
utilized.  They  are  related  by  the  expression 

♦n(“)  “  J  Rn(0.0.Oe1wTdx  (2) 

_  ao 

where  t  is  a  time  delay  and  and  n  are  longitudinal  and  lateral  spatial 
separation  variables,  respectively. 

Although  many  phenomenological  models  of  Rn  have  been  proposed  to  represent 
experimental  observations  none,  even  the  most  elaborate  or, as ,  have  been 
entirely  successful.  Since  no  reliable  low  frequency  data  now  exist  only  a  very 
simple  model  for  *n(u>)  will  be  here  proposed.  Later,  a  more  elaborate  model  can 
be  developed  if  the  data  warrant. 

Accordingly,  the  trial  function  for  Rn  will  be 
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RnU,n,t)  -  exp  [-  M-^MDexpCxtU-VcTl]  (3) 

vce 

which  yields  a  low  pass  power  spectrum  that  is  deemed  valid  at  high  frequencies 
only.  Here  Vc  Is  a  mean  convection  speed  and  o  Is  a  mean  eddy  lifetime. 

Before  proceeding  a  few  remarks  about  the  above  functional  form  for  R,,  are 
In  order.  When  the  above  function  is  compared  with  various  experimental  obser¬ 
vations  several  points  can  be  made.  First,  when  R^o, n,o)  versus  n/6*  Is  fit  to 
most  wind  tunnel  data  the  parameter  <2  *  2/6*.  Second  vrfien  Rf,U,o,n/Vc)  is  fit 
to  longitudinal  data  one  finds  that 

0  V./6*  *  18 

and  (4) 

Vc/V.  *  0.8 

at  least  for  i/6*  sufficiently  large,  i.e.  £>O.8V.0.  Further,  most  data  exhibit 
sharp  peaks  in  I^U.o.S/Vc)  for  close  observation  points  only  and  the  curves 
become  broader  as  i  increases.  Also,  Vc  is  not  strictly  constant  but  varies 
from  0.58V*  for  very  close  observation  points  and  gradually  approaches  0.80V. 
as  1;  is  increased.  However,  when  one  fits  power  spectra  data  to  measured  values 
a  "best  fit"  usually  obtains  for  Vc  =  0.8V«.  Eq.  3  on  the  other  hand  exhibits 
sharp  peaks  in  RpU.o.S/Vc)  for  all  i.  Finally,  there  are  discrepancies  in 
absolute  pressure  levels  of  roughly  2  to  3  dB  due  to  finite  transducer  size 
effects  at  higher  frequencies  which  is  deemed  not  so  serious  as  discrepancies  in 
convective  velocity  predictions.  The  ratio  of  the  rms  surface  pressure 
<p2(x,y ,t)>^  to  the  dynamic  velocity  head  1  n  p0V„2f  where  p0  is  the  mean  fluid 
density,  is  a  constant  whose  value  can  be  determined. 

To  introduce  the  convective  waveguide  model  I  refer  to  the  Kolmogoroff 
theory  of  turbulence  which  asserts  that  energy  for  sustenance  of  the  fluctuating  " 

I 

I 
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part  of  the  flow  Is  given  directly  to  a  group  of  eddies  (here  called  energy¬ 
bearing  eddies)  In  a  rather  definite  range  of  sizes  through  Inertial  transfer 
from  the  mean  flow.  A  very  complete  discussion  of  this  process  can  be  found  In 
Batchelor*  and  the  theoretical  and  experimental  results  here  utilized  are  taken 
from  chapters  6  and  7  of  his  book. 

On  the  other  hand,  energy  is  removed  from  the  turbulent  flow  at  the  same 
rate  by  viscous  dissipation  which  Is  strongest  for  eddies  whose  scale  Is  much 
smaller  than  the  scale  size  Le  of  a  typical  energy  bearing  eddy.  The  range  of 
eddy  sized  L  «  Lg  Is  called  the  universal  equilibrium  range  because  the  motions 
of  these  eddies  1$  uniquely  determined!  on  a  statistical  basis  by  only  two 
parameters,  the  energy  per  unit  mass  Eg  supplied  by  the  fluctuating  flow  and  the 
kinetic  viscosity  v  of  the  fluid. 

It  Is  possible  to  relate  the  sustaining  energy  supplied  to  the  energy 
bearing  eddies  to  a  peripheral  velocity  Ve  and  a  scale  length  via  the 
relation 


Then  this  theory  predict  a  time  scale  for  the  energy  bearing  eddies  of  roughly 


In  order  for  Eqs.  5  and  6  to  be  valid,  it  is  required  that  the  inertial 
forces  predominate  over  the  viscous  forces,  hence  the  Reynolds  number 

VpLp 

Re  -  »  1  (7) 

for  a  universal  equilibrium  range  to  exist.  The  mean  value  of  the  rms 
fluctuating  velocity  Ve  is  approximately  0.05  V.. 


These  numbers  yield 
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ee  * 


*  20 


Jil 

v. 


(8) 


which  Is  In  close  agreement  with  the  result  of  Eq.  4  and  which  In  turn  were 
based  on  the  correlation  function  of  Eq.  3.  Hence,  we  readily  Interpret  the 
parameter  0  In  Eq.  3  as  a  "mean  eddy  lifetime".  The  parameter  iq  in  Fig.  3  that 
best  fit  most  experimentalist  data  in  wind  tunnels  is  <i  a  1/6*  and  hence  xj  is 
Interpreted  as  an  "inverse  mean  eddy  size". 

By  taking  the  Fourier  transform  of  the  normalized  trial  temporal  correlation 
functions 

Rn*(0,0t)  =  exp[-<1Vc  |t  |]  (9) 

the  normalized  power  spectrum  becomes 


2  »h*n*(*) s  v . 2 '  T  u°) 

vc 

where  wh  *  Vc/6*  is  interpreted  on  the  upper  cut-off  frequency  of  the  power 
spectrum  and  numerically  is  in  excell  ant  agreement  with  high  frequency  spectral 
data.  The  normalized  frequency  a  is  given  by  Eq.  1. 

However,  at  the  low  end  of  the  power  spectrum  there  is  much  more  uncertainty 
and  ambiguity.  By  viewing  the  turbulent  boundary  layer  as  a  convective  wave 
guide  we  will  set  the  scale  of  the  largest  energy  bearing  eddies  as  Lm  *  6  where 
o  is  the  boundary  layer  thickness.  For  smooth  boundaries  *  =  86*  on  the  basis 
of  the  7th  power  law  velocity  profile.  If  Le  is  the  scale  size  of  the  eddy  and 
Ve  is  the  peripheral  eddy  velocity,  then  the  time  it  takes  for  the  eddy  to  cycle 
once  is 

il  Le 

Tp  —  =  iiup  (11) 

Ve 
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where  from  Eq.  4  we  take  ee  *  186*/V«.  Then  for  the  largest  energy  bearing 
eddies  their  characteristic  frequency  will  be 


or 


U»A 


2fl 

T* 


2V»  = 
186** 


0.1111— 

6* 


3*  =  0.1111. 


(12) 

(13) 


Thus,  we  take  u»A  as  the  estimate  of  the  lower  frequency  half-power  point  in  the 
power  spectrum  for  incompressible  fluid  flow. 

As  an  aside  it  is  interesting  to  observe  several  other  matters.  The  scale 
of  the  energy  bearing  eddies  now  should  be 

6  >  Lq  >  6*/8 

while 

^cmitl  ~  a  o*)  (14) 

Thus,  the  center  of  mean  energy  bearing  eddy  is  located  an  average  distance 
3/2  6*  from  a  smooth  wal 1 .  By  the  same  token 


''emit,  ’  0.58V.  »  V(^  »*)  (16) 

which  is  characteristic  of  the  smallest  eddies  in  the  energy  bearing  range. 

Note,  they  are  much  closer  to  the  wall. 

We  also  speculate  that  any  variation  in  Vc  with  spatial  separation  as 
observed  by  various  investigators  can  be  attributed  to  the  more  rapid  decay  of 
energy  in  the  high  frequency  eddies.  Since  these  are  small  eddies  that  lie 
close  to  the  boundary  they  tend  to  lower  the  average  drift  velocity  as  observed 
at  small  separations.  However,  at  greater  separation  distances  say  i  >  0.8VJ3 
or  C/6*> 14  only  the  low  frequency  eddies  retain  any  coherence  from  point  to 
point.  This  observation  also  may  account  for  any  broadening  of  the  peaks  in  the 
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space-time  correlations  at  larger  separations.  Thus,  we  are  aware  of  variations 
In  the  average  drift  velocity  with  sensor  element  spacing  which  could  be  a 
troublesome  point  when  interpreting  coherence  data  obtained  from  an  array  of 
sensors  spaced  a  finite  distance  apart. 

Low  frequency  estimate  of  *n(w)  based  on  the  convective  waveguide  model. 

The  general  shape  of  the  normalized  measured  power  spectrum  for  <3>0.Q2  is 
shown  in  Fig.  1.  We  note  that  it  peaks  at  some  frequency  up  and  exhibits  a  high 
frequency  cut-off  (-3  dB  point)  at  The  drop-off  date  for  a»©h  is  _g  per 
octave.  However,  frequencies  below  flp  the  normalized  spectra  does  not  decrease 
rapidly  with  decreasing  frequency  but  appears  to  approach  a  constant  for  0<<0p 
characteristic  of  some  static  pressure  level  associated  with  the  fluid  flow. 

A  number  of  attempts  have  been  made  by  others  to  reproduce  this  spectrum 
theoretically  at  low  frequencies  but  all  require  a  large  number  of  adjustable 
parameters.  Moreover,  these  parameters  are  difficult  to  interpret  physically 
and  just  as  difficult  to  measure  experimentally.  A  somewhat  different  approach 
will  be  presented  next  which  seems  to  incorporate  most  of  the  salient  features 
of  «n(u)  and  hopefully  is  accurate  at  low  frequencies,  well  below  <3p.  Also,  any 
parameters  that  control  <t>n(u)  for  u<up  should  be  readily  measureable,  preferably 
independent  of  spectral  measurements. 

The  first  task  is  to  incorporate  the  estimate  of  5^  Eq.  12  into  the  nor¬ 
malized  power  spectra.  For  this  purpose  the  convective  waveguide  will  be  viewed 
as  a  bandpass  filter  and  a  typical  low-pass  to  band-pass  frequency  transfor¬ 
mation  of  Eq.  10  will  be  utilized.  However,  that  transformation  by  itself  will 
not  be  sufficient  since  the  estimated  spectra  will  drop  off  too  rapdily  below 
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flp.  In  order  to  Insure  a  proper  low  frequency  asymptotic  behavior,  a  cut-on  In 
wave  number  it  some  low-frequency  u»0  will  be  hypothesized  for  a  turbulent 
boundary  layer  of  finite  depth.  We  could  estimate  a0  *  3  and  this  may  well 
suffice.  However,  Inltally  I  would  prefer  to  sa>  that  the  cut-on  occurs  when 
the  perimeter  of  the  largest  eddy  permitted  by  the  guide  2*a  equals  the  eddy 
wavelength  Vc/f0  where  f0  =  u»0/2w .  Here  2a  is  the  physical  width  of  the 
waveguide. 

Thus,  we  have  for  the  cut-on  wavelength 


f- *  2*a 
To 


yielding  a  cut-on  wave  number 


kf>  =  U  =  . 


a) 


b) 


(16) 


The  convective  waveguide  now  will  be  dispersive  at  low  frequencies  and  the 
convective  ridge  obeys  the  relation 

l£--UJ  k  >  kg 


s  k^Vr2-wo^  k  >  kn  a) 


or 


(17) 


Note  for  the  sufficiently  large  u>,  k  ♦  .  However,  as  u>  ♦  0,  k  ♦  and 

*  c  vc 

(18) 


V. 

_  .  0 

vc  w  +  h  “ 


Thus,  the  power  spectrum  must  be  modified  to  agree  with  Eq.  17  which 
requires  that  the  normalized  frequency  3  be  transformed  to  S2+(-j*-)2  *  32+u2 
where  o<u<l  since  0<o**u. 

Before  proceeding  with  the  frequency  translation  required  by  Eq.  17  let  us 
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examine  the  phase  velocity  of  spectral  components  defined  as 

Vp  5  £  -  Vc//1  ♦  -o 2/J  (19) 

Clearly,  Vp  cannot  exceed  the  convection  speed  of  the  fluid  Vc  but  will 
approach  Vc  as  Note,  Vp*0  as  w-*-0  In  this  formulation.  More  about  the 

behavior  of  Vp  In  a  subsequent  tract. 

Further^the  group  velocity  Vg  defined  as 


(20) 


fas  characteristics  that  enable  it  to  exceed  Vc,  and  below  a  critical  frequency 
can  exceed  the  speed  of  sound  c. 

Before  making  the  required  frequency  transformation  a  few  definitions  are 
required.  First,  we  define  a  half-power  bandwidth  for  the  turbulent  wall 
pressure  spectrum  as 


<*>b  =  “*h  "  w£  (21 ) 

where  *  0.8  V„/6*  and  u»A  *  0.1111  V./o*.  Also  we  define  a  band  center 
frequency 

wc  =  f  (22) 

When  these  parameters  are  normalized 


Wb  *  0.689 


a) 


(23) 


and  uic  =  0.298  b) 

Note,  that  the  predicted  spectrum  actually  will  peak  somewhat  below  due  to 
the  frequency  translation. 

Following  standard  procedures  we  now  rewrite  Eq.  10  in  factored  form 


1 

2 


wh  *n(w) 


rv  "| 

rv-  i 

7"  (i  -)  ♦  i 

7"  (-i  -»)  ♦  i 

[_  c  J 

L  c  J 

(24) 


and  then  replace 
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and 


V  X2  -  (A2  ♦  U2) 

L y  )  ( i  a)  by  — = - 

i  / 132  +  u2  aK 


<2  2  -  (<22  +  u^) 


a) 


(25) 


ITTJM  «)  by  - - 

l~l/  w  + 


b) 


2  \  - 
M  )  u)b 


in  order  to  obtain  a  normalized  spectrum  that  reduces  to  Eq.  10  for  5j2  »  u2  and 

has  the  requisite  low  frequency  behavior  for  uS2  <<  w2. 

The  transformed  spectrum  now  becomes 

<2b2  (“^  + 

5  9~  •>  T. 

.4 


(26) 


ry  0*5  o  0  O  O 

(w  +  u  )  ♦  (a  -  2SC  )(a  +  u  )  + 

which  exhibits  the  proper  high  and  prescribed  low  freqency  behavior.  Note  that 

1  <2u2 

7  wh  *n  (“)  ■-  ~2~  •  (27) 

a  ♦  ®  a4- 

Thus  at  high  frequencies,  a2  >>  wb2,  £2  >>  u2f  Eq.  27  predicts  the  -6  dB  per 

doubling  of  frequency  in  agreement  with  most  reliable  wind  tunnel  data.  The 
spectrum  peaks  at  ap  =  7ac2  .  y2  as  required  but  will  not  exhibit  a  -3  dB  point 
below  2C  unless  u2  is  sufficiently  small. 

At  low  frequencies  u2  <<  wc2 ,£2<<u2 ,  Eq.  26  reduces  to 


2  wh  ^ n  (y) 


X 


u2  -c2  2 

1  +  .  2  ^  '  u2i 

uib  u 


(28) 


which  predicts  a  finite  normalized  level  (less  tha  0  dB)  for  finite  values  of  u. 
Note,  the  power  spectrum  will  be  flat  for  a<ac  in  the  special  case  u>c=m. 

It  is  convenient  and  plausible  to  approximate  a  by  6 ,  the  boundary  layer 
thickness,  on  the  basis  of  the  convective  wavegui ie  theory  applied  to  a  wind 
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tunnel  or  Mater  tunnel.  If  a  «  6,  then  u  *  6*/ 6  *  0.125  for  a  1/7  poMer  shear 
velocity  profile.  Note,  numerically  that  u  *  which  probably  Is  not 
fortul  tous* 

Eq.  28  predicts  at  low  frequencies  that 

10  log10Lu.h  *nM25)J»  "2.4  dB  (29)  . 

in  excellent  agreement  with  reliable  low  frequency  wind  tunnel  spectral  data, 
see  Fig.  1.  Most  wind  tunnel  data  level  off  between  -0.5  and  -2.5  dB  with  about 
a  ±  1/2  dB  spread. 

Note,  as  u  ♦  0,  the  value  of  the  low  frequency  asymptote  «n*(u)  becomes 
progressively  lower.  For  y=0 

^  “h  *n(a)  +  7^  A2  (30) 

3c 

A  plot  of  Eq.  30  also  is  shown  in  Fig.  1  for  comparison  purposes. 

A  few  remarks  about  the  experimental  curves  are  in  order.  Normalized  power 
spectrum  measurements  of  M.  K.  Bull  performed  at  Southhamptom  in  1963  in  a  low 
noise  wind  tunnel  are  shown  by  the  shaded  area.  This  normalization  was  done 
by  Dr.  0.  H.  Tack  in  1965  who  also  attempted  to  reproduce  the  spectrum  theoretically 
using  other  techniques,  namely  a  modification  of  Eq.  3  which  agreed  more  closely 
with  Bull's  spatial  correlation  data.  Bull's  data  still  are  regarded  as  among 
the  best  and  most  complete  available  for  wind  tunnels.  His  pressure  sensors 
were  miniature,  0.76  mm  diameter  lead-zi rconate-titinate  discs  flush  mounted  on 
the  wall.  His  experimental  parameters  were  flow  speeds  from  100  to  164  m/s  and 
displacement  thickness  between  1.5  and  4.4  mm.  The  individual  measurement 
points  are  not  shown,  only  the  range  of  spread  after  normalization.  The  data  do 
normalize  exceedingly  well  in  the  range  0.02  as  shown  by  the  shaded  area. 
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Also  shown  for  other  comparison  purposes  are  low  frequency  wall  pressure 
fluctuation  spectra  obtained  from  glider  experiments.  These  data  do  exhibit 
rather  wide  scatter  but  level  off  at  about  10  dB  below  their  peak  at  low  fre¬ 
quencies  whereas  Bull's  data  leveled  off  considerably  higher.  If  the  asymptote 
Is  taken  to  be  -10  dB  re  maximum,  then  from  Eq.  28  u  3  0.042  utilizing  the  same 
scaled  values  of  fflb  and  oc,  requiring  a  somewhat  smaller  boundary  layer  ratio 
for  6*/a.  Currently  a  search  is  being  made  to  find  good  low  frequency  spectral 
data  utilizing  a  range  of  boundary  layer  thiccness  ratios. 

This  formulation  of  the  power  spectrum  do»s  provide  relatively  simple  quan¬ 
titative  relationships  between  the  normalized  power  level  and  the  boundary  layer 
profile  as  measured  by  6*/6,  Eqs.  26  and  28.  This  relationship  is  worthwhile 
exploring  in  more  detail  since  a*  and  6  can  be  measured  Independently  of  *(w). 
(See  appendix). 

However,  what  is  especially  appealing  about  this  formulation  Is  that  fairly 
accurate  values  for  both  u>b  and  u>c  are  available  from  high  frequency  experimen¬ 
tal  data  combined  with  a  simple  hypothesis  about  the  scales  of  the  energy 
bearing  eddies. 

Finally,  a  table  of  parametric  values  of  6*,  u>b»  “>£»  and  <*>c  ^or  the  water 

tunnel  facility  at  a  number  of  flow  speeds  ranging  from  3  to  18  m/s  is  attached. 
The  estimates  of  f^  =  ^  range  from  30  to  217  Hz.  Hopefully,  these  values  are 
within  the  range  of  good  spectral  measurements  in  our  facilities. 
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Appendix 

This  appendix  will  consider  several  Implications  of  the  results  and  may 
suggest  sone  additional  experiments. 

For  example.  If  one  assumes  a  ^  power  low  for  the  shear  flow  near  the  wall 
then  the  parameter 


* 


which  leads  to  the  result  that  y  +  0  as  n  ♦  *.  The  value  of  u  has  to  do  with 
the  shape  of  the  shear  boundary  layer  flow  profile  which  in  turn  is  related  to 
the  compressibility  of  the  flow,  the  geometry  of  the  wall,  the  surface  rough¬ 
ness,  etc. 

Also  it  can  be  shown  that 

n  s  «  v  A2 

f  00 


where  VT  is  a  friction  (shear  stress)  velocity  and  a  is  a  constant  that  can  be 
evaluated  numerically  under  certain  hydrodynamics  assumptions  and 
approximations. 

For  small  values  of  y  the  normalized  power  spectrum  approaches  the  value 

1  Vt  1  1 

2  uh  *n(M)  >  *2  (  —  J2  (.  —  -  —  J2  ,  S«y2  A3 

M  ♦  0  »  (3^  uJh 

revealing  that  VT  becomes  an  important  parameter  in  setting  the  low  frequency 
power  level.  These  predictions  suggest  that  controlling  y  via  boundary  layer 
shape  could  be  an  important  experimental  procedure. 

In  1970  Aupperle  and  Lambert  reported  in  a  study  supported  by  NASA  that  the 
surface  roughness  could  alter  the  low  frequency  power  level  by  as  much  as  lOdB, 
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the  rougher  the  surface  the  higher  the  level.  We  Mere  able  to  normalize  the 
effects  of  surface  roughness  on  the  pom er  spectrum  by  employing  the  parameters  a 
and  VT.  The  changes  In  measured  values  of  y  also  Mere  consistent  Mlth  the 
changes  In  observed  values  of  «n(Q,w).  However,  another  Independent  scale  para¬ 
meter  Is  required  to  normalize  the  frequency  scale  properly  for  rough  surfaces. 
This  Is  the  earliest  study  I  know  of  where  a  deliberate  attempt  to  roughen  the 
surface  was  made  to  enhance  the  effect  of  wall  pressure  fluctuations  at  low  fre¬ 
quencies.  The  signal  to  ambient  acoustic  noise  ratios  were  of  the  order  of  + 
25dB  over  most  of  the  frequency  range  of  interest.  The  studies  were  conducted 
down  to  200  Hz. 

As  far  as  I  can  tell  all  of  the  predicted  results  in  this  tract  are 
compatible  with  experiment  and  in  no  significant  way  conflict  with  theoretical 
power  spectrum  predictions  by  Blake  (1984),  Chase  (1980)  and  perhaps  Williams 
(1982).  In  the  situation  discussed  by  Williams  one  can  see  that  *n(u>)~u>2  at  low 
frequencies  i.e.  a2«fflp2  and  a2«y2,  only  in  the  limiting  case  y  0.  The 
glider  experiment  seems  to  be  as  close  to  that  situation  as  any  experiment  known 
to  me.  An  evaluation  of  Williams  work  will  be  made  later. 
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VI  Estimate  of  the  Low  Wavenumber/ frequency  Spectrum 
of  Wall  Pressure  Fluctuations 


Robert  F.  Lambert 


Introduction 

Recently,  in  tract  V,  I  proposed  a  new  method  for  estimating  the  normalized 


power  spectrum  «n(ui )  for  wall  turbulence  that  was  consistent  with 


experi¬ 


ment  especially  at  low  frequencies  and  had  a  fairly  simple  physical  interpreta¬ 
tion  in  terms  of  measureable  parameters.  The  turbulent  boundary  layer  was 
viewed  as  a  convective  waveguide  that  cuts-on  at  a  frequency  of  about  u0*Vc/6 
above  which  the  fluid  motion  could  be  regarded  as  incompressible.  Here  Vc  is 
the  convective  fluid  velocity  and  6  is  the  thickness  of  the  boundary  layer. 

Below  a  corresponding  wavenumber  the  pressure  fluctuations  are  governed  by 
compressible  wave  motion,  i.e.  sound  waves.  The  transition  wave  number  was 
estimated  to  be  k0  *  u>0/Vc  -  1/6  at  low  convective  mach  numbers. 

An  analytical  expression  for  a  normalized  power  spectrum  *n(w,u)  was 
developed  which  yielded  results  consistent  with  most  experiments  and  whose  shape 
was  controlled  by  three  parameters,  u»^=Vc/6*  =  0.8  Vw/6*,  =  V,„/96*,  and 

u=6*/6.  It  predicted  a  low  frequency  asymptote  that  was  dependent  on  these 
three  parameters  according  to 


the  relationship 
l/2u»h  *n(°.w) 


2  -  2 

U c  U)r  c  0 

1+— / 
a b  n 


where  wjj  and  u»c  are  related  to  and  Here  i3=u)6*/V08  is  a  dimensionless  fre¬ 
quency  variable.  In  the  limit  as  p  + o,  corresponding  to  perfectly  incompressible 
flow,  «n(w,o)  approaches  zero  as  u^asu  tends  to  zero.  The  general  results  seem 
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to  be  consistent  with  theoretical  predictions  by  other  workers  except  for  a 
recent  paper  by  Williams  (1982).  More  on  that  matter  later. 

The  purpose  of  this  tract  is  to  extend  the  convective  waveguide  scenario  to 
predictions  of  the  wavenumber/frequency  spectrum  PCk^.kg.u)  where  kj  and  k£  are 
longitudinal  and  lateral  wavenumbers,  respectively,  and  u>  is  the  radian  fre¬ 
quency  in  question.  To  do  this  from  first  principles  is  a  long  arduous  task 
and  rarely  done,  so  it  is  expedient  to  make  use  of  some  of  the  results  developed 
by  Williams  (1982)  and  then  adapt  them  to  a  convective  waveguide  model.  I  deem 
William's  general  results  to  be  valid  only  in  the  normalized  frequency  range 
( w6/ V»)^  >>  1.  Thu  does  not  exclude  the  wavenumber  region  k*w/c  so  some  of  his 
results  on  compressible  flow  can  be  utilized  directly.  The  emphasis  will  be  on 
developing  suitable  analytical  forms  for  the  wavenumber/ frequency  spectrum  valid 
for  low  longitudinal  wavenumbers  and  low  frequencies  at  low  mach  numbers. 
Analysis 

Briefly,  Williams  (1982)  derived  an  expression  for  the  wavenumber/frequency 
spectrum  by  a  modification  of  the  Corcos  model  so  as  to  "account  for 
compressibility  of  the  fluid  at  low  wavenumbers  while  retaining  the  apparently 
good  agreement  with  experiment  at  high  frequencies  and  at  higher  values  of  wave- 
number  where  the  spectrum  is  immune  from  the  influence  ot  compressibility".  In 
that  sense  his  contribution  is  valuable.  Without  going  into  any  detail  his  pre¬ 
dicted  wavenumber/ frequency  spectrum  has  the  general  form 


p<kl'k2-“>  *  “o2v3s3GoOfto[1+kl^Bo(k2t) 

xla0(5>2  *  al”2  *  a2H4*n<o)“J<tlt>2-H2:l1 


or 


P*(k1,k2,w)  =  Pq(T,u)  +  PjfT,*)  +  P2(k,u) 


a) 


b) 


(1) 


consisting  of  three  separate  additive  components.  In  his  view  the  component  P0 
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is  due  exclusively  to  incompressible  flow  while  Pj  and  ?2  which  are  mach  depen¬ 
dent,  M=V/c,  are  contributions  from  compressible  fluid  flow.  The  parameter 
60(<*»6/V)  is  the  frequency  spectrum  due  to  an  integrated  fluctuating  velocity 
source  function  that  tends  to  a  constant  as  ^  tends  to  zero.  While  I  have 
no  serious  qualms  about  that  particular  result  some  of  the  details  are  not 
clear  from  his  work.  The  functions  A<£  ]  and  BJ[  ]  are  longitudinal  and 
lateral  wavenumber  dependent  factors  that  “agree  with  experiment".  The  other 
parameters  a0.  a^,  a2»  are  scaling  constants  and  R  is  a  "large  radius  on  a  flat 
surface  over  which  the  turbulent  field  is  homogeneous".  The  function  U.X  ]  is  a 
singularity  function  that  arises  from  a  spatial  transform  of  the  lossless  wave 
equation  and  predict  a  singularity  in  pjtT.u)  at  k=u/c. 

Williams  proceeds  to  compute  a  power  spectrum 

oo 

*(“>)  =  JJ  P(k1,k2,u))dk1dk2  (2) 

—  00 

which  yields 

*(U)  =  POV3O02GO(-J“){a  +  tsM2  +  YM4*n(^)}  a)  (3) 

or 

«(w)  =  ®Q(w)  +  +  4>2  (  oj  )  b) 

where  <x,e,  and  y  are  constant  of  integration. 

Williams  results  then  reveal  that  all  three  terms  in  $(u>)  corresponding  to 
both  incompressible  and  compressible  flow  tend  to  zero  as  as  ui  tends  to  zero. 
He  also  notes  that  the  shape  of  the  spectrum  is  independent  of  M,  only  its  level 
is  mach  dependent.  I  do  question  the  validity  of  these  general  results  and 
there  is  no  credible  evidence  to  support  his  predictions  at  least  for  u)6/V.<l. 
His  predictions  are  at  odds  with  power  spectrum  predictions  by  others  and 
results  predicted  by  myself  in  tract  V.  I  do  have  reason  to  believe  that  $i(u>) 
and  <t>2(‘*>)  do  exhibit  an  ^  behavior  at  low  frequencies  but  not  for  the  reasons 
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stated  by  Williams.  The  leading  term  $q(u)  should  not  exhibit  an  u?  behavior 
unless  6*/ 6=0  according  to  the  convective  waveguide  model  of  tract  V. 

However,  I  can  bring  the  two  formulations  into  comptability  (to  within  a 
constant  multiplier)  by  noting  that  the  frequency  parameter  in  tract  V  can 

be  rewritten  as 

a2  *  *  ■  (¥?)2  *  (f)2  *  ‘*2(fe2u^)2ci*(S2]t  (4) 


In  this  analysis  I  will  call  the  frequency  variable 

'i— 

V  6  »  *  wo' 


<t>.  ■  <£>/£&>*  * f  C 

C  C  '  »  o 

f-c  ,  A)2  « 


1  oj6; 


a) 


b) 


(5) 


and  treat  rr-  and  6*  as  constant  parameters. 

*  oo 

If  the  inequality  in  Eq.  5b  holds  then  my  results  agree  with  William's  predic¬ 
tions  but  that  condition  does  not  hold  at  low  frequencies  which  are  of  special 
interest. 

Next,  I  will  hypothesize  the  cross  spectral  density  function  r0U,n,w)  has 
the  new  similarity  form 

i^ 


r*0U.n,u>)  =  «*(w)A[  (^)^]B0C 


(6) 


where  4>*(u>)  =  r*0{o,o,w)/A(o)B(o)  is  a  known  function  and  compatible  with  Eq. 
24  in  tract  V. 

We  now  can  compute  the  new  P0(  k i » k2 a  Pure^y  formal  way  and  the 
result  can  be  expressed  as 

V<kl-k2>“>  =  A0*U-(^).ki]B0*C(^).k2]  (7) 

where 

ao 

A*(tJ  »  4  /  A(a)e‘Tlc4dcs  (a) 

.00 
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CD 

B*(t2}  =“2il  J  BCiO^dU  .  (b) 

We  now  can  follow  Williams  in  detail  and  re-express  Eq.  la)  in  the  form 

P*(kj,k2.»)  -  !>ovVgo0A*C  ]B*[  ] 


x(ao(^).2|t2  *  alM,2  *  a2H.4‘"<i)“^(y>.2  -  M,2]) 

where  the  two  mach  dependent  terms  now  are  frequency  de[ 
relation 


(9) 

according  to  the 


M*2  -  2  • 


rZ'u>  '*  ¥»  o 

i  — \c 


V»  . 

l+(— )' 

Ud 


:“>V 

(— )2 
*  tu6 

Vc2 

(JLv2 
»  'o>6' 

c2 

and  k2  =  kj2  +k£2  in  Eq.  9. 

The  location  of  the  singularity  at  k*w/c  remains  unchanged  by  this  trans¬ 
formation,  ie  by  utilizing  Eq.  5.  At  this  point  we  do  not  have  functional  forms 
for  -  cher  A*£  ]  and  B*[  ]  in  Eq.  9  that  "agree  with  experiment"  so  we  can 
test  only  specific  integrals  in  Eq.  9. 

To  complete  the  comparisons  and  highlight  differences  with  William's 
results  I  will  now  recompute  the  <j»( u> )  and  in  particular  the  $0(w)  term. 

Thus,  by  specific  integrals  v(u)  now  becomes 

•„*(“)  =  +  W*  +  ^4*n(|)}  (a)  (li) 


vVc'*  o'V*' 


w  i 

9  («)  3  ♦0*(u)  +  f^*(u))  +  <#2*(w) 
where  the  integration  constant 

*  =  a  Jj  A*[  ]B*[  ]k2(-^-)2  d2[k(-^)  ]. 


u)  *  u>  * 


(a)  (11) 

(b) 

(a)  (12) 


^  h.*'*  ^ i  A  '  ’  *_*^I fcV  - 
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i  =  aj  jj  A[  ]B  [  ]d\  k^)^].  (b) 

and  *  * 

■y  =  a2it  A  [  1  -  M*]B  (o)  (c) 

The  low  frequency  behavior  of  Eq.  11  now  should  be  compared  with  the  low 
frequency  behavior  predicted  by  Eq.  24  in  tract  V.  The  normalized  spectrum 


corresponding  to  <$Q  (w)  reduces  to 


V*<->  - 


S„(^)  .  <3 2  »  1 


a*  o'  V 


A_  2/u6>2r  ,u)6.  A 2  , 

6*  <V  >  (vj  Go(V  }  *  u)6  <<  1 


The  normalized  spectrum  now  approaches  a  constant  as  o>  tends  to  zero  which 
is  the  desired  result  for  the  leading  4>0*(id)  term  and  is  now  compatible  with 
Eq.  24  in  tract  V.  We  could  factor  out  the  term  (u>2+u^) from  Eq.  24 
in  tract  V  and  identify  the  G0(yLa)  factor  more  explicitly  but  that  is  not  neces¬ 
sary  to  make  our  point. 

However,  we  also  can  say  that  the  leading  term  in  Eq.  9 

fV<kl-k2-“>  *  VoWgoO  A't‘-kl  (14> 

oo 

could  be  used  for  an  initial  estimate  of  the  low  wavenumber  spectrum  at  low  mach 

numbers  and  low  frequencies  provided  a0.a^,A*£  ]  and  B*[  ]  can  be  determined 

from  experiments.  Certainly  this  will  be  the  case  for  k^  >>  — k/Vc2.  A  sinyu- 

=  M  d° 

larity  can  be  expected  at  but  th<lt  Point  probably  will  fall  out¬ 

side  the  range  of  our  experiments  at  low  mach  numbers. 

There  are  some  difficulties  with  making  numerical  estimates  of  P*(ki»k2*u)) 
over  a  wide  range  of  wavenumbers.  For  example,  if  we  choose 

<(“-)  d  *  e-“#'vc)*l<l  (a) 


B-23 
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BC(“-),n]  •  e-“2<t»|n| 


(») 


where  aj  and  are  empirical  constants  and  proceed  to  calculate  A*[  ]  and  B*[  ] 

from  Eq.  8  then  there  are  difficulties  at  high  wavenumbers.  The  functions  in 
Eq.  15  yield,  for  axial  wavenumbers  kj  and  lateral  wave  numbers  k£  *  0,  the  wave¬ 
number/frequency  spectrum  (see  Appendix  A) 


P  *(k.,o,«)  *  a 


UmX  ,21 


“  “0  „2  1 


JA — 


°2'Ckl-(^).]S202* 


*0^  1^12  f°r  kl  “ 


a  (■^i)  (-£)2  for  k  »  (-J-) 

0  n2  'a'  vu»  \  1 


(a)  (16) 

(b) 

(c) 


where  G(«)  »  P0V3o3Gq(^). 

CD 

Note,  Eq.  i6c  does  not  exhibit  the  proper  high  wavenumber  behavior  and  this 
difficulty  has  been  pointed  out  by  Chase  who  also  attempted  to  modify  the  Corcos 
forum  to  fit  his  theoretical  results. *  Eq.  16b,  however,  does  appear  to  exhibit 
the  proper  low  wavenumber  behavior  as  observed  by  Chase  and  others. 

On  the  other  hand,  if  one  selects 


and 


U> 

K(p,a  -  CM®, 

BC(f-)  1]  *  e  -«2(^),l"l 

c  * 


(a)  (17) 


(b) 


*The  spectrum  of  Eq.  16a  does  peak  at  k^  =  (~)^  but  the  high  wavenumber 

asymptote  is  only  14dB  below  the  peak  value  for  =  0.2.  It  should  be  at 
least  20dB  below  the  peak. 
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then  for  axial  wavenumbers  one  obtains  the  spectrum 


Po*(klto,co) 


ag^La. _  V* 

“2 

1  *c  »c 


2  2 

(flj  --l_k-L-  (\2  fop  k  «  (ji_j 
o  ,2  a2  (l+a^)  }*  *  *1  V* 


(b)  (18) 


3 

2ao^  - for  ki »  (t>. 


which  has  the  proper  wavenumber  behavior  at  high  wavenumbers  but  a  somewhat  dif¬ 


ferent  parametric  dependence  at  low  wavenumbers  than  predicted  by  Eq.  16b. 


As  of  this  date  I  have  not  made  any  numerical  comparisons  between  the 


results  of  either  Eq.  16  or  Eq.  18  and  Chase's  predictions  of  P(Tc,u>)  so  will  not 


comment  further  on  that  matter.  I  plan  to  do  so  in  the  near  future. 


However,  several  other  comments  regarding  William's  predictions  of  «(w) 


with  respect  to  the  compressible  flow  terms  and  can  be  made.  In 

the  ♦!*(<*»)  contribution,  the  tiM*2  factor  now  reveals  an  ^  dependence  for 


(y-)  <  1.  This  is  easily  seen  from  Eq.  10.  For  (— )  <<  1  the  contribution 


reduces  to  Williams  prediction.  Hence,  at  low  frequencies,  (-fr  >>  1,  the  *.(u>) 

U)0  1 


will  be  greatly  diminished  in  comparison  with  4>0(u>)  and  probably  not  of  experi¬ 


mental  significance  especially  at  low  mach  numbers,  i.e.  M  <<  1. 

In  the  $2*(w)  contribution  the  factor  now  reveals  an  J *  depen¬ 


dence  for  (y~)  <.  1.  Hence  $2*{w)  should  also  yield  an  extremely  small 


contribution  to  $*(w)  at  all  frequencies  at  low  mach  numbers  for  realistic 


values  of  R. 


V rx  -J’JV’ 


Conclusions 


Thus,  I  conclude  that  the  two  terms  in  corresponding  to  and 

are  not  experimentally  significant  on  the  basis  of  a  convective  waveguide 
model  at  all  frequencies  at  low  mach  numbers.  This  conclusion  satisfies  my 
intuition  and  agrees  with  all  experimental  results  that  I  have  seen  to  date. 

The  *on(w)  term  as  predicted  by  Eq.  24  in  tract  V,  therefore,  should  yield  good 
agreement  with  experiment  for  all  J$g>,^an d  all  frequencies  a  >  0  at  low  mach 
numbers.  Eq.  14  is  promising  as  an  initial  estimate  of  the  low  wavenumber/ 
frequency  spectrum  of  wall  pressure  flutuations  at  all  frequencies  at  low  mach 
numbers  on  the  basis  of  the  convective  waveguide  model.  Synthesis  of  the  proper 
functional  forms  for  AQ*[  ]  and  BQ*[  ]  still  await  experimental  results.  I 
intend  to  compare  numerical  predictions  of  P*(ki»k2.<*>)  based  on  Eqs.  16  and  18, 
with  predictions  based  on  the  Chase  model  for  a  range  of  parameters. 
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Appendix  A 


Some  numerical  comparisons  of  low  wavenumber  spectra 


For  this  purpose  we  must  represent  the  analytical  form  of  both  the  longitu¬ 


dinal  and  lateral  coherence  functions  Aq  and  B0.  On  the  basis  of  experiment  and 


Chases'  work  I  have  selected  them  in  the  convenient  forms 


AC(f-)^]  -  Cl  ♦  I^^Vlle  1  Vc 


(A-l) 


SC(^)*n]  -  e  Vc  * 


where  the  parameter 


U1 

°1  "r  *  r  >  0 


(A-2) 


and  r  is  a  mixture  coefficient.  If  the  pressure  fluctuations  are  dominated  by 


the  mean  shear  then  r  <<  1  and  Both  and  must  found  experiment¬ 


ally. 


An  expression  for  PQ* now  is  found  via  Eqs.  8  and  14  to  be 


Po*(k1,k2,u.)  =  ^  G(u»)(k12  +  k22)i2 


(““)  a  (a  +  a  )  ♦  [k  -  (^)  ]  (a  -  o  ) 

V..  *  i  i  i  i  V-  *  i  i 


ir  i  /—  '  t2  .  2 /  03  .  2  i  2r  r,  2  .  2 1  uj  ,  2-i 

tCki  -  (v  )#]  +  ai  (v  )if  l^k2  a2  Vr  * 

c  c  t  ^ 


(A-3) 


where 
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In  the  second  case  (case  B) 


12 


3  0  or  r  >  - 


Then  the  normalized  spectrum  corresponding  to  =  0  becomes 


2  2 


L 


PoB^1)  (x-l)2^2 
Final  1; 

poC^x^  “l 

PqC^1^  =  (x-1)2^2 


(A-10) 


Finally,  the  Corcos  spectrum  (case  C)  yields 
2 


( A- 1 1 ) 


Note,  only  one  parameter  a1  is  required  to  plot  this  ratio  (in  all  cases) 
versus  the  dimensionless  variable  x.  This  makes  the  approach  adopted  here  very 
appealing. 

To  plot  these  results  I  initially  will  choose  =  0.20,  a  nominal  value 
based  on  wind  tunnel  experiments.  In  this  case  the  Corcos  spectrum  levels  off  at 
about  -14dB  for  x<<1.  The  three  cases  are  plotted  in  the  attached  figure  for 
x<l.  There  are  substantial  differences  in  the  three  curves  for  all  x  below 
about  0.8.  There  will  be  some  differences  for  x>l  but  that  matter  is  not  of 
immediate  concern. 

There  are,  of  course,  many  other  possibilities  for  choices  of  A  and  B.  I 


deem  it  wise  to  wait  until  measured  spectra  are  available  and  then  base  any 
further  analytical  work  on  those  results.  A  numerical  comparison  with  the  Chase 
spectra  will  be  the  next  endeavor. 


VII  Spectral  Properties  of  the  Convective  Waveguide  Model 
of  Wall  Pressure  Fluctuations 

Robert  F.  Lambert 


Introduction 

The  purpose  of  this  tract  is  to  make  some  numerical  predictions  of  the 

wavenumber/ frequency  spectrum  level  employing  the  convective  waveguide  model  and 

compare  the  results  with  the  Chase  model.  Because  of  the  dispersive  froperties 

of  the  waveguide  model  at  low  frequencies  some  significant  differences  are 

expected.  In  particular,  I  woulc  like  to  compare  the  effects  of  frequency  on 

the  location  of  the  convective  peak  which  occurs  at  the  longitudinal  wavenumber 

uj/v  in  both  the  Chase  model  and  other. high  frequency  models.  It  turns  out  that 
0 

2  2 

at  sufficiently  high  frequencies,  i.e.  id  /u>q  »  1,  that  all  models  agree  where 
id  s  Vc/6  in  the  convective  waveguide  model.  All  other  frequencies  will  be  con¬ 
sidered  as  low  frequencies. 


Analysis 


In  order  to  make  a  fair  comparison  with  the  Chase  model  I  have  selected  the 


logitudinal  and  lateral  coherence  functions  to  be  of  the  special  forms 


respectively,  where  the  decay  parameters  >  o  and  >  o  are  treated  as 
constants  and  must  be  found  from  experiments.  From  tract  VI,  ( vc )*  s 


The  parameters  and  o 2  are  associated  with  a^  and  a^  by  the  formulas 


2 


Oj  =  c»j/(l+r) 

°2  *  c2^1_r)/(1+r) 


22  3  22 

°1  °1 


(1-r) 


where  r  is  a  mixture  coefficient  that  appears  to  lie  between  0  and  1.  If  r  «  1 
the  pressure  fluctuations  are  dominated  by  the  mean  shear. 

In  order  to  determine  a  self  consistent  set  of  numerical  values  for  0j,  <*£ 
and  r  I  employed  the  coherence  data  of  M.  K.  Bull  and  attempted  a  good  fit  of 
Eqs.  la  and  lb  to  his  data  for  several  choices  of  5/6  and  n/6. 

I  was  not  able  to  obtain  a  "best"  fit  in  any  mathematical  sense  but  did 
establish  by  trial  that  the  choices  a1  *  0.2,  s  1.0  and  r  =  0.6  provided  an 
acceptable  fit  to  Bull's  data.  The  resulting  ratio  a2/0j  =  5  seems  to  be  con¬ 
sistent  with  choices  made  by  others  including  Chase  who  selected  a^/a^  =  7.  For 
this  set  of  parameters  the  values  of  ACO,  5/6]  and  B[0,n/6]  are  both  less  than 
unity  for  all  values  of  Vo  and  n/o  greater  than  zero.  For  the  choices  made 

both  A  ^  and  B  ^  ,  -0  peak  at  the  origin. 

_c  J  ]_c  6J 

It  is  interesting  to  note  that  the  coherence  functions  in  Eq.  1  may  or  may 
not  exhibt  a  maximum  value  for  ^  and  un/V  greater  than  zero  depending  on  the 

V  _  C  7""  - 1 


the  origin. 


relative  values  of  the  other  parameters  involved.  For  example,  if  A 


X  i( 

LVc  ’  *J 


plotted  versus  7T  employing  5/0  as  a  parameter  it  will  exhibit  a  maximum  for 
vc 

|w5/Vc |  >  0,  provided 

aj  <  Oj  a) 


«1  ox 


>2  >  <^)2<t)2 


For  M.  K.  Bull's  data  that  did  not  appear  to  be  the  case.  Moreover,  Eq.  2 


b  -32 


requires  that  o,  <  for  r  >  o. 

From  the  definitions  in  tract  VI  we  now  have  by  direct  inteyration  of  Eq.  1 


The  above  expressions  are  then  incorporated  into  the  leading  term  for  the 
wavenumber/ frequency  spectrum  (Eq.  9  of  tract  VI)  namely 


*  *  Vr  k  p 

Po(kl’  k2*  =  aoG(u)A  B  ■'  „  ! )*  a) 

where 

s(»)  ■  P0VC3  s3G0(ff)  b)  (5) 

and 

k2  -  kx2  +  k22  c) 

This  formulation  yields  a  power  spectrum  due  to  incompressible  fluid  flow  of 
the  form 

*0(“)  =  a(v“)*2G(w)  (6) 


where  a  is  a  constant  of  inteyration  and  in  the  notation  of  tract  V  the  spectrum 
parameter 
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Go0 


i [0  0 
0 


.  2 
“b 


(f)404  *  (S2-2ac2)'f)2(^)S 


(7) 


where  the  parameters  are  related  to  the  power  spectrum  center  frequency  o)c  and 
the  power  spectrum  bandwidth  via  definitions  in  tract  V. 

By  setting  1^  in  Eq.  5  to  zero  the  longitudinal  wavenumber/ frequency 
spectrum  reduces  to 


P0( ki *  °*  =  w2 


kl^.  2  ag+q, 

v  u  /★  '  2  ' 


(8) 


al  (al+ai)  + 


,klvc 

(„)*  -1 


{Oj  ♦ 


<„>.  -i 


2,2 


Note,  aQ  is  an  un<nown  paramter  in  this  formulation  and  also  both  and 


are  involved  in  the  prediction. 


Vcki 

The  peak  occurs  at  =  1  which  yields 


U> 


^  2, 


(y  )*#  °.  «*» 
c 


— .  d  (  )*  G(u>) 

0  10  (  ) 


2  2 


(9) 


U1  a2 


which  ould  afford  an  evaluation  of  aQ  if  the  other  parameters  are  known  from 


independent  measurements. 

Since  at  present  we  don't  have  values  for  all  the  parameters  it  is  expedient 
and  convenient  for  comparison  purposes  to  normalize  the  spectrum  with  respect  to 
its  peak  value,  Eq.  9.  After  normalization  the  longitudinal  wavenumber/fre¬ 
quency  spectrum  level  takes  the  form 
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10  log1Q 


P(x*) 


P(l) 


10  log10 


■,4  * 


-here  x,  >  (k,*c/«).  -  <  „  .  2- 


L  [a,2  ♦  (x„-l)2]2  J 


1 


x*2}dB 


(10) 


Numerical  Predictions 


y 

The  above  expression  can  be  plotted  either  as  a  function  of  x*(or  x  *  if 
“o/u  is  taken  as  a  parameter)  for  both  positive  and  negative  values  of  x*.  It 
will  peak  at  x*=l  and  approaches  a  constant  as  x*-*».  The  value  of  this  asymp¬ 
tote  is  important  and  is  quite  sensitive  to  the  difference  aj-Oj.  If  a^Oj, 
i.e.  r=o,  then  P(x*)+0  as  x*+®. 

Typical  plots  of  Eq.  10  for  the  parameters  a,a0.2  and  r=0.6  are  shown  in 
Figs.  1  and  2.  Note,  in  particular  tnat  Eq.  10  levels  off  at  about  -20  dB  for 
both  x*  >>  1  and  x*  «  1  for  this  choice  of  parameters. 

The  same  general  type  of  plot  can  be  made  for  the  lateral  spectrum 
P(0,  k2«  <**)• 

After  normalization  the  lateral  wavenumber/frequency  spectrum  level  takes 
the  form 


10  log1Q 


£1x1 


P(l) 


=  10  log^g 


4  +  2 

al  al  ai+0\ 


L  (a,2  +  l)2  J 


*2  *  *2 


L  (a  2  ♦  y*2)2 


yj}d& 


xk?VC 

where  y*  =  (-£“*), 


(ID 


Eq.  11  is  symmetrical  in  y*  and  a  typical  plot  is  shown  in  Fig.  3.  Note 

again,  the  function  tends  to  level  off  at  about-20  dB  for  y*  >>  1.  These 

results  suggest  that  I  have  selected  a  self  consistent  set  of  parameters  c^, 

a2,  o2  on  the  basis  of  Bull's  measurements. 

At  this  point  I  am  satisfied  that  a  meaningful  comparison  between  the  level 

predictions  of  Eq.  10  and  Chases'  predictions  can  be  made.  Chase  also  selected 

the  parameters  in  his  numerical  work  from  Bull's  data.  Several  comparisons  are 

2  2 

shown  in  the  next  set  of  figures.  It  should  be  obvious  that  for  “*o  /w  «  1, 

kiVr  k1Vf  w0 

(')*  ~  and  the  comparison  is  direct.  On  the  other  hand  for  >1  that 

will  not  be  the  case  and  effects  of  waveguide  or  other  frequency  dispersion 

should  be  evident. 

A  plot  of  the  Chase  spectrum  is  shown  in  Fig.  4  employing  frequency  as  a 
parameter.  In  all  of  the  comparisons  o=lcm,  V^lOcm/s  and  x=k^vc/uj.  In  the 
Chase  model  all  of  the  curves  peak  at  x=l  independent  of  id  but  there  are  some 
dispersive  effects  at  frequencies  below  500  Hz  which  are  very  evident.  The  nor¬ 
malized  level  increases  with  decreasing  frequency  quite  significantly  at  low 
wavenumbers.  Note,  above  1000  Hz  there  is  no  significant  change  in  level  with 
increasing  frequency  at  all  wavenumbers  for  the  parameters  chosen. 

kiVc 

Plots  of  Eq.  10  versus  x=  on  the  other  hand,  exhibit  a  significantly 

different  behavior  if  frequency  is  employed  as  a  parameter.  This  is  quite  evi¬ 
dent  in  Fig.  5  where  uio=v*/o=1000  and  u/u  is  employed  as  the  parameter.  The 

function  peaks  at  k  =  at  high  frequencies,  i.e.  “o^/id^  <<  1,  as  expected  but 

1  c 

the  peaks  shift  to  higher  values  of  k^  as  u  is  decreased  especially  in  the  range 
u)  s  <dq.  Note,  the  predicted  level  decreases  at  low  wavenumbers  as  frequency  is 
decreased  on  this  scale.  Fig.  6  shows  the  same  results  but  on  a  logarithmic 


wavenumber  scale  which  is  easier  to  use  at  low  wavenumbers  for  comparison  pur¬ 
poses. 

It  also  is  informative  to  plot  the  Chase  spectrum  level  versus  x*=(kivc/w)* 
for  the  parameter  wo=1000  and  then  employ  frequency  as  a  parameter.  This 
characteristic  is  shown  in  Fig.  7.  Note  the  shift  in  wavenumber  peak  is  now 
downward  as  frequency  is  decreased  and  there  exists  a  correspondingly  greater 
increase  in  normalized  level  with  decreasing  frequency.  This  general  result  is 
to  be  expected  and  such  trends  should  be  observable  experimentally. 

The  last  two  figures  show  the  effects  of  varying  about  a  nominal  value  of 
0.2.  In  the  results  shown  in  Fig.  8  r=0.6  is  fixed.  Thus,  one  can  see  that 
experimental  errors  in  finding  and  can  be  significant,  the  greater  the 
decay  rates  the  higher  the  predicted  level  at  both  low  and  high  wavenumbers. 

Fig.  9  shows  the  extreme  case  where  r=o.  Then  and  the  levels  are 

greatly  diminished  at  both  low  and  high  wavenumbers  for  the  same  choice  of  c^. 
The  corresponding  wavenumber/ frequency  spectrum  is  very  sharp.  Hopefully,  such 
trends  also  will  be  observable  experimentally. 
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Figure  Legend 

I.  Normalized  longitudinal  wavenumber/frequency  spectrum  level  versus 

ki  V 

^or  *+  >  0  and  °i=0*20,  r=0.6.  Convective  waveguide  model. 

II.  Normalized  longitudinal  wavenumber/frequency  spectrum  level  versus  x*  for 


x*  <  o  and  0^=0. 20  and  r=0.6.  Convective  waveguide  model. 


III.  Normalized  lateral  wavenumber/ frequency  spectrum  level  versus  y*2  (  L 

"  (jJ  '  " 

for  aj=0.20,  a=1.00,  and  r=0.6.  Convective  waveguide  model. 

ki  Vc 

IV.  Normalized  longitudinal  waveguide/frequency  spectrum  level  versus  x— -*•  v 

0J 

employing  frequency  as  a  parameter  for  6=lcm.  Chase  model. 

ki  V 

V.  Normalized  longitudinal  wavenumber/frequency  spectrum  level  versus  x- 

U> 

employing  “/u^  as  a  parameter.  Here  u)Q=1000r/s  corresponding  to  6=lcm  and 
V«=10  cm/s.  Linear  wavenumber  scale.  Note  the  predicted  wavenumber  shift 
due  to  waveguide  dispersion. 

VI.  Same  plot  as  figure  V  but  x  is  shown  on  a  logarithmic  scale.  Note  pre¬ 
dicted  wavenumber  shift  due  to  waveguide  dispersion. 

ki  Vc 

VII.  Chase  spectrum  when  plotted  versus  x*=(-  ^  )*.  Note,  the  wavenumber  shift 

due  to  dispersion. 

VIII.  Sensitivity  calculation  showing  the  effects  of  variations  in  uj  in  the 
longitudinal  wavenumber/frequency  spectrum.  Here  r=0.6  is  held  fixed. 

IX.  Sensitivity  calculations  showing  the  effect  of  lowering  the  mixing  para¬ 


meter  r  to  zero. 
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VIII.  Properties  of  the  Spatial  Coherence  Functions 
Wall  Pressure  Fluctuations 
by 

Robert  F.  Lambert  and  Keith  Alcock 


Introduction 

The  purpose  of  this  tract  is  to  focus  attention  on  the  functional  form  and 
evaluation  of  numerical  parameters  for  the  longitudinal  and  lateral  spatial 
coherence  functions  that  were  employed  in  tracts  VI  and  VII  to  predict  the 
spectral  characteristics  of  wall  pressure  fluctuations  in  boundary  layer  flow. 
Also  in  Appendix  B  a  new  experiment  is  suggested.  In  tracts  VI  and  VII  we 
employed  these  functions  in  the  special  forms 


of  the  model.  The  problem  Is  how  to  optimize  subject  to  constraints  the  choice 


of  aj,  Oj,  Og,  and  a In  particular,  their  values  depend  on  the  ratio  V«/Vc 
via  Eq.  2.  Initially,  It  was  assumed  that  ^  «  0.80  a  value  characteristic  of 
the  main  energy  bearing  eddies  whereas,  in  actuality,  the  ratio  is  frequency  (eddy 

III  ^  ^ 

size)  dependent.  At  extremely  low  frequencies,  (i.e.,  y“  <  0.2,)  it  is  usually 
y  oo 

assumed  that  Tp  approaches  unity.  At  extremely  high  frequencies,  i.e.,  y"  >  5, 

Vco  v» 

the  measured  value  is  closer  to  0.6.  (See  Appendix  B). 

A  method  for  measuring  Vc  and  determining  values  for  both  A  and  B  was 

described  in  tract  I.,  Eqs.  6  and  7. 

In  tract  VII  the  fit  of  Eq.  1  to  coherence  measurements  by  M.  K.  Bull 

was  reasonable  but  not  entirely  satisfactory,  especially  for  small  values  of 

and  y0,  .  It  appeared  that  one  could  obtain  a  better  fit  of  Bull's  data  if 
vc  c 

Vc/V.  *  0.60  (see  tract  V.)  The  observed  drift  velocity  is  smaller  at  small 
spatial  separations,  i.e.,  small  values  of  5/6  and  n/6. 

It  was  decided  to  relax  the  requirement  that  and  be  positive  numbers 
and  to  reformulate  the  A  and  B  functions  as 


and 


L 

<5 


1  ■  1 
c 


]  =  OojJOJle 


a) 


(3) 


b) 


and  then  to  attempt  a  new  fit  for  ^  =  0.60  utilizing  the  values  of  5/6  and  u/6 

00 

employed  by  Bull. 


Numerical  Predictions 

The  values  of  the  four  numerical  parameters  obtained  by  a  numerical  curve 
sitting  algorithm  are: 

<ij  =  0.166  =  0.786 

Oj  =  0.074  o  £  =  -.0842 
which  yields  a  mixture  coefficient  r  =  1.24. 

Plots  of  A,  Eq.  3a,  versus  u)5/V employing  ratios  of  5/5  =  0,  1.23,  and  2.52 
are  shown  in  Fig.  1.  The  data  of  M.  K.  Bull  for  5/5  =  1.23  and  2.52  are  shown 
for  comparison. 

Plots  of  B,  Eq.  3b,  versus  u>n/Vc  employing  n/5  =  0,  0.103,  and  0.627  are 
shown  in  Fig.  2.  Bull's  data  points  for  the  same  ratio  are  shown  for 
comparison. 

The  agreements  now  are  more  satisfactory  for  close  spatial  separations  but 

not  so  good  for  values  of  oin/V  >  1  as  before.  Again  maximum  values  of  both  A 

c 

and  B  occur  at  the  origin.  The  peak  value  for  5/6  *  0  and  n/6  *  0  is  unity. 
Several  characteristics  of  the  predicted  cross  spectrum  AxB,  where 


am 

V 

c 


00 


are  shown  in  Figs.  3,  4,  5  and  6  for  several  parametric  choices  of  5/S  and  n/5. 
Note,  the  surfaces  are  similar  {but  not  identical)  in  shape  to  cross  spectrum 
measurements  obtained  in  the  water  tunnel,  save  for  differences  in  amplitude  and 
spatial  scales  (See  Fig.  9  in  the  April  1984  -  October  1984  program  progress 
report. ) 

In  the  last  three  figures  (7,  8  and  9)  are  plots  of  normalized  wavenumber/  . 
frequency  spectra  P(kj,0,u>)  versus  (k^,Vc/uj)*  for  a  range  of  both  positive  and 
negative  kj  and  PCO^.a)  versus  (k^y^/u)*  for  a  range  of  positive  ( k2 Vc/u> ) * . 
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These  curves  are  to  be  compared  with  those  in  Figs.  1,  2,  and  3  of  tract  VII. 

As  can  be  seen  the  two  sets  of  predicted  spectral  characteristics  are  for  all 
practial  purposes  the  same.  This  is  evident  by  noting  that  the  dominant  term  in 


Eq.  10,  tract  VII,  at  both  low  and  high  wavenumbers  is  proporational  to 
2  ai”°i 

a\  (a  +c  )»  a  factor  that  did  not  change  by  more  than  half  a  decibel.  The 

1  1  a2‘°2 
same  result  is  obtained  for  Eq.  11  where  the  a  2(— -)  term  changed  by  at  most 

(.  a2  °2 

1  decibel.  The  conclusion  is  that  any  changes  in  the  cross  spectrum  AxB  would 


have  to  be  quite  drastic  to  change  the  functions  by  5  decibels  or 

more  at  both  low  and  high  wavenumbers.  Even  the  peak  value  given  by  Eq.  9, 
tract  VII,  did  not  change  by  more  than  0.5  decibels  in  this  example.  It  should 


be  possible  to  improve  the  level  predictions  by  assuming  an  appropriate  fre¬ 


quency  dependence  for  Vc/V*  and  repeating  the  evaluation  of  oj,  o1#  a2  and  o2. 
But  that  alone  is  not  likely  to  change  the  spectral  level  of  P(kj,  k2,  <*>)  by  a 
significant  amount.  However,  it  will  change  the  predicted  location  of  the  con¬ 


vective  peak  in  P(klt  o,  w)  (see  Appendix  B). 
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Appendix  A  -  Curve  fitting  procedures 


In  this  application  the  problem  Is  to  convert  two  sets  of  data  points 


(xll *yll ) *  a^d  ^x21,^21^>***  *^x2n*^2n^  are 

characterized  by  parameters  k  =  alfa2,  respectively,  into  a  specified  function 


Aw 


f(x,k)  =  [l+o  /xSk"  3e"WA  T*  (A_1) 

with  a  and  o  selected  to  best  fit  the  data  points.  The  approach  used  is  a 


variation  of  the  method  of  least  squares  by  Gauss. 


The  first  step  is  to  construct,  by  hand,  smooth  curves  through  each  set  of 


data  points.  These  curves  can  correct  any  data  errors  or  deviations  that  do  not 


appear  to  be  representative  of  the  entire  data  set. 


The  initial  curves  so  contributed  are  then  sampled  at  intermediate  point 


into  two  new  sets  of  data  points  (*n*.yu*) ,•••  »(xlm*.ylm*)  and  •”* • 

if  it 

(x2n  ,y2n  )  which  then  replace  the  initial  sets.  Such  samples  can  be  taken 


either  evenly  spaced  or  weighted  (i.e.,  close  together  where  an  accurate  fit  is 


very  important)  and  thereby  provide  discrete  data  to  be  entered  into  a  computer 


routine  and  then  compared  to  f(x,k)  for  choices  of  a  and  a. 


The  next  step  is  to  write  a  program  that  will  generate  values  of  a  and  o  and 


calculate  f(x,k)  using  the  values  *xln*  and  x21*’x22*’’ *  * ,x2n*  as 


generated  above  for  domain  values  and  k  =  a^^  for  appropriate  curves.  The 


difference  in  the  range  values  yn  »  y^2  »*”,ylr 


y21  *y22  ****  ,y2n 


the  values  of  f(x,k)  is  a  measure  of  the  accuracy  of  the  value  of  a  and  o.  For 


each  o  and  a  all  values  of  f(x,k)  are  calculated  for  each  domain  value, 


subtracted  from  the  range  value  (of  the  sampled  point),  squared,  and  summed. 


That  is 


in  «  n  a 

VAR  Cf(xli*,k1)-yli*r  ♦  I  [f(x2i  ,k2)-y2i*r  (A-2) 


B  -52 


6 


which  Is  analogous  to  the  variance  or  area  between  f(x,k)  and  the  data  points. 

The  problem  now  reduces  to  finding  values  of  a  and  o  that  minimize  VAR. 

Best  fit  values,  a.*  and  o.  are  recorded  and  the  function  f(x,k)  then  is 
of  Dt 

plotted  for  continuous  x  as  a  comparison.  If  the  fit  to  the  data  point  is  not 
satisfactory  the  original  curve  can  be  re-sampled  with  greater  weighting  at  the 
segment  of  greatest  variance  and  the  steps  repeated.  One  could  also  conclude 
that  the  function  f(x,k)  is  not  appropriate  for  the  data  set,  which  does  not 
appear  to  be  the  case. 


A-VALUE 


W*Zato/Vo 


B-VALUE 


u 


a> 

x 

< 


(kj  Vc/w) 

Fi9-  7 


Appendix  ^-Suggested  Experiment 


It  was  pointed  out  in  tracts  VI  and  VII  that  the  convective  waveguide 
model  of  boundary  layer  pressure  fluctuations  predicted  a  frequency  dispersion 
not  found  in  the  Chase  model. 

The  wave  guide  model  predicts  that  P(lq,  o,  o>)  will  exhibit  a  peak  at  the 
wavenumber  value  kp  given  by 


whereas  the  Chase  model  predicts  a  peak  et  w/Vc.  If  at  the  frequency  in 
question,  u>»u>0  *  V.,,/5,  then  the  two  predictions  are  about  the  same. 

However,  at  low  frequencies,  i.e.  u<<j0,  a  significant  difference  is 
expected  and  this  suggests  that  we  try  to  verify  this  prediction  by  experiment. 

Note  that 

<B-2> 

and 


=  1  /Jfax 

5  W  ' 
c 


>♦  &)2<f¥ 


( B— 3 ) 


Thus,  in  the  convective  waveguide  model  the  peak  in  P(k^,  o,  u>)  should  occur 
at 


The  ratio  V«/Vc  is  not  constant  but  should  vary  with  frequency  according  to 
the  empirical  relation 


1+0. 67e 


-MWJ 


(B-5) 


which  predicts  that  V  *V.  at  wSVV^O  and  V  +0.60  V*  as  U6+/V*-*—  in  good 

v  v* 

agreement  with  results  from  wind  tunnel  studies. 
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On  this  basis  one  can  predict  how  kp  will  behave  with  frequency  if  6*/6  is 
known.  That  prediction  is  shown  on  the  attached  figure  for  6*85*.  The  upper 
curve  is  a  plot  of  Eqs.  4B  and  58  for  the  range  of  u)6*/Va,  of  greatest  interest. 
The  lower  curve  is  a  plot  of  a>6/Vc.  The  two  curves  merge  for  w6*/Vo#>0.4  but  do 
diverge  significantly  for  oj6*/V«<0*1  •  This  should  be  well  within  the  frequency 
range  of  the  bouyant  body  facility  if  S*=.001m. 

One  of  the  initial  experiments  with  the  bouyant  body  facility  should  be  to 
estimate  the  location  of  the  peaks  in  P ( *  o,  <*>)  as  a  function  of  frequency. 
The  spectral  estimators  should  work  very  well  in  that  region  of  k. 
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IX.  Pressure  Phase  Convection  Speed  Characteristics 
Wall  Pressure  Fluctuations 
by 

Robert  F.  Lambert 


Introduction 

Th  purpose  of  this  tract  is  to  focus  attention  on  convection  speed  charac- 
terist  cs  and  make  some  predictions  based  on  the  convective  waveguide  modol .  Ti 
particular,  there  is  great  interest  in  the  behavior  of  the  convection  speeds 
with  frequency  and  how  to  interpret  convection  speeds  based  on  phase  measure¬ 
ments..  In  my  judgement  they  may  have  been  misinterpreted  in  the  past. 

Analysis 

The  relative  phase  between  two  points  a  distance  £  apart  in  the  streamwise 
direction  based  on  the  convective  waveguide  method  is 


4>*  =  k*5 


)2  + 


n_  r/ 


U)+U>r 


(1) 


A  L  C  J 

where  k*  is  a  streamwise  wave  number,  m  is  the  frequency  in  question,  Vc  is  the 
boundary  layer  convection  speed,  A  is  the  half  width  of  the  flow  channel,  and 
the  cut-on  frequency 


O) 


o 


(2) 


In  general  Vc  and  hence  u)Q  will  be  frequency  dependent.  As  will  be  pointed 
out  Vc  cannot  be  determined  directly  from  phase  measurements  and  must  be 
inferred  from  other  data. 

For  convenience  I  will  interpret  the  boundary  layer  thickness  6  to  be  iden¬ 
tical  with  A. 


Eq.  1  can  be  rewritten  in  terms  of  normalized  parameters  and  expressed  in 
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the  dimensionless  form 


where  V.  Is  the  free  stream  or  center-line  flow  speed,  6*  is  the  displacement 
boundary  layer  thickness,  and  £5  =  ai6*/Vo.. 

To  interpret  phase  measurement  we  now  define  a  wall  pressure  phase  convec¬ 
tion  speed  as 

VP  i  t  (♦> 


which  is  a  parameter  that  can  be  determined  from  relative  pressure  phase  (or 
time  delay  Tj)  measurement  made  by  point  transducers  spaced  a  known  distance  £ 
apart  via  the  relations  or  Vp*S/Td*“£/$*. 

Eq.  4  yields  a  pressure  convection  phase  speed  ratio 


It  is  evident  from  Eq.  5  that 


Vc  5 

provided  )2(— )2  at  all  frequences  of  interest.  Hence,  phase  measure¬ 

ments  at  high  frequencies  will  yield  good  values  for  Vc.  However,  at  lower  fre- 

w  ★ 

*  C  o  6  p 

quencles  where  0  <(T7”)  (~)  this  is  no  longer  the  case  and  in  the  limit  as 

*OD 

approaches  zero 


(6) 
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^0  , 6  ,  u>5 

V  - ►  -  w“. 

’•  fl+Q  5  • 


(7) 


Thus,  Me  expect  Vp/V,  to  vanish  at  zero  frequency  and  the  slope  of  the 
characteristic  will  be  5/6*  there. 

We  further  expect  the  value  of  Vc/V.  in  Eq.  5  to  be  well  behaved  in  terms  of 
frequency  so  the  above  conclusions  should  be  valid. 

However,  to  make  numerical  predictions  of  Vp/V*  based  on  Eq.  5  I  have 
decided  to  approximate  Vc/V,  by  the  semi -empirical  relation 


V. 


l+0.67e~a 

1.67 


(8) 


which  predicts  that  as  This  characteristic  satisfies  my  intuition 

* 

about  the  drift  velocity  of  large  low  frequency  eddies  and  satisfies  experimen¬ 
tal  observations  at  high  frequencies  where  Vc“Vp,  Eq.  6.  It  is  essentially  the 
same  assumption  made  by  Chase  in  his  numerical  predictions  of  Vp/V„.(l) 


Numerical  Predictions. 

I  will  assume  for  prediction  purposes  that  6*/5=0.125  which  is  a  fairly 
nominal  value  for  this  ratio.  Employing  Eq.  8  for  Vc/V*  the  predicted  value 
for  Vp/V,*  is  shown  in  Fig.  1.  Note,  it  falls  to  zero  rapidly  as  3>o  and  peaks 
to  a  value  of  about  0.84  at  3=0.30.  For  higher  values  of  3  the  characteristic 
levels  off  at  the  values  required  by  the  !/<./!!„,  characteristic.  This  prediction 
is  in  remarkable  agreement  with  pressure  phase  speed  values  measured  by  Blake 
(2).  However,  his  measured  values  do  depend  somewhat  on  the  relative  spacing 
C/6*  used  In  the  measurements.  But,  like  all  phase  (or  time  delay)  measurements 
employing  random  signals  and  finite  size  transducers  there  are  some  errors  to  be 
expected.  Also,  I  believe  that  the  observed  behavior  of  Vp/V*  with  C/6*, 


-*  ,s  £-  ->0. 


¥  14 
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Including  measurements  reported  by  Aupperle  and  Lambert  (3),  has  as  much  to  do 
with  the  finite  size  of  the  transducers  at  close  spacings  as  with  the  physics  of 
the  flow  field.  Such  measurements  are  not  highly  repeatable  and  have  never  made 
good  sense  to  me.  The  results  predicted  by  Eq.  5  are  independent  of  spacing. 

So,  In  practice  we  actually  measure  Vp  and  the  value  of  Vc  must  be  inferred 
from  Eq.  5  taking  into  account  the  expected  frequency  dispersion  at  low  frequen¬ 
cies.  A  simple  manipulation  yields 
Vc  _ a _ 

*  QO 

as  an  algorithm  for  computing  Vc/V,,  from  measured  values  of  Vp/V*  for  known 
values  of  a  and  6*/6. 

Finally,  a  few  remarks  about  the  fluctuating  pressure  group  velocity  Vg, 
which  usually  is  defined  as 

<io> 

In  most  dispersive  media  (where  energy  losses  can  be  neglected)  Vg  is  a 
measure  of  the  speed  at  which  the  energy  associated  with  the  pressure  fluc¬ 
tuations  at  the  wall  propagated.  It  also  is  frequency  dependent. 


The  convective  waveguide  model  predicts  that 


The  values  of  Vg/V*  predicted  by  Eq.  11  will  be  greater  than  the  value  Vc/Vw 
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at  all  frequencies.  In  particular,  the  group  velocity  can  exceed  the  speed  of 
sound  c  below  some  frequency  fls  at  all  mach  numbers.  I  expect  fis  to  be  fairly 
low  at  low  mach  numbers  and  hence  approximate  Vc  by  V„  for  prediction  purposes. 
The  resulting  expression  for  ffls  then  becAMjey 


or 


M. 


a) 


b) 


(12) 


where  M»*V"/c«l. 

The  value  of  ws  predicted  by  Eq.  12  should  mark  the  lower  edge  of  the  con¬ 
vective  energy  regime  assuming  incompressible  fluid  flow  since  Vg<c  for  0>0S  and 
Vg>c  for  o«3s. 

In  conclusion,  I  recommend  in  the  future  that  all  speed  and  spectral 
measurements  also  be  interpreted  in  the  light  of  the  convective  waveguide  model 
of  boundary  layer  pressure  fluctuations.  Of  special  note  are  low  frequencies 

o  p  p  6  p 

i.e.  <  (v )  (5“)  »  where  dispersion  is  likely  to  be  most  pronounced.  The 
characteristic  shown  in  Fig.  1  is  a  direct  result  of  such  dispersion  and  the 
assumed  innocuous  behavior  of  Vc/V„.  I  find  the  high  frequency  model  of  Core** 
and  even  the  more  elaborate  one  due  to  Chase  somewhat  misleading  and  incon¬ 
sistent  with  certain  experiments.  Fortunately,  Chase's  model  has  a  large  number 
of  adjustable  parameters  so  that  many  experimental  characteristics  can  be  fit  by 
judicious  choices  of  these  parameters.  Of  course,  if  our  spectral  measurement 
at  low  frequencies  do  not  follow  the  predicted  dispersion  requirements  of  the 
convective  waveguide  model  I  will  modify  (or  abandon)  it  quickly. 
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ERRATA  i  OK 


FINAL  TECHNICAL  REPORT 
RESEARCH  PROGRAM  IN  HYDROACOUSTICS 

p.  25,  line-  2:  instead  "of  "monogram",  should  reau  allograph"  . 

p.  26,  Eq .  5:  instead  of  "  ] e  “ ,  should  read  "  ]e  w  ". 

p.24,  Fig.  V2 . 7 :  instead  of  "Experimental  plot  of  the  location  of  the  convective 
peaks  versus  the  normalized  frequency  w  for  the  bouyant  body 
facility.  The  solid  line  is  a  fit  of  Eqs .  5  and  7  for  the 
parameters  given  in  the  legend."  should  read  "Comparison  of  the 
phase-speed  prediction,  Eq .  4,  with  experimental  data  of 
Blake . " 

p.28,  Fig.  V2 . 8 :  Instead  of  "Comparison  of  the  phase-speed  prediction,  Eq .  4, 
with  experimental  data  of  Blake."  should  read  "Experimental 
plot  of  the  location  of  the  convective  peaks  versus  the 
normalized  frequency  w  for  the  bouyant  body  facility.  The 
solid  line  is  a  fit  of  Eqs.  5  and  7  for  the  parameters  given  in 
the  legend." 


